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This report will only provide the second part of a total research which consists 
of three projects which are shown as follows: 

Project 1: Evaluation and progression of the space feeding concept from early 

missions to the present stage. 

Project 2: Projected nutritional criteria for the Space Operations Center. 
Project 3: Projected nutritional criteria for the Closed Ecological Life 
Support System. 
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Based on experience from past missions, there are known physiological as well 
as psychological stresses existing in crev/ members as a result of space flight, 

It is also known that stress can contribute to the increased less of calcium and 
nitrogen from the body. Since high density foods will be made with high calorie 
nutritionally complete liquid, it will be able to protect the body's loss of 
nitrogen (4, 5), 

After considering all of these needs for the crew in an extended mission, 
the researcher looks for these projects to develop new high density foods, 
OBJECTIVES 

a. To achieve the nutritional need in LEO Space Station. 

b. To bring varieties in the diet of crews for extended mission. 

c. To provide .highest cjlorie in concentrated form. 

d. To provide complete nutrients, 

e. To study the shelf-life of the products produc'^d, 

BACKGROUND AND REVIEW OF THE RELATED LITERATURE 

The space flight missions are always supported with nutritious foods for 
the, crews of the space flight. Research and technological developments in 
this field are going on as the man's desire to stay longer and longer in space 
is becoming increasingly inevitable. The dehydrated, frozen, irradiated and 
thermostabilized foods have played a very important role in space flights 
supplying necessary ingredients needed for supporting life in orbit and in 
space. 

In this preview in mind, the researcher is interested in the preparation of 
high density food v^hich will meet the criteria of space food. Following facto's 
should be reviewed before applying any theoretical concepts foods into practice, 

1. Nutritional Va^ ue. The Nutritional value of the food should be formulated 
to meet the need o"^' the human in a relatively simple product. Nutritional 
availability must te fully supported by the food in quality and quantity. 


2. Safety. The food intake should be safe in a given environment, ^ood consumed 
in a closed environment should take into account the special needs, if any, for 
safety as to its side effects, especially in regard to microbiological growth 
as well as its toxicity levels. 

Convenience . Food items should be convenient to use by the crews with con- 
siderations &s to their age, sex, health status and also to the personal 
preferences of the individual. 

4. Food Acceptance . Very often food habits are very much of personal prefer- 
eftces given that variety of foods are available for consumption. Hence factors 
must be considered for their acceptability. 

With these objectives in view, the researcher has tried to evaluate, the 
high density liquid foods designed for enteral feeding by the Organon Food 
Company. Products evaluated were "Magnacal ," "Sumacal" and "Microl ipids. '' 

These products are used for hospital patients who may suffer from protein 
calories undernutrition, and of Cardiac Cachexia. 

Magnacal is a lactose-free formula designed for both total and supplemental 
nutritional support in those Clinical conditions that demand caloric and 
nutritional requirements which standard liquid f(;edings cannot satisfy The 
exceptionally high caloric density (2,0 cal/ml) makes Magnacal especially 
suitable for patients who cannot tolerate high volume feedings but who require 
high caloric levels. Magnacal is useful in the dietary management of patients 
with a wide variety of hypermetabol 1c conditions resulting from such causes as: 

1. Major Burns 

2. Severe Trauma 

3. Multiple Fractures 

4. Major Sepsis 

5. Major Surgery 


6. Advanced or Metatastic Cancer 

7. Severe Anorexia 

8. Hyperactive States or Prolonged Muscular Hypertonicity 

Magnacal is also effective as an oral supplement for patients requiring 
additional calories (15). 

Microlipid is useful in the dietary management of patients requiring 
caloric supplementation. Microlipids high concentration of calories (4,5 cal/ml) 
is derived solely from fat and can provide energy for metabolism without an 
undue increase in the volume of the diet. Microlipids are rich in linoletic 
and other unsaturated fatty acids. Microlipids may be useful in the following 
conditions (16): 

1. Fatty Acid Deficiencies 

2. Uremia 

3. Diseased and Debilitated States 

4. Postoperative Convalescence 

5. Preoperative Preparation of Malnourished Patients 

6. Severe Burns 

7. Major Sepsis 

Sumacal is designed for nutritional support in those clinical conditions 
that demand a Supplemental Source of carbohydrate calories Unflavored Sumacal 
can provide additional carbohydrate calories in combination with other Organon 
nutritional products nutritionally complete formulas. Sumacal contains no fat 
or protein which makes it especially useful as a supplement tor patients on fat 
and protein restricted diets, Sumacal may be beneficial in :ne dietary manacement 
of patients with the following condition (23): 



1. Severe Burns 

2. Major Sepsis 

3. Severe Trauma 

4. Acute or Chronic Renal Disease 

5. Major Surgery 

6. Cystic Fibrosis 

7. Anorexia 

8. Pre- ind Postoperative Care 

Dr. Steven B. Hemsfilsd in his paper on Cardiac Cachexia has pointed out 
some of the uses of these HDL products (22). 

The enteral route of refeeding is selected first unless severe intestinal 
problems limit intake. If the full dose of the prescribed diet cannot be 
orally ingested, then the balance is administered through a fine-bore naso- 
enteral cannula. The enteral formula selected for tube feeding should have 
the following properties; high caloric density (2,0 preferred over 10 Kcal/cc), 
low sodium (<2 g/day), and Kcal/gmN ratio of about 150:1. A commercially pre- 
pared solution is available which meets these criteria (i,e., MAGNACAL^^\ 

Organon Nutritional Products), or the formula can be prepared by the dietitian 
from readily available modules . If only part of the daily requirements can 
be met by enteral feedings, the balance of daily requirements are supplied by 
peripheral vein fluids (i,e,, P-900^, Intralipid^ etc. Recall that the 
undernourished cardiac patient is unduly susceptible to repletion heart 
failure, and therefore requires careful clinical monitoring and gradual 
increasing doses of nutrients. Dr, Horwitz, et, al , observed the use of the?e 
HDL foods in cases of Ascitises, Cirrhosis and Cachexiac 

ACC patients with minimal elevations in serum bilirubin and plasma ammonia, 
often resistant to standard enteral formulae, can be effectively repleted using 
modular or fixed composition high density solutions (11), 


Enteral and perenteral hypersli mentation can renouri*‘h the severely wasted 
malabsorbing diabetic as shown in a study conducted by Dr. T. Hersh. et al (10), 

Lean body mass can now be maintained or repleted by the enteral route in 
contrast to the parenteral route in patients with a wide variety of gastro- 
intestinal and other diseases. The applicability of enteral hyperalimentation 
as a primary form of therapy in inflammatory bowel diseases, in reversing the 
cachexia of advanced cirrhosis, in protecting irradiated cancer patients from 
onteric and homotologic toxicity and in the treatment of acute renal failure, 
hold promises but must await the outcome of carefully controlled clinical 
studies (21), 

George L. Blackburn, et al,, have pointed out in their study about the 
nutritional and metabolic assessment of the hospitalized patients that the 
nutritional and metabolic profile indicates the causes of malnutrition, namely 
insufficient protein and calor‘e-> and establishes individual guidelines for the 
patient's nutritional therapy (2) 

At present enteral hyperalimentation are used in hospital patients only. 

The researcher feels that these products can be used for the conventional diets 
with further study and modification to improve physical characteristics. 

The Pill'ibury Company also has made some remarkable contribution in the 
preparation of the high density bar to supplement the food supply for the 
Skylab 4 mission which was extended from the nominal mission of 56 days to a 
mission of 85 days in duration (18), 

Jack R. Durst conducted research about stable food pieces, has oointed out 
that compact food units designed to provide a completely balanced diet provided 
by incorporating essential nutrients into a unitary, flexible food piece con- 
taining a continuous external phase of a hydrophilic film former, water and an 
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edible humectant. The external phase encapsulates an internal phase of minute 
fat globules. The food pieces are resistant to physical, chemical, and bacterio- 
logical degradation (13). 

A protein fortified jod bar of a controlled calorie content is composed of 
several crisp wafers stacked one above the other with a creamy filling between 
them. Each wafer is composed of 10 parts of flour, about 6 parts of added 
protein, e.g., calcium caseinate or a lactalbumin-casein coprecipitate, about 
0.8 parts oil as a release agent and a small amount of chemical leavening. The 
filling layered between the wafers consists of 10 parts shortening, about 10 
parts of a finely divided protein such as milk protein, about 6 parts sugar, a 
minor amount of flavor if desired and a vitamin and mineral moisture, if desired, 

A confectioners coating can be applied over the composite bar if desired. 

Jay E. Morgan, in his research about high protein food bar, has discussed 
that (12). Betty L. Brooking conducted similar research as J. Durst about pro- 
cess for forming stable food pieces has found out that foods remain flexible 
during storage and do not lose oil when molded or otherwise formed into pieces (1), 

Food bar and method for making, another research conducted by Jack R, Durst, 
showed that a compact, solid food unit characterized as possessing hardness, 
frangibility and water dispensible characteristics is imparted by utilization 
of a binder which provides a structural matrix for discrete edible food particles. 
The binder is comprised of a discontinuous phase of fat globules encapsulated 
by a continuous phase of water and a hydrophilic film former (14). 

EXPERIMENTAL DESIGN AND PI^OCEDURE S 

The research activity was conducted in different phases to formulate high 
density foods. A standard recipe (#300) for Scones/Biscuits adapted from the 
^ recipe book of Culinary Arts Institute (7). 



In the experimental phase, the researchers tried to make Scones/Bi scui I' 
changing variables sucn as replaci;-q milk with various high density liquids pro- 
duced by the Organon Food Company. Five different recipes were studied. In 
each tet ipe the milk component \/as replaced by the HDL product. 

The first experimental high density Scones/Biscuits were made with 
Magnacal , a high density liquid product produced by Organon Food Comoany, 
replacing milk from the controlled recipe keeping all other ingredients intact 
as in the original recipe. 

The second high density Scones/Biscuits were made by replacing milk with 
Sumacal t all other components unchanged from the controlled recipe. 

The third high density Scones/Biscuits v;ere made, replacing milk with 
Microlipids. Since microlipid is a fat emulsion product, shortening was also 
deleted from the experimental recipe. 

The fourth high density Scones/Biscuits were made by replacing milk with 
Microlipids. 

The fifth high density Scones/Biscuits were made by replacing milk with 
'*nj^’^.al and Microlipids by doubling the controlled recipe in every respect 
except baking powder. 

#300 is the control recipe for Scones/Biscuits and #400 is the control 
recipe for muffins. 

The researcher also made Mu<"flns a standard recipe adapted from the bock of 
Culinary Arts Institute. Muffins were also made in an experimental basis by 
replacing milk with Magnacal. 

The recipes generated due to the experiments are appended in the appendix 
for further evaluation in order to produce and culture more about high density 
foods in the Space Program. 
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DISCUS SION 

Sensory methods in which palatability is evaluated by a panel of judges are 
essential to most food experiments because they answer all the important ques- 
tions of how a food should taste, smell, look, and feel. It is axiomatic that 
only a food which is prepared in a reproducible manner is worth the effort of 
scoring for it is only when the conditions of preparations are carefully con- 
trolled and defined that differences in quality can be a tribute to known 
variables. Careful sampling of the food is also necessary (20). 

All of the food samples were judged by trained panelists (JSC Food Laboratory). 
All of the foods were labeled and coded carefully, and all foods were evaluated 
by the sensory method on the basis of the following characteristics: 

a. Appearance 

b. Color 

c. Odor 

d. Flavor • 

e. Texture 

f. Overall 

A nine-point hedonic scale was used to evaluate the products: 



9-Like Extr*::nely 
8-Like Very Much 
7-Like Moderately 
6-Like Slightly 
5-Neither Like or Dislike 
4-Dislike Slightly 
3-Dislike Moderately 
2-Dislike Very Much 
1 -Dislike Extremely 


Shelf-life of space food must be considerably longer than most conventional 

foods to allow for tine to permit testing, evaluation and packag f-.g of flight 

I 

foods. These are time consuming operations and require that sufficient tine be 
allocated before each mission. 

The experimental Scones/Biscuits and Huff ins required a shelf-life study, 
and stored at three different temperatures for three different time lengths. 

The products are stored at a temperature of 40°, 70° and 100° Farenheit for 15, 
30, and 90 days. 

After shelf-life study, the high density Scones/Biscuits and Muffins will 
again be tested by the sensory evaluation panel. The products are also tested 
for any kind of visible changes, development of mold, etc. 

. Nutritional analysis of the high density foods are also in progress in the 
JSC Biomedical Laboratory. The results of the analyses will be completed within 
the next two months. 
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RESULTS 


Table I. Scones/Biscuits vn‘ th High Density Liquid 





Magnacal 

Sumacal Micro- 
lipid + 
Magnacal 

Hicro- 

lipid 

Micro- 

lipid 

with 

Shortening 

Control 

Overal 1 
Rating 

7.0-1 .2 

5. 8-1. 3 6. 2-0. 9 

6.5^1.! 

5. 2^1. 7 

6. 5^1. 3 

Range 

4-8 

3-8 5-8 

5-9 

2-8 

4-9 

- 

T^ole 

II. Muffins with High 

Density Liquid 



Control Magnacal Muffin 


Overs 1 1 
Rating 

6. 7-0.9 


6. 8-1. 2 


Range 

5-8 


4-8 



A shelf-life study will be conducted only about the product which is rated 
above five in the nine-point hedonic scale. 
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SUMt-IARY AND CpNCLUS ION 

Space flight poses special problems in food technology. Calorie-dense, 
nutritious and palatable, balanced diet in a suitable fonn jire necessary for 
space flight. Food must influence with packaging, water supply, waste and 
residue disposal systems and with food preparation and serving equipment. Food 
must be adequately protected against deterio»*ation, chemical, physical and 
biological. Also, all food must be adapted to consumption at Zero-gravity 
( 3 . 8 , 9 ). 

The conclusive report about high density foods -- particularly for Space 
Operations Center -- is yet to come. 
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MUFFINS 


2 cups sifted flour 

4 teaspoons baking powder 

1/2 teaspoon salt 

1/4 cup sugar 

1 egg, beaten 

1/4 cup melted shortening 

1 cup milk 

M ethod of Preparation 

Sift dry ingredients together. Mix egq, shortening and milk together thoroughly. 
Combine mixtures, stirring just enough to dampen flour. Fill greased muffin 
pan 2/3 full. Bake at 400°F. 25 min. Makes 12-15 muffins. 


MAGimCAL_ MUFFjj^S 

2 cups Bisquick 
1 cup Magnacal 
1 egg 

1/4 cup brown sugar 
1 teaspoon cinnamon 
1/2 cup )ple (finely chopped) 

1 tablespx.jn cooking oil 
dash of salt 

Method of Preparation 

Heat oven to 400°F. Mix ingredients vigorously for 30 seconds. Fill muffin cups 
2/3 full. Bake 15-18 min. (until golden brown). 

) 

NOTE: This is an experimental recipe. In this experiment milk component was 

substituted with high density fluid Magnacal. 


SCONES/BISCUITS 


2 cups sifted flour 
1 tablespoon sugar 
1/2 teaspoon salt 
4 teaspoons baking powder 
1/4 cup shortening 
1/2 to 2/3 cup milk 
1 egg well beaten 
Melted butter 
Sugar 

Method of Preparation 

Sift dry ingr 2 dients together; cut in shortening. Add milk to egg; arid to flour 
mixture gradually, adding more milk if necessary. Knead lightly on a floured 
surface. Roll into a round 1/2" thick. Cut into wedges, place on baking sheet, 
brush with melted butter, dredge with sugar. Bake at 400°F 15-10 min. Makes 

15 (4") Scones. 


NOTE: This is a standard recipe for making Scones. In each experimental recipe, 

Scones were made by replacing milk with Microlipid, Sumacal and flagnacal. 


| K 82 231 


CARDIOVASCULAR ADAPTATIONS IN WEIGHTLESSNESS: 

The Influence of In-flight Exercise Programs on the 
Cardiovascular Adjustments during Weightlessness 
and Upon Returning to Earth 


Charles H. Bennett, Ph.D. 

Department of Qnlogy 
Kentucky State University 


TABLE OF CONTENTS 


CONTENT page 

Table of Contents i 

List of Figures ii 

List of Tables iii 

Introduction 1 

Statement of the Problem 5 

Methods 5 

Cardiovascular evaluation during in-flight and postflight 

periods: Skylab M171 Experiment-Metabolic Activity .... 6 

Skylab M092 Experiment-Lower Body Negative Pressure 8 

Results 12 

Skylab Crewmembers Preflight Exercise Training 13 

Skylab Crewmembers In-flight Exercise Training 14 

Physiological Changes in Muscle Strength and Volume 

During Weightlessness 18 

:-,ylab Astronaut Cardiovascular Responses Recorded During 
Exercise Before, During, and Following Weightless. less: 

M171 Experimental Results 21 

Skylab Astronaut Cardiovascular Responses Recorded During 
the Application of Lower Body Negative Pressure(LBNP) : 

M092 Experimental Results 27 

Conclusions 31 

Acknowledgments 32 

References 33 


i 



LIST OF FIGURES 


'"igiire Page 

1. Levels of lower body negative pressure and time of individual 

phases of the lower body negative pressure protocol .... 10 

2. Skylab Bicycle Ergometer 15 

3. Skylab “MK-I" exerciser and several exercise methods 16 

4. Skylab Treadmill Exerciser 17 

5. Exercise related quantities on Skylab missions 18 

6. Average strength changes, arms: all Skylcb missions 19 

7. Average strength changes, legs: all Skylab missions 19 

8. Average leg volume change, postflight: all Skylab missions . . 20 

9. Heart rate responses for SKylab Crews, Level 3 exercise. ... 23 

10. Diastolic blood pressure. Level 3 exercise: all Skylab 

missions 24 

11. Changes in cardiac output and stroke volume in Skylab 

astronauts following periods of weightlessness 25 

12. Cardiovascular responses of the Skylab 4 Pilot during 

lower body negative pressure(LBNP) test 35 days 

prior to flight 27 

13. Cardiovascular responses of the Skylab 4 Pilot during 

lower body negative pressure(LBNP) test on mission 
day 10. Presyncopal symptr~s led to termination of 
the test 27 


i i 



LIST OF TABLES 


Table Number Table Content Page 

I. Skylab Medical Experiments 4 

II. M171 Experimental Protocol 7 

III. M171 Performance Summary - Three Skylab Missions 7 

IV. M171 Physiological Measurements 8 

V. M092 Performance Summary - Three Skylab Missions 9 

VI. Preflight Heart Rates Among Skylab Astronauts 13 

VII. In-flight Quantitative Personal Exercise Summary 15 

VIII. Daily Personal Exercise - Skylab 4 Crewmembers 18 

IX. Skylab 4 M171 Data Summary 22 

X. Plasma Volume - Percent Change from Premission Value ... 26 

XI. Differences Between Mean In-flight Values for Heart Rate 
and Blood Pressure of t.'ie Skylab 4 Crewmen during Rest 
and 50 mm Hg Phase of Lower Body Negative Pressure from 
Corresponding Mean Values during Preflight Tests .... 30 


iii 



CARDIOVASCULAR ADAPTATIONS IN WCIGHTLFSSNESS: 

The Influence of In-flight Exercise Programs on the 
Cardiovascular Adjustments during Weightlessness 
and Upon Returning to Earth 

Charles H. Bennett, Ph.D. 

Department of Biology 
Kentucky State University 

£1 

The American Space Program has provided an environment for the bio- 
medical investigation of physiological systems in the absence of gravity. 
Through use of the first long-term space station, the United States has 
acquired giant achievements in knowledge concerning the physiological 
effects of increasingly extended periods of exposure to the space flight 
environment and in det*?rmining how well man can function while performing 
tasks in space. 

Space medicine studies using experimental animals were initiated 
p ior to 1959. Limited medical studies and observations on men in space 
were initiated in the United States with the Project Mercury Program 
(1961-1963). This program served to dispel many basic concerns regarding 
the frailties of the human space explorer. It was shown that man could 
operate effectively during the acceleration periods of launch and entry, 
and that he could adapt to the weightless environment and perform useful 
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tasks. The biomedical studies conducted during the Gemini Program( 1965- 
1966) were oriented toward evaluating the magnitude of flight-related 
changes first noted following the Mercury flights. Emphasis was placed 
upon evaluation of the cardiovascular system, since the principal changes 
observed during the Mercury flights involved alterations in cardiovascular 
reflexes that regulate blood flow in the presence of the Earth's gravita- 
tional field. 

Several in-flight measurements or experiments were accomplished on 
the Gemini missions, as well as preflight and postflight studies. These 
investigations confirmed the postflight orthostatic intolerance observed 
in Met'cury and extended the findings to include moderately decreased post- 
flight exercise tolerance and red cell mass, minimal loss of bone calcium 
and muscle nitrogen, and the high metabolic cost of extravehicular activ- 
ity. The few changes which occurred were moderate but reversible. 

The Apollo Program( 1968-1973) originally included the conduct of 
a series of medical studies for the early orbital missions. After the 
Apollo ?04 accident, the decision was made to delete the medical studies 
and to dedicate all resources to the complex lunar landing program. Con- 
sequently, medical studies were primarily conducted with the Apollo crew- 
men before and after each flight. Generally, they confirmed the Gemini 
findings of postflight dehydration and weight loss, postflight decrease 
in orthostatic tolerance, and postal ight reduction in exercise capacity. 

In addition, the decreased red cell mass and plasma volume noted in 
Gemini were confirmed, but were less pronounced in Apollo. 

The three principal objectives of the Skylab Program( 1973-1974) were 
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the study of man, his Earth, and his Sun. Skylab, at its inception called 
the Apollo Applications Program, was a natural and necessary follow-on to 
the Gemini and Apollo Programs. Development of medical experiments, initi- 
ated in the mid-1960's, included a decision to design the experimental 
program along classical lines of medical and physiological research. The 
results from the Gemini and Apollo flights influenced the planning and 
placement of emphasis for the new program. The experimental protocols 
developed to study the cardiovascular, musculoskeletal, hematologic, 
vestibular, metabolic and endocrine systems in the body, with few excep- 
tions, remained unchanged throughout the Skylab Program. 

Twelve Skylab medical experiments(Table I) were designed to provide 
an indepth study of individual body systems and at the same time provide 
an overlap to give coniprehensi ve understanding of man's reaction to long- 
term weightless flight. The inclusion of major in-flight medical experi- 
ments provided the capability to study physiological responses during 
exposure to weightless flight as opposed to the pre- and postflight 
studies Ls carried out in the Apollo and Gemini flights. This paper 
will investigate the cardiovascular changes recorded during in-flight 
exercise and metabolic studies as well as during the application of 
lower body negative pressure and compare these responses with pre- 
and postflight values. ?ostflight cardiovascular responses following 
the extended duration flights of Skylab will also be examined in order 
to determine if the in-flight exercise programs selected by the Skylab 
crewmembers reduces the cardiovascular deconditioning observed during 
previous space flights. 
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TABLE I. SKYLAB MEDICAL EXPERIMENTS 


1 

Number 

Experiment 

M071 

Mineral balance 

M073 

Bioassay of body fluids 

M074 

Specimen mass measurement 

M078 

Bone mineral measurement 

M092 

Lower body negative pressure 

M093 

Vectorcardiogram 

MllO 

Hematology/iiiununol ogy 

M131 

Human vestibular function 

M133 

Sleep monitoring 

M151 

Time and motion study 

M171 

Metabolic activity 

M172 

Body mass measurement 


Medical evaluations after Gemini and Apollo flights demonstrated 
reduced orthostatic tolerance in virtually all crewmen specifically tested. 
Orthostatic tolerance describes the ability of the cardiovascul ar system 
to maintain an adequate level of arterial blood pressure, which provides 
the 'driving force* for maintaining organ system blood flow, during changes 
in body orientation in the earth's gravitational field. During changes in 
body position on earth, namely following movement from a supine to a 
vertical heads-up orientation, neurohumoral cardiovascular adjustments 
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occur to preserve arterial blood pressure through alteration in cardiac 
output(determined by the product of heart rate and stroke volume) and 
the peripheral resistance to blood flow. This 'short-term' regulation 
of cardiovascular homeostasis is further supplemented by the 'long- 
term' renal regulation of circulating fluids and electrolytes deter- 
mined in part by the release of antidiuretic hormone(ADH) , regulated by 
the Henry-Gauer reflex, as well as the renin-angiotensin-aldosterone 
regulatory system. The diminished ability of the cardiovascular system 
to function effectively against gravitational stress following exposure 
to weightlessness, while usually mild and not operationally significant 
during previous missions, sometimes resulted in pronounced increases in 
heart rate and decreases in pulse pressure during orthostatic testing. 

The magnitude of this postflight loss of orthostatic tolerance showed 
no clear correlation with flight durations ranging between 4 and 14 days. 
This caused concern for the potential postflight orthostatic intolerance 
in the crewmen of the extended Skylab missions, whose durations were 
planned for 28 days(Skylab 2), 59 days(Skylab 3), and 84 days(Skylab 
4 ). 


Statement of the proble m 

Postflight orthostatic intolerance and reduced physical work capacity 
have been reported following periods of weightlessness, bed rest, and 
water immersion. A quantitative assessment of the effects that in-flight, 
as well as preflight, exercise programs has on the cardiovascular responses 
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recorded during Skylab missions and following return to the earth's 
environment has not been performed. This report will attempt to deter- 
mine if cardiovascular differences exist between Skylab astronauts 
following the application of physiological stress(e.g., lower body 
negative pressure, exercise, orthostatic stress) that can be related 
to the crewmembers pre- and in-flight training programs. 

Methods 

Cardiovascular evaluation during in-flight and postfliqht periods : 

§ HyJ ^ 1 L Experiment - Me tabolic Activi ^ 

Twenty of the 27 Apollo crewmen tested following spaceflight ex- 
hibited a statistically significant decrease in their tolerance for 
exercise. Although this response was reversible within 1 to 2 days, it 
became obvious that man could not be conmitted to long-duration space 
flight until the magnitude and time course of these changes could be 
established and the underlying physiological mechanisms understood. The 
primary objective of the M171 Metabolic Activity experiment was to deter- 
mine whether man's metabolic effectiveness while performing mechanical 
work was progressively altered by exposure to the space environment. The 
secondary objective was to evaluate the M171 bicycle ergometer as an in- 
flight crew personal exerciser. 

A detailed description of the experimental hardware utilized for the 
M171 experiment has been reported elsewhere(ref . 1). The M171 Experiment 
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Protocol is shown in Tabic* II, while the M171 Performance Summary is 
listed in Table III. 


TABLE II. M171 EXPERIMENTAL PROTOCOL 


Time 

Exercise protocol 

5 minutes 

Rest 

5 minutes 

25 Percent of Maximun (level 1) 

.^2 

5 minutes 

50 Percent of Maximum V (level 2) 

.°2 

5 minutes 

75 Percent of Maximum V (level 3) 

O 2 

5 minutes 

Recovery 


Performed by each of three crewmen: Five times in 28-day 
mission; Eight times in 59-day mission; Twelve times in 84- 
day mission . 


TABLE III. M171 PERFORMANCE SUMMARY - THREE SKYLAB MISSIONS 


Skylab 

mission 

Prefl i ght 

In-flight 

Postflight 

Total 

2 Commander 

5 

6 

8 

19 

Scientist Pilot 

5 

6 

7 

18 

Pilot 

5 

7 

9 

21 


15 

19 

24 

58 

3 Commander 

7 

9 

8 

24 

Scientist Pilot 

7 

9 

8 

24 

Pilot 

7 

9 

8 

24 


21 

27 

24 

72 

4 Comniander 

8 

12 

8 

28 

Scientist Pilot 

8 

12 

8 

28 

Pilot 

8 

12 

8 

28 


24 

1 36 

24 

84 

Total: 3 Missions 60 

82 

72 

214 
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Results from this experiment, determined during pre-, in-, and post- 
flight periods, provided cardiovascular and pulmonary data(Table IV) which 
reflect adaptation as well as exercise-induced responses during weight- 
lessness. 


TABLE IV. M171 PHYSIOLOGICAL MEASUREMENTS 


Raw data 

Derived data 

Ergometer Work Level (watts) 

Respi ratory: 

Respi ratory : 

Respiratory exchange ratio. 

Oxygen Consumption. 

Carbon Dioxide production. 
Minute volume. 

Oxygen pulse. 

Pulmonary efficiency(Vr/V_ ). 

t 

Cardiovascular: 

Mechanical efficiency (V q /watt). 

Cardiovascular: ^ 

Mean arterial pressure. 

Total peripheral resistance. 
Arterial-Venous O 2 difference. 

ECG/VCG. 

Systolic blood pressure. 
Diasystolic blood pressure. 
Cardiac output.* 
Vibrocardiogram.* 

Carotid pulse.* 

Stroke volume. 


*P>"eflight and postflight only. 
ECG, Electrocardiograph. 

VCG, Vectorcardiograph. 


Skylab M0 9 2 Experiment - Lower Body Negative Pressure(LBNP) 

The objective of the lower body negative pressure experiment, designa- 
ted M092, was to determine the extent and the time course of changes in 
orthostatic tolerance during the weightlessness of space flight and to 
determine whether in-flight data from the experiment would be useful in 
predicting the postflight status of orthostatic tolerance. The equipment 
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design for supporting these studies has been described in detail (ref. 2 ). 
Preflight baseline data were acquired from the Skylab crewmen over a 4^^ 
to 5 month period prior to launch, with three tests during the last 6 
weeks prior to launch. In-flight tests were conducted usually at 3- or 
4-day intervals while postflight tests were carried out daily at first 
and then at increasing intervals of time over a period of approximately 
2 months. Table V lists the number of lower body negative pressure tests 
conducted for the three Skylab missions. 


TABLE V. M092 PERFORMANCE SUMMARY - THREE SKYLAB MISSIONS 


Skylab 

mission 



Preflight 

In-fl ight 

Postfl ight 

Total 

2 

Commander 


6 

7 

8 

21 


Scientist 

Pilot 

G 

7 

8 

21 


Pi lot 


6 

8 

8 

22 




18 

22 

24 

64 

3 

Commander 


6 

16 

8 

30 


Scientist 

Pilot 

5 

17 

8 

30 


Pilot 


5 

16 

8 

29 




16 

49 

24 

89 

4 

Commander 


6 

22 

9 

37 


Scientist 

Pilot 

7 

22 

9 

38 


Pilot 


7 

23 

9 

39 




20 

67 

27 

114 

Total : 

3 Missions 


54 

138 

75 

267 


Basic measurements during all tests with the lower body negative 
pressure device included blood pressure at 30-second intervals, from an 
automatic system which detected and analyzed Korotkoff sounds, heart rate 
continuously from one component of a Frank lead vectorcardiogram, and 
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percentage change in calf volume continuously from capacitive plethysmo- 
graphic bands encircling the legs. Prior to positioning these bands, a 
manual measurement was made of circumference of the largest portion of the 
calves, which also corresponded to the position where the left band was to 
be placed. The right band served to measure changes in capacitance due to 
alterations of temperature and humidity within the negative pressure device. 

The lower body negative pressure protocol was identical to that adopted 
for Apollo studies. The first and last 5 minutes of the 25-minute test v/ere 
at ambient atmospheric pressure to provide data from resting control and re- 
covery periods, respectively. The 15-minute stress period consisted of five 
distinct levels of negative pressure applied sequentially: 8 and 16 mm Hg 
negative pressure for 1 minute each, 30 mm Hg for 3 minutes, and 40 and 50 
mm Hg negative pressure for 5 minutes each. This profile is shown in fig. 1. 
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Figure 1. Levels of lower body negative pressure and time of individual 
phases of the lower body negative pressure protocol. 
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The hemodynamic changes induced by lower body negative pressure, 
like the orthostatic stress produced by tilting, depend upon the pooling 
of blood in the caudal part of the body. Depending upon the magnitude of 
the negative pressure and the time of exposure, 300 to over 800 ml of 
blood may be removed from the effective circulating blood volume. Central 
venous pressure is reduced, and the return of blood to the right heart 
(venous return) is diminished. Stroke volume and cardiac output, compensa- 
ted by an increase in heart rate, is reduced. Lower body negative pressure 
has been reported to cause a smaller increase in heart rate than upright 
tilting for equivalent reductions in cardiac output. The absence of stimu- 
lation of carotid baroreceptors due to postural change has been offered 
as an explanation of this difference. A decrease in blood flow to the 
upper extremities and concomitant increase in venous tone of the upper- 
extremity veins has been reported in both situations to compensate for the 
contracted effective circulating volume. 

A minor reduction in orthostatic tolerance, as observed following 
space flight, is accompanied by an excessive increase in heart rate, an 
excessive narrowing of pulse r>ressure, and a fall in systemic arterial 
blood pressure. The failure of cardiovascular compensation to gravity leads 
to the so-called vasodepressor reaction, the manifestations of which are 
presumably due to an overwhelming increase in parasympathetic nervous 
system activity(ref . 3). This reaction is characterized clinically by 
pallor, nausea, dimming of vision, sweating and eventually loss of con- 
sciousness, arising from a precipitous fall in arterial blood pressure 
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secondary to bradycardia(decreased heart rate) and a decrease in peri- 
pheral vascular resistance. 

There is no doubt that the decrease of blood volume and reabsorption 
of fluid transudate from tissues of the lower extremities during exposure 
to weightlessness would diminish orthostatic tolerance, since a decrease 
of blood volume in a .normal active individual from any cause, such as 
blood loss or dehydration, will result in a strain being placed on normal 
mechanisms required to maintain cardiovascular integrity during orthostatic 
stress. Orthostatic intolerance is still present in astronauts returning 
from space even after blood volume returned to preflight levels, suggest- 
ing that cardiovascular adaptation during weightlessness is influenced 
by more than one physiological regulatory system. 

One important question addressed in this report is the effect of 
inflight exercise, as well as preflight conditioning, upon orthostatic 
tolerance, lower body negative pressure responses, and exercise capacity 
changes resulting from space travel. Since skeletal muscle loses tone, 
strength, work capacity, and mass during periods of inactivity such as 
weightlessness, inflight attempts to reduce these loses might assist in 
maintaining cardiovascular conditioning and therefore help reduce ortho- 
static intolerance during return to a 1-g environment. A quantitative 
assessment of Skylab data related to these points has not been performed. 

R esul ts 


This section will evaluate the results from the Skylab missions in 
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the areas of crewmen pre- and in-flight training methods and exercise 
levels, and how training may have affected the cardiovascular responses 
observed following lower body negative pressure and exercise during and 
following spaceflight. 

Sky lab Crewmembers Preflig ht Exercise Training . 

Preflight exercise methods among Skylab astronauts varied consider- 
ably, and the training effects are difficult to document. One can evalu- 
ate the state of cardiovascular conditioning, on a first order basis, by 
measuring preflight resting and stressed heart rates among the crewmembers. 
These are shown in Table VI, with the stressed heart rates determined 
following application of -50 mm Hg lower body negative pressure{64 pre- 
flight determinations an.^ng all Skylab astronauts). 


TABLE VI. PREFl.IGHT HEART RATES AMONG SKYLAB ASTRONAUTS 
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Resting heart rates presented in Table VI were determined by aver- 
aging six to eight preflight measurements on each individual crewmember. 
Based on these data, one would predict that the Skylab crewmembers had 
undergone significant preflight conditioning, regardless of the method 
used. 

The preflight exercise programs selected by individual crewmembers 
of the Skylab missions can be separated into one of three closses, al- 
though some overlap occurs. Class 1 includes a training program based 
primarily on distant running, Class 2 applies to crewmembers who pri- 
marily exercised by playing racquetball , and Class 3 included the crew- 
members who combined the Class 1 and 2 programs, although at lower levels 
of intensity, with weights, swimming, gymnastics, etc. Again, there is 
considerable overlap among some crewmembers, but where possible, those 
who can be placed into one of the three clases will be identified. 

S kylab Crewmembers In-flight Exercise Training . 

Although in-flight exercise profiles and work levels varied among 
crewmembers. Table VII shows that following normalization of work data 
to the crewman's body weight, the Skylab 3 crewmembers exercised about 
107 percent more than the Skylab 2 crew and the Skylab 4 crewmembers 
exercised 130 percent more than the Skylab 2 crew. The effects that the 
increased in-flight exercise load had on postflight recovery will be 
discussed later in this section. Special emphasis will be placed on 
measuring the rate of return of cardiovascular responses to preflight 
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TABLE VII. IN-FLIGHT QUANTITATIVE PERSONAL EXERCISE SUMMARY 


Skylab 

mission Crewman 

{l)Total 
(watt min) 

(2)Daily avg 
(watt min) 

(3)Daily avg 
(watt min/kgm 
body weight) 

Avg 

2 Commander 

■19 

2855 

47 


Scientist Pilot 


1618 

21 

31.3 

Pilot 

55795 

1993 

26 


3 Commander 

228581 

3874 

58 


Scientist Pilot 

214645 

3638 

62 

65.0 

Pilot 

386193 

6545 

75 


4 Commander 

349210 

4108 

62 


Scientist Pilot 

469420 

5523 

80 

72.3 

Pilot 

414760 

4879 

75 



(1) Includes M171 experiment tests and personal exercise 

(2) Based on 28-day Skylab 2 mission, 59-day Skylab 3 mission, and 84- 
day Skylab 4 mission. 

(3) Based on mean in-flight body weight. 


Skylab 2 crewmembers primarily utilized the M171 bicycle ergometer(fig. 


2) for in-flight exercise. 


Fig. 2. Skylab Bicycle Ergometer 
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When their postflight return of exe^'cise tolerance to preflight status 
exceeded the time for recovery observed in Apollo crewmen, it was decided 
to add a mass-produced commercial device, called Mini Ciym, which was ex- 
tensively modified and designed "MK-I"(fig, 3). This device was developed 
to provide adequate exercise to arms, trunk and legs. A series of exer- 
cises utilizing "MK-I" are presented in figure 3. 



A second device, designated "MK-H", consisted of a pair of handles 
between which up to five extension springs could be attached, allowing 
maximum forces of 25 pounds per foot of extension to be developed. These 
two devices were first flown on Skylab 3, and food and time for exercise 
was increased in-flight. 

Although the bicycle ergometer is an excellent machine for aerobic 
exercise and cardiovascular conditioning, it could not develope either the 
type or level of forces to maintain strength for walking under 1-g. Be- 
cause of the leg muscle strength and volume loses shown following Skylab 
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2, work wds started to develope an in-flight treadmill exerciser(fig. 4). 



Fig. 4. Sky lab Treadmill Exerciser 

This device a lowed walking and running, plus other isotonic training, 
under forces equivalent to gravi ty( 1-g) . Four rubber bungees provided an 
equivalent weight of 80 kilograms and were attacned to a shoulder and 
waist harness. By angling the bungees, an equivalent to a slippery hill 
is presented to the subject who must climb it. High loads were placed on 
some leg muscles, especially the calf, and fatigue was rapid such that the 
device could not be used for significant aerobic work. 

On Skylab 4, the crew used the bicycle ergometer at essentially the 
same rate as Skylab 3, and the "MK-I and II" exercisers. In addition, they 
typically performed 10 minutes per day of walking, jumping, and jogging on 
the treadmill. Food intake had again been increased. Daily exercise protocols 
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for the Sky lab 4 crewmembers are shown in Table VIII. 


TABLE VIII. DAILY PERSONAL EXERCISE - SKYLAB 4 CREWMEMBERS 


Exercise 

Commander 

Scientist 

Pilot 

Pilot 

Leg Ergome try (watt min). 

5000 

8337 

6000 

Minigyin(Total repeti tions) . 

100 

200 

200 

Springs(Total repetitions). 

75 

0 

120 

Torso Isometrics(Total 

20 

0 

20 

repeti tions) 




Treadmil 1 ; 




Walk (min) . 

10 

0 

0 

Run(mi n) . 

1 

0 

0 

Springs (Repeti tions) . 

300 

1000 

100 

Too Rises(Repetitions). 

200 

200 

75 


S bxMV3^th_^ 

Figure 5 presents a suimiary of Skylab exercise levels as well as the 
changes in body weight, leg strength and leg volume following weightlessness. 
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Fig. 5. Exercise related quantities on Skylab missions 
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From figure 5, it is apparent that the increased exercise loads, along with 
the increased caloric intake from dietary supplements, allowed the Skylab 
4 crewmembers to undergo a significantly smaller decrease in body weight, 
leg volume and leg strength. One crewmember of Skylab 4. the commander. 

was only the second American astronaut to gain weight in space, even after 
a spaceflight of 84 days duration. 

Even prior to muscle testing following re-entry from space, it was 
obvious that the Skylab 4 crew was in surprisingly good condition. They 
stood and walked for long periods without apparent difficulty on the day 
after recovery in contrast to the earlier missions. Results of the muscle 
testing, which will be presented later in this section, confirmed a sur- 
prisingly small loss in leg strength after almost three months in weight- 
lessness(fig. 5). A summary of the exercise and strength testing shown in 
average values for the three missions is presented in figures 6 and 7. 



Fig. 6. Average strength changes, 
arms; all Skylab Missions. 



Fig. 7. Average strength changes 
legs: ail Skylab Mission 
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Relatively small losses in arm strength(fig. 6), as compared to 
leg strength{fig. 7), is shown in all missions. This reflects the fact 
that in space ordinary work provides loads for the arms that are relative- 
ly much greater; the legs receive virtually no effective loading . The 
increase in leg extensor(antigravi ty) muscle strength observed on Skylab 
4(fig. 7) undoubtedly contributed to their early postflight recovery. 

A common measure of muscle condition is muscle size. Figure 5 shows 
that muscle size of the Skylab 4 astronauts was reduced only one-half the 
volume of the shorter flights. Recovery of leg volume to preflight status, 
shown in figure 8, illustrates that the Skylab 4 crewmembers displayed 
faster recovery rates than the Skylab 2 and 3 crewmen. 



Fig. 8. Average leg volume change, postflight; all Skylab missions 

It is interesting to note that the amount of exercise performed in- 
flight was inversely related to the length of time required postflight for 
return to preflight status. Except for some isolated individual responses 
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in cardiac output and stroke volume data, all other cardiopulmonary and 
work capacity parameters retunred to normal in approximately 18 - 21 days 
for the Skylab 2 crewmemebers , 7 days for the Skylab 3 crew, and 4 days 
for the Skylab 4 crew. Paradoxically, the Skylab 3 crewmen, who displayed 
leg volume and strength decrements equivalent to those recorded for the 
Skylab 2 crewmembers, returned to preflight status at a rate almost equal 
to the Skylab 4 crew. 

Skylab Astronaut Cardiovascular Responses Recorded During Exercise Before, 
During, and Pol l owj ng Weightlessness; M171 Experiment . 

Table IX presents the cardiovascular responses of the Skylab 4 crew 
collected from the M171 experiment. The measured parameters included heart 
rate, systolic blood pressure, and diastolic blood pressure collected at 
rest, during the final 3 minutes of 5-minute level 3 exercise, and the re- 
covery values were those obtained for the second minute during the 5-minute 
recovery period. 

All tliree crewmembers of the Skylab 4 mission experienced a signi- 
ficant decrease in in-flight recovery heart rate. It would appear that 
following exercise in weightlessness, heart rate values return toward 
resting values more rapidly, and to a lower level, than 1-g responses. 
In-flight heart rate increases during level 3 exercise were not signi- 
ficantly different from those recorded preflight. When the Skylab 4 
crewmembers returned to e th, they all exhibited significantly elevated 
heart rate values in comparison to pre- and in-flight recorded resting 
heart rate values. This was a common observation among a majority of 
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Skylab astronauts following return to earth. The heart beats at a faster 
rate upon returning to earth in an attempt to maintain cardiovascular 
homeostasis induced by a reduction in cardiac output. 

TABLE IX. SKYLAB 4 M171 DATA SUMMARY 


Variables 



Heart Rate(BPM) 
Rest Level 3 Recovery 

Rest 

SBP(mm 

Level 

Hg) 

3 Recovery 

Rest 

DBP{mm Hg) 
Level 3 

Rec( 

CDR 










Prefl ightOO 

66 

157 

112 

96 

192 

179 

67 

71 

6i 

In-flight(X) 

66 

152 

87* 

97 

195 

181 

59* 

56* 

6 

Postfl ight(X) 

76* 

163 

109 

106 

195 

170 

72 

67 

7: 

SPT 










Prefl ight(X) 

64 

164 

104 

127 

204 

186 

84 

55 

6i 

In-flight(X)^ 

62 

166 

92* 

119* 

200 

174 

74* 

52 

6 

Postfl ight(X) 

74* 

167 

102 

123 

198 

189 

78* 

51 

6: 

PT 










Prefl iqht(X) 

54 

147 

114 

115 

204 

188 

72 

60 

61 

In-fl ight(XJ 

53 

147 

91* 

115 

200 

186 

64* 

51* 

61 

Postfl ight(X) 

65* 

156* 

118 

125* 

213 

204 

74 

59 

6! 


*0utside the preflight 95% confidence limit. 

BPM, beats per minute; SBP, Systolic Blood Pressure; DBP, Diastolic Blood 
Pressure; CDR, Commander; SPT, Scientist Pilot; PT, Pilot. 


Table IX also indicates that spaceflight significantly lowers diastolic 
blood pressure during resting conditions in all crewmembers. This para- 
meter is also significantly reduced in two crewmen during level 3 exer- 
cise periods. Diastolic blood pressure changes will be discussed later 
in this section. 


Generally, the in-flight and postflight responses to exercise by the 
crews of Skylab 2, 3 and 4 were similar. In-flight, some subtle, isolated 
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differences were seen. There were no trends observed which would indicate 
a degradation in the exercise responses of the crews. The Skylab 4 crew 
exhibited a significant in-flight decrease in recovery heart rate but not 
in resting(si tting position) heart rate. The Skylab 2 .crew, on the other 
hand, exhibited decreases in both parameters while the Skylab 3 crew ex- 
hibited no changes in either. Figure 9 displays the heart rate responses 
of all Skylab crewmembers during pre-, in-, and postflight level 3 exer- 
ci se. 


t St§nd»rd 



Fig. 9. Heart Rate Responses for Skylab Crews, Level 3 Exercise. 

The data presented in figure 9 indicate that postflight recovery heart rates 
tend to return to preflight(and in-flight) values with longer duration flight 
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All Skylab 2 crewmembers had significantly higher heart rate values 
following postflight exercise, while two Skylab 3 and only one Skylab 4 
crewmember displayed this response. 

Figure 10 displays the diastolic blood pressure responses for all 
Skylab astronauts during level 3 exercise periods conducted pre-, in-, 
and postflight. 



*Signi(ic»nt ditterence fioin biiebne ip < 0 06) 


Fig. 10. Diastolic blood pressure, level 3 exercise: all Skylab missions, 


Figure 10 shows an interesting trend in diastolic blood pressure among the 
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Skylab astronauts during spaceflight. Five of six Skylab astronauts from 
the longer duration missions(59- and 84-days) showed significantly lowered 
diastolic blood pressure during heavy exercise{level 3) while in the 
weightless condition. 

One potentially dangerous cardiovascular adjustment noticed in all 
Skylab astronauts upon retunring to earth is a significant reduction in 
cardiac output and stroke volume(fig. 11). The exact explanation for this 
adjustment is not clear, but may reflect a reduced circulating blood volume 
(Table X) which therefore limits venous return of blood to the heart. Posi- 
tive chronotropic(increases in heart rate) and inotropic{increases in 
myocardial contractility) responses of the heart may also be reduced 
following extended periods of weightlessness. 


CARDIAC OUTPUT CHANGES IN SKUAB ASTRONAUTS STROKE VOLUME CHANGES IN SKYLAB ASTRONA 




Figure 11. Changes in Cardiac Output and Stroke Volume in Skylab 
Astronauts Following Periods of Weightlessness. 


f 


One final observation in the area of cardiovascular adaptations during 
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weightlessness that awaits further study is the apparent loss of circula- 
ting plasma volume and what can be done to prevent this loss(Table X). 
This loss of circulating blood volume undoubtedly affects cardiovascular 
performance during spaceflight as well as following return to earth. As 
shown in Table X, longer duration missions produce greater losses in 
plasma volume. However, since the Skylab 4 astronauts were capable of 
retruning carviovascular parameters to preflight values at a rapid rate 
(4-days), their loss of plasma volume does not appear to limit cardio- 
vascular performance following return to earth. 


TABLE X. PLASMA VOLUME - PERCENT CHANGE FROM PREMISSION VALUE. 


First Skylab Mission 
Premission volume(ml) 

R + 0* 

R + 13 
R + 42 
R + 67 

Commander 

1 

St .ienti s t 
Pilot 

Pilot 


Mean 

3042 

Percent 

-2.5 

-1.2 

+8.5 

-5.7 

3506 

Percent 

-10.3 

-5.6 

-1.1 

3472 

Percent 

+2.6 
+14.1 
+18.6 
+ 17.0 

Percent 

(**-12.3) 

(**-2.5) 

(**+1.4) 

(**0.0) 

Percent 

-8.4 

-3.1 

+4.9 

-2.3 

Second Skylab Mission 






Premission volume(ml) 

3157 

2798 

3885 




Percent 

Percent 

Percent 


Percent 

R + 0* 

-18.4 

-9.1 

-11.8 


-13.1 

R + 14 

+0.1 

+ 14.7 

+2.0 


+5.6 

R + 45 

+2.8 

+ 11.7 

+6.8 


+7,1 

Third Skylab Mission 






Premission value(ml) 

3067 

3620 

3195 




Percent 

Percent 

Percent 


Percent 

R + 0* 

-15.7 

-19.2 

-12,9 


-15.9 

R + 14 

+8.6 

+7.4 

+ 13.0 


+9.7 

R + 31 

+5.4 

+ 17.7 

+5.9 


+10.0 


*R +, Recovery + day(s). **% change calculated using R + 67 day value. 
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ORtQINAL PAGE 13 
OF POOR QUALITY 


Skylab Astronaut Cardiovascular Responses Recorded During the Applicati on 


of Lower Body Negative Pressure(LBNP) - M092 Experiinent. 


An explanation of the changes in cardiovascular responses to lower body 
negative pressure in a weightless environment requires understanding of the 
manner in which the systems of the body and their functions adjust to the 
weightlessness of space flight. A more comprehensive understanding of thes?? 
adjustments will, in turn, require the correlation of massive volumes of 
data from all the Skylab experiments, a task of monumental proportions 
which has barely been started. 

Cardiovascular responses of the Skylab 4 Pilot during lower body nega- 
tive pressure test 35 days prior to flight(fig. 12) and on mission day 10 
(fig. 13) are shown for comparison. The mission day 10 test was terminated 
early due to presyncopal conditions developed by the Pilot. 



Fig. 12. Cardiovascular responses of 
the Skylab 4 Pilot during LBNP test 
35 days prior to flight. 


Fig. 13. Cardiovascular response: 
of the Skylab 4 Pilot during LB! 
test on mission day 10. Presync 
symptoms led to termination of 
test. 
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One area of interest that this report will not completely cover is 
the ability to predict the factors responsible for the presyncopal reactions 
demonstrated by Skylab crewmembers. No crewmen seemed inimune to this response 
but some seemed more susceptible than others. One factor which awaits further 
analysis is the influence of physical training on the develoment of pre- 
syncopal conditions during lower body negative pressure. Recent studies 
have indicated that distant runners may be more susceptible to developing 
presyncopal responses during lower body negative pressure than individuals 
who train using other form of exercise, as weightlifting or swimming(ref . 4). 

Otfier factors certainly influence the cardiovascular responses to 
lower body negative pressure. For example, none of the 13 instances of early 
termination of tests due to presyncopal symptoms occurred during tests con- 
ducted in the morning or within 7 hours of arising from sleep, although 
approximately one-third of all in-flight tests were conducted within that 
period. These 13 tests were associated with larger than usual calf volume 
increases(see e.g. , fig. 13) during the lowest negative pressure phases of 
the tests. This suggests that the pooling of large volumes of blood during 
the first few minutes of the test may so alter the effectiveness of compen- 
satory cardiovascular mechanisms as to render them incapable of adequate 
responsiveness to the greater stress later in the test. 

Among the three astronauts, the Commander of Skylab 2, the Pilot of 
Skylab 3, and the Commander of Skylab 4, whose in-flight tests were nomi- 
nally scheduled in the morning, only one presyncopal episode occurred and 
this was in the Commander during a test performed in the afternoon over 9 
hours after arising. These crewmembers also could be grouped into class 2 
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based on their preflight exercise training. 

The Pilot and Scientist Pil it of Skylab 4 exhibited a larger number of 
presyncopal responses during lower body negative pressure than a majority 
of the other Skylab astronauts. These two crewmembers, especially the 
Scientist Pilot, could be placed in Class 1 according to preflight exer- 
cise training. The Pilot had taken up distant running several years before 
flight and had experienced a significant weight loss during this period. 

When one examines the calf volume changes induced by lower body negative 
pressure, in-flight increases greatly exceeded those which occurred during 
preflight testing. Calf volume increases of the Skylab 4 Commander in pre- 
flight tests were relatively small in comparison to those of the Scientist 
Pilot and Pilot of that mission, which were larger than usually seen. This 
same pattern of difference continued throughout the mission. Whereas pre- 
flight calf volume increases at the end of the 50 mm Hg negative pressure 
phase averaged between 3 and 4 percent, during in-flight tests they reached 
values usually between 8 and 11 percent, in the Scientist Pilot and Pilot. 

In the Commander, calf volume increases averaged 2.4 percent preflight 
but usually reached levels ranging between 5 and 7 percent in-flight. 

The volume of blood pooled in the lower extremities during in-flight 
tests did not seem to correlate from test to test with the magnitude of heart 
rate or blood pressure change during lower body negative pressure. In general 
those of the nine Skylab astronauts with the greatest increase in calf volume 
during in-flight tests also showed the greatest increases in heart rate and 
changes in blood pressure. These data, for the Skylab 4 crewmembers, are shown 
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In Table XI. In addition, although correlation of heart rate increases with 
leg volume increases was not evident in preflight tests, the crewmen whose 
calf volume increases were greatest in preflight tests usually also showed 
the largest increases during the in-flight tests. 


TABLE XI. DIFFERENCES BETWEEN MEAN IN-FLIGHT VALUES FOR HtART RATE AND 
BLOOD PRESSURE OF THE SKYLAB 4 CREWMEN DURING REST AND 50 mm HG 
PHASE OF LOWER BODY NEGATIVE PRESSURE FROM CORRESPONDING MEAN 
VALUES DURING PREFLIGHT TESTS. 



CDR 

SPT 

PT 

Resting control: 




Heart ra te(bprn) 

+ 7.8* 

+13.3* 

+11.3* 

Systolic blood pressure(mm Hg) 

+ 1.8 

+ 2.5 

+ 1.0 

Diastolic blood pressure(mm Hg) 

-5.6* 

-2.0 

-2.5 

Pulse pressure (mm Hg) 

+ 7.4* 

•+4.5* 

+ 3.4 

Mean arterial pressure(njii Hg) 

-3.2* 

-0.5 

-1.3 

Stressed, -50 mm Hg: 




Heart rate(bpm) 

+12.2* 

+36.7* 

>■26.7* 

Systolic blood pressure(mm Hg) 

+6.2* 

-10.9* 

-1.6 

Diastolic blood pressure(ii«n Hg) 

-1.0 

-4.8 

+ 1.3 

Pulse pressure(mn) Hg) 

+ 7.1* 

-6.1 

-2.8 

Mean arterial pressure(mm Hg) 

+ 1.5 

t 1 1 

-6.8* 

+0.3 


CDR, Commander; SPT, Scientist Pilot; PI, Pilot; 
♦Significant to 0.05 level by Student's paired t-test. 


The in-fl:ght mean increase in heart rate during 50 mm Hg lower body negative 
pressure over resting rates for all nine Skylab crewmen averaged 20.4 beats per 
minute, a highly significant difference. The increase in heart rate during ortho- 
static sir- '>s has generally proven to be the best single index in the assess- 
ment of orthostatic tolerance. The pattern of fluctuations of resting and 
stressed heart rates and blood pressure varied for each crewman, but their 
magnitude and frequency were greater during the first part of the mission than 
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later. Their nature and significant is unknown, but, in the case of the Skylab 
4 crewmen, their prominence and duration appeared to decrease as cardiovascular 
responses to in-flight lower body negative pressure stabilized. 

The Skylab Lower Body Negative Pressure experiment demonstrated that loss 
of orthostatic tolerance had already developed by tiie time of the first tests 
after 4 to 6 days of flight. Cardiovascular responses to lower body negative 
pressure showed the greatest instability and orthostatic tolerance thi reat- 
est decrement during the first three weeks o.f flight. After approximately 5 
to 7 weeks, cardiovascular responses became more stable and evidence of im- 
proving orthostatic tolerance appeared. 

In-flight data from the lower body negative pressure experiment proved 
to be useful not only in predicting the early postflight status of ortho- 
static tolerance, but also in the in-flight assessment of crew health status. 
In-flight lower body negative pressure presented a much greater stress to the 
cardiovascular system than the same levels of negative pressure during pre- 
flight tests. 

At rest, in-flight mean resting heart rates, systolic blood pressures, 
and pulse pressures were typically increased while diastolic and mean arterial 
pressures decreased compared to preflight values in the majority of Skylab 
crewmembers. Differences in in-flight responses to lower body negative pres- 
sure stress from preflight responses included greater heart rate and leg volume 
increases in all crewmen and, in most, higher diastolic pressures and mean 
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arterial pressures and lower systolic blood pressures and pulse pressures. 

The successful completion of the 28-, 59- and 84-day Skylab missions 
showed that man can perform submaximal and maximal aerobic exercise in the 
weightless environment without detrimental trends in any of the physiologic 
data. Exercise tolerance during flight was unaffected. It was only after 
return to Earth that a tolerance decrement was noted. 

The rapid postflight recovery of orthostatic and exercise tolerance 
following two of the three Skylab missions appeared to be dirertly related 
to total in-flight exercise as well as to the graded, regular program of 
exercise performed during the postflight debriefing period. 

The postflight orthostatic intolerance and diminished exercise capacity 
are both related etiological ly to a decreased effective circulating blood 
volume at one-g, with consequent decreased venous return and cardiac output. 
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INTEBMATIONAL AEROSPACE ENGINEERING: 

NASA SHUTTLE and EUROPEAN SPACELAB 
By 

Roger E. Bilstein* 

The successful mission of the shuttle Columbia in April 1981 
brought a new sense of achievement to NASA and to America. The 
flight also sent an expectant tremor through the European Space 
Agency ESA, who plans to have a significant role in future 
shuttle flights. ESA designed and built the Spacelab, a major 
payload item for shuttle missions to come. Spacelab provides a 
shirt'Sleeve environment for non-astronauts , or "mission** 
specialists," to conduct a host of laboratory research projects 
in space. It offers a flexibility never available to scientists 
before the Shuttle/Space lab appeared on the scene. More import- 
antly, Spacelab gives western Europe, for the time, a leading 
role in manned space flights, eind represents a major turning-point 
for ESA and its contractors. A spokesman for ERNO, the European 
prime contra^'^tor, summarized the changed situation: "Award of 

the contract to build the space laboratory and its complex sys- 
tems meant that ERNO and its partners would tread completely new 
paths," the European contractor observed. "Their space activities 
had so far been concentrated on the satellite and launch sector. 
Now new procedures and methods had to be established and tackled 
for the first manned space flight system. Spacelab exceeded all 
the familiar dimensions,"^ 

"The opinions and conclusions are the author’s only, and do not 
reflect official NASA policy. This paper is part of a longer 
study still in progress, and much of it is still tentative; 
contact the author for permission to cite or quote. 
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Shuttle/Spacelab operations also represent a major change 
for NASA. The American space agency traditionally managed its 
own contractors in the development of space hardware. Although 
NASA occasionally launched satellite payloads for foreign cus- 
tomers, including ESA, these were essentially discrete projects 
with minimal integration problems. Spacelab changed all that, 
requiring major collaboration on the design and operation of the 
Shuttlc/Spacelab combination. Interaction between Europeans and 
Americans, against a background of long-established NASA pro- 
cedures, demanded diplomatic give-and-take involving top-level 
administrators, program managers, and engineers — from chancellories 
all the way down to the shop floor. This paper emphasizes 
Shuttlc/Spacelab developments through the Memorandum of Under- 
standing in 1973. 

IDEAS AND ORGANIZATIONS 

At least a decade before Sputnik in 1957, the idea for a 
reusable space transportation system appeared in American popular 
and technical articles alike. The complexities of such schemes 
proved discouraging; in any case, the drive to perfect ICBM 
weapons offered a reasonable alternative to boost scientific 
payloads into space. Sputnik shocked the U.S. into a heavy in- 
vestment in space research. The Russian- American space race that 
followed evolved in the spotlight of international attention and 
undersvx>red space research as a measure of national prestige. 

In advUtion, an active space program seemed to generate consider- 
able technological fall-out. Some leaders in Western Europe 
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began to worry about decline in both prestige and technical 

abilities unless they developed their own space program. France 

became particularly active in this respect and put its own 

Asterix satellite, aboard its own Diamant three-stage booster, 

in orbit by 1965. Britain bec 2 une interested in the possibilities 

of using its Blue StreaJc rocket as a first-stage booster for a 

multi-staged European launch vehicle. These and other strands 

of European activity finally resulted in the creation of the 

European Launch Development Organization (ELDO) and the Europe 2 m 

2 

Space Research Organization (ESRO) by the mid-60s. 

The European space organizations hoped to get into the space 
business in a way that would benefit Continental business and 
industry. ESRO, with 10 members, emphasized scientific payloads 
and Ssitellit^a development. ELDO concentrated on a larger booster 
that would end dependence on NASA rockets launched from Cape 
Kennedy. ELDO's principals included Britain, France, Germany 
(builders, respectively, ; of first, second, and third stages of 
a launch vehicle) , and Italy (coordinator of a satellite payload) . 
By 1969, after eight years of effort, Europe's united bid for an 
independent space capability failed to coalesce.^ The stage was 
set for European space research to move in different directions. 
Two lines of activity emerged; renewed attention to manned, i 
earth-orbital vehicles that eventually led to the Spacelab con- 
cept; an "alliance of convenience" with NASA. 

In the U.S., plans for a reusable Space Shuttle had already 
begun to finalize. Mercury, Gemini, and Apollo programs provided 
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incxoeslng understanding of manned space vehicles. High-speed 
aerodynamics and thermal problems became better understood 
through the X-15 and other experimental aircraft programs. A 
series of "lifting bodies” yielded design information on high- 
speed descent and landing. By 1970 » based on this experience » 
as well as prototype shuttle studies » the U.S. was committed 
to a shuttle transportation system. Plws called for a vehicle 
to ferry up to 12 passengers back and forth between Earth and 
a space station, itself built from modules carried in a shuttle 
cargo bay and assembled in orbit. ^ Through space transport 
studies similar to these, the idea for a European Spacelab be- 
gan to emerge. 

During the mid-1960s, European thinking about space explor- 
ation began to sharpen its focus through a 150-member industrial 
consortium called Eurospace. As early as September 1964 a 
Euroi-j'ace report summed up the group's thinking about a reusable 
spacecraft — an aerospace transporter — that had been under informal 
discussion for several years. The concept envisioned piloted 
booster and orbiter stages, the latter to have a two-man crew 
and to carry 2.5 tons into a 180-mile orbit. During a visit to 
the U.S. in the spring of 1965, Dr. Erhard Loewe, a vice-president 
of Tele funken, and Eurospace spokesman, discussed Europe's steike 
in the aerospace transporter. The U.S. and Russia had already 
committed themselves to a manned lunar landing, he observed. 

Europe did not have the technology for such liiner aspirations, 
but the aerospace transporter promised Europe a reasonable role 

ORIGINAL PAGE .3 
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in the near-Earth space environioent. Hore than that# the trans- 
porter would serve as a springboard for Europeim survival in a 

technical world. "We are • . • inf luenced by the conviction that# 
tin 

should we be^<able to go ahead with the realization of our ob- 
jective# we should fall far behind in the technical field#” 

Loewe empahsized. At the same •time# he left the door open for 
discussing joint U.S. -Europe space ventures. Loewe stressed the 
importance of space stations for the true exploitation of space# 
and noted that space stations could offer a very useful option 
for joint efforts.® 

Eurospace consideration of the aerospace transporter appar- 
ently included rather extensive paper studies and scale models 
of at least one proposed design. Junkers# the German aviation 
firm, fielded a design teeun including Boelkow-^twicklungen KG, 
and another consortium known as ERNO. Eurospace sources insisted 
that the transporter idea was not too ambitious for them, and 
that they could develop the idea by 1975-80. Talk of a modest 
space station alsc continued. If all this was to be achieved, 
many European aerospace experts realized that the current organi- 
zation for space exploration would have to be replaced by a 
different scheme. ESRO and ELDO continued to work independently. 
They shared laudable liaison connections, but still had no 
mechanism for joint planning. With no comprehensive policy, 
various European governments remained reluctant to fund programs. 
Christopher Ashton, an assistant technical director for ELDO, 
advocated more liaison at a high policy-making level, especially 
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a government-’to-governinent basis, to establish a coherent Euro- 
pean space program. Eurospace thinking tended to agree. . In 
short, ideas began coalescing for a single space agency — a Euro- 
pean NASA.^ 

In retrospect, thh activities and mood of European space 
interests in the Spring of 1965 assumed. considerable 

significance in light of subsequent developments. Design studies 
for the aerospace transporter and space station started European 
aerospace thinking about the details involved in the launch, 
operation, emd return of orbiting spacecraft; the life support 
systems and operational requirements for the manned transporter 
emd manned space station in orbit; the need for a unified space 
organir.ation to accon^lish all of this. The eventual creation 
of ESA (the European Space Agency) , and development of the Space- 
lab reflects a direct lineage to this period. Moreover, EBMO 
emerged as the principal Spacelab contractor. 

Over the next few years, a series of developments drew 
Europe and the U.S. together as aerospace partners. Even though 
France put a number of modestly sized satellites into orbit, 
beginning in 1965, & vigorous European program failed t ter- 
ialize. Coordination of efforts remained indecisive, and sub- 
stantial funds for extensive RSD never appeared. NASA spurned 
suggoKlions for joint research on the aerospace transporter as 
too much duplication of advanced U.S, work on the Shuttle. In- 
stcavl, NASA opened a different door for cooperation. In Febru- 
ary during a visit of West Germany's Chancellor, Ludwig 

Erhard, to the U.S., President Lyndon Johnson suggested collab- 

in a major space exploration effort. The U.S, suggested 
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proboi to the Sun and to Jupiter. Johnson possibly made the 
offer to siue-step the issue of West German control of nuclear 
weapons » but .Erhard, at least, appe£u:s to have encouraged thu 
space venture. Shortly thereafter, NASA Administrator, Jeunes 
Webb, and a high-level NASA executive group left for Europe to 
pursue the issue. 

The nature of subsequent talks seems to have shifted to- 
wards manned missions. NASA, during the late 1960s, began con- 
sidering 4 7ientific modules to dock with space stations. By 
1969, European interest in joint NASA ventures began to con- 
centrate on the "Research and Applications Module" (RAM) , an 

0 

aspect of the American space station studies. During the same 
period, ERNO became involved in studies to enhance Sky lab oper- 

9 

ations, thus adding to EKNO's growing understanding of manned 
missions. By the Autumn of 1969, a series of post-Apollo studies 
resulted in more positive steps for U.S. -Europe joint efforts. 

In Sc|'ton±»er 1969, the Space Task Group advised President Nixon 
to pin :;uc both manned and unmanned efforts. In the case of the 
former, STG recommended a reusable space transport system, and 
called specifically for international cooperation in U.S. space 
programs. In October, NASA Administrator Paine and a select 
staff toured Europe (as well as Australia and Japan) , describing 
NASA proprams and officially inviting their participation. Europe 
quickly reacted to the possibilities, planning a permanent ELDO/ 
ESRO l.inison in Washington, D.C. The European representative 
would meet with NASA officials on a weekly basis and make periodic 
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trips to NASA centers and to contractors.^^ 

By June 1970, ESRO earmarked funds for preliminary ppace 
station module studies, by September 1971, ELDO began collab- 
orating with NASA on shuttle technical studies. In February 
1972, a Joint Technical Experts Group of NASA and the European 
Space Conference endorsed emphasis on "Sortie Lab" studies. 

Urged by the U.S. to increase this aspect of investigation, the 
European Space Council eventually approved a comprehensive sortie 
lab study by ESRO on December 20, 1972. Several considerations 
motivated American and European commitment to the sortie lab, 
eventually re-baptized as the Spacelab. The sortie lab idea, 
originating as a scientific module for space station operations, 
ran afoul of tight budgets. When NASA finally gave up its dreams 
for a permanent space station, the sortie lab became cui attrac- 
tive payload for the Shuttle, providing a flexible and resuable 
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space station in miniature for near-Earth orbit. But NASA 
found itself strapped for funds in 1972, having barely enough to 
continue the Shuttle project!* let alone Spacelab. According to 
Werhher von Braun, former Director of Marshall Space Flight Center, 
and then Deputy Associate Administrator for planning at NASA 
Headquarters, this seemed a highly opportune time to invite 
European participation. Von Braun wrote that Europe seemed eager 
for an active space role, but continued to have major problems 
with KhlV's multi-stage booster, Europ e. In any case, since 
Europe was a "throw-away" booster, and the U.S. Shuttle was re- 
usable, Europe’s large- thrust booster already seemed outmoded. 



9 


OF POOR 

Europ* needed the neans to establish a secure beachhead In manner'l 
spaceflight; NASA needed a moneyed partner; one might say it was 
a marriage made in heaven (or at least in near-Earth orbit) • 

As Europe gravitated towards acqu i^ion of the Spaceled) 
effort, the Continent's organization for space activities became 
realigned and consolidated. In 1968, a special advisory committee 
to the European Space Conference recommended more centralization 
of auUiority in order to facilitate coordinated plans and estab- 
lish a cost-effective approach to European space efforts. The 
ESC would remain as the essential body at the ministerial level, 
with a standing committee to carry out decisions and set policies 
for both R&D and for commercial projects. The new organization 
would replace ESRO, ELDO, and CETS (European Council on Tele- 
communications by Satellite) . Four more years passed until 
the European Space Agency began to assume substantial outlines. 
During 1972, the Communications Satellite Corporation (Cdmsat) , 
dominated by the U.S., continued to block European efforts to 
put up competing satellites. Growing European frustration 
prompted the ESC to take film steps towards organizing the 
Europe.nn Space Agency (ESA).^^ Even so, two more years passed 
before ESA's charter went permanently into force. At an ESRO 
meeting on 15 February 1973, representatives approved a new ESA 
charier, and planned to open it for signatures between March 
and September that year. Final details did not get ironed out 
until 15 April 1975, and ESA became a de facto entity on 31 May 
1975 .^' 
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Against this background, the diplomatic and contractual 
details for Spacelab moved ediead. NASA and European representa- 
tives reached a series of agreements even before ESA became a 
reality. As noted earlier, the European Space Conference en- 
dorsed an ESRO study on the Sortie Lab late in 1972. During 
March, 1973, ESRO issued a statement "concerning execution for 
the Space Laboratory Program," signed by Belgium, Denmark, 

Federal Republic of Gemany, France, Italy, Netherleuids , Spain, 
Switzerland, United Kingdom, and later, Austria. This state- 
ment involved financial commitments and budgetary estimates. 

In many ways, it can be considered the beginning of the Spacelab 
program in Europe. During the Spring, contractor definition and 
cost studies proceeded; ESA began to form, and on 24 September 
1973, European and American representatives signed a Memorandum 

of Understanding (MOU) that established the respective roles of 

18 

Europe and the U.S. in the Spacelab program. 

Obviously, Europe did not plunge blindly into Spacelab 
development. Paper studies and negotiations with ESRO and ESC 
extended over many years prior to the MOU of 1973. In addition, 
European space firms had already received an impression of how 
NASA and its contractors approached complex manned space projects. 
The U.S, had its own sortie lab studies, and European firms were 
already involved in them, 

TECHNICAL INTERCHANGES 

Within NASA, Marshall Space Flight Center (MSFC) seized an 
early position in prototype "Spacelab" studies. In the early 
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days of Marshall's investigation# planners referred to the pro- 
ject as the "Sortie Can." JIC# with principal responsibilities 
for the Space Shuttle# received a major MSFC briefing on the 
Sortie Can (as a likely payload} in April 1972. "As most of 
you know#" the JSC coordinator wrote to an invited audience# ' 
"the Sortie Can has received a considerable amount of attention 
by Marshall and NASA Headquarters# and this will be our first 

occasion to hear all the details directly. concept also 
used the acronym# RAM# meaning Research and Application Module. 

MSFC's activities included an extensive contract with General 
Dyneunics/Astronautics for RAM studies# which "analysed the Pay- 
load/Shuttle interface to a considerable depth. General 
Dynamics submitted a lengthy document on 12 May 1972. The study 
itself holds interest as an early, comprehensive Spacelab proto- 
type overview; as a clue to early European and U.S. contractor 
collaboration; and for description and illustration of early 
module studies. 

Although General Dynamics led the contractor team# support 
caune North American Rockwell Space Division (later called Rock- 
well International) , Bendix# and TRW Systems Group. The latter 
eventually wound up as a sub-contractor to the European effort. 

In addition# General Dynamics had international assistance from 
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SAAB-Scania# Kesserschmitt-Bolkhow Blohm# and from ERNO. The 
latter two both became competitors on the Spacelab project. In 
fact# EHNO's Spacelab studies# and its competition with MBB# 
predated its U.S. collaboration. In March of the 


same year# a 
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European conecrtium known as MESH launched a Spacelab feaal^ 
bllity study* with ERNO as the project leader. Competing pro- 
posals came from the "Cosmos” consortium (headed by MBB) and 
from "Star” (headed by British Aerospace Corporation.) ERNO 
and MBB thus accumulated considerable capability concerning 
Spacelab uses and design details. The General Dynamics report In 
which they collaborated projected several RAM variemts for dif- 
ferent missions. There were two basic mission profiles: RAMS 

carried Into orbit and still attached to the orbiter for seven- 
day missions; several RAMS delivered Into orbit and joined to- 
gether to create a modular space station. The individual RAM 
units could be varied for different missions. A pressurized 
sortie RAM* attached to the Shuttle Orbiter* provided a basic 
space laboratory with shirt-sleeve environment for mission 
specialists ;^RAM support module provided more room for addition 
personnel; the RAM payload module was a "stripped-down" version. 
Although pressurized* It was basically a cargo space and had a 
minimum of autonomous syfcsy stems, depending on the Orbiter or 
another RAM unit to function. The study also discussed RZ^ pal- 
lets* semi-circular structures with control and data hook-ups, 

mounted in the open section of the Orbiter *s payload bay to 

23 

expose experiments to the space environment. This detailed re- 
view possessed striking similarities to eventual Spacelab modules 
and p.nllet configurations. European participation* including MBB 
as well as ERNO* represented an unusually early and extensive 
stage of collaboration* well before Europe 
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■ and America signed the formal MOU in the Autumn of 1973. It 
would also seem that the experience provided early opportunities 
for the eventual association between the European Spacelab con- 
tractors and U.S. firms like TR^7 and Rockwell, 

MSFC, in the meantime, pursued work on the project, which 
received an official change in designation, from "Sortie Lak." 
to "Sortie- Ceui. " MSFC's effort included fairly detailed studies 
for laboratory equipment inside the Sortie Can, stressing off- 
the-shelf material as much as possible. Marshall technicians 
also established guidelines for identifying flammable and toxic 
materials, temperature and pressure environments, and the acous- 

0 A 

tic environment. Planning intensified during the summer of 
1972, as preparations for a NASA-wide "Shuttle Sortie Workshop" 
began to take shape. The workshop held at Goddard Space Flight 
Center from 31 July to 4 August engaged personnel from Head- 
quarters and Centers alike, and the agenda convered a broad spec- 

25 

trum of management issues and technical questions. The Shuttle 
Sortie Workshop also highlighted controversy about the size and 
utilization of modules. One camp advocated a well-equipped lab 
modulo, supporting two or more specialists for seven days. A 
second group argued for a minimum crew, filling the module with 
automated research equipment and allowing missions beyond seven 
days.*" Some observers felt the whole meeting was too elaborate. 
"For mo, the entire week was far too much time to spend, relative 
to the results," one JSC representative reported. "I believe 
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the equivalent could have been accon^llahed in about 2 days.* 

The JSC correspondent also reported MSFC's desire to build 
the Sortie Can in-house to save money. *It sounded like the 
clock had been turned back several years when HSFC went from 
the wet workshop to the dry workshop,” he wrote, referring to 
Sky lab. 

Justified or not, the remark suggests a logical consequence 
of Spacelab's heritage. MSFC spent considerable time and effort 
on the Sky lab program, and its expertise obviously provided 
reference points in Spacelab's evolution. European firms must 
have picked up useful information from their Skylab association, 
as well as firm participation in RAM investigations. At the 
S 2 ime time, a succession of NASA activities strongly influenced 
Spacelab's basic configuration. Many of the preliminary design 
problems nnd technical issues became ironed out as a result of 
contract studies like the General Dyncimics RAM study, continuing 
in-house work at MSFC, and agency-wide workshops such as the one 
at GSFC. 

As Spacelab's outlines emerged into sharper focus and NASA's 
negotiations involving European participation neared formal 
agreement during 1973, exchanges at the technical levels were 
already well advanced. By the Spring of 19 73, MSFC received the 
designation as "lead center" for the Sortie Lab. JSC questioned 
the decision, and Dale Myers, Associate Administrator for Manned 
Space Flight reiterated the decision in a letter to Christopher 
Kraft. Even though Marshall remained the lead center, many 
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issues involved the Shuttle 0rblter,and'corre8pondence often 
beceuno a three-way affair, involving JSC, MSFC, and European 
space officials. As the technical discussions focussed on in~ 
creasingly specific details, the first areas of technical dif- 
ferences began to crop up. 

Constituted more as an international association, ESRO 
originally relied on contractors within its member states for 
technical contributions. But the organization recognized the 
usefulness of an indigenous and independent facility with tech- 
nical expertise to perform certain in-house engineering pro- 
jects and to enhance management of assorted contractors. For 
this reason, ESRO established the European Space Research and 
Technology Center (ESTEC) at Noordwijk^ near Amsterdam, the 
Netherlands. With a staff of several hundred, ESTEC developed 
facilities for research, development, and testing of space 
hardware and space vehicles. At this stage of Spacelab's 
evolution, ESTEC personnel became more enmeshed in technical 
issues. In May 1973, Heinz Stoewer of ESTEC voiced concerns 
to Bob Thompson of JSC, following a visit to Houston by a Space- 
lab team. Stoewer worried about the center of gravity and the 
landing weight of the Shuttle/Specelab combination. After news 
of possible changes in center of gravity constraints popped up, 
"there has been considerable activity, as you know, on both sides 
of the Avi.mtic," Stoewer wrote, "and the results are far from 
encouraging." The arrangement preferred by ESTEC kept the Space- 
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lab snugged up to the Orblter's forward bulkhead. Mounting the 
Spacelab further aft meant a longer access tunnel, more weight, 
more crew transfer problems, and less space for using pallets, 
Stoewer pointed out. He also suggested more flexibility for 
landing weight, targeted at 20,000 lbs. empty. If weight in- 
creased, it meant reducing payload or redesigning a smaller 
Spacelab.^® Thompson wrote back, alluding to many issues alea 
ready settled, and referring to documentation underway to en- 
hance the Spacelab effort in Europe. As for the center of 
gravity issue, Thompson said that the first two Orbiters would 
have to stay at 2% envelopes; a modification kit for Orbiter 13 
would give ESTEC the 3% they wanted. As for the landing weight, 

the payload limits had been raised, allowing for both Spacelab 

31 

growth and payload capability as well. 

In addition to periodic visits andcontinuing correspondence, 
personnel assignments to European sites assisted ESRO/ESTEC in 
evaluating contractors and in negotiating the Memorandum of 
Understanding, pending in the Autumn of 1973. Eui.opean design 
of the Spacelab made good progress, but they Slso recognized the 
wisdom in tapping accumulated experience of NASA experts. By 
the Spring of 1973, NASA already stationed two liaison officers 
at ESTEC, with additional people ^n route during June. As Myers 
explained the situation to the NASA Administrator: "In response 

to a si'ccific ESRO request for technical assistance during the 
next 4-6 w'ceks, we are detailing to ESTEC from our Centers six 
technical experts in the fields of systems engineering, payload 
requirements, environmental control, safety, flight operations. 
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and ground operations." Their job called for working with the 

E5R0 team looking at European contractors and evaluating their 

capabilities. ESTEC also placed liaison person in residence at 
32 

MSEC. NASA Headquarters evidently acquiesced in the dispatch 
of NASA experts to evaluate European contractors because the 
ESRO request appeared to be unique. As ESTEC built up its tech~ 
nical staff# there would be less demand for U.S. technical sup- 
port over a lengthy period. On the other hand, NASA foresaw the 
need to send large observation teams to Europe as ESTEC reached 

certain milestones for systems and subsystems reviews, as well 

33 

as for other key events. 

In fact, a major NASA deputation attended a "Phase B2" 
review during the Summer covering contractor definition and cost 
studies. Including personnel from Headquarters, JSC, MSEC, and 
KSC, 22 NASA executives and engineers assisted their ESTEC counter- 

I^arts in evaluating the competing presentations for the Spacelab 

/ 

/bontract. MBB and ERNO each made their pitches at this time, 
backed by respective subcontractor teams from other nations. The 
NASA emissaries came away impressed by what they saw and heard. 
"Beth contractor briefings were very professional and related 
the in-depth technical abilities of the entire team," a JSC 
delegate reported. Eventually, ESTEC and NASA teams both sum- 
marised th'>ir findings for ESRO management. The information from 

I * 

these deliberations became part of the information used by 
European negotiators for structuring the Memorandum of Under- 
standing in September, and for the selection of the Spacelaib 
contractor by Spring of the following year. ^ 


These events not 
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only underscored the collaboretion in assessing many Issues be- 
fore signing the Memorandum of Understanding; NASA's action in 
sending experts from the various Centers to review the Europeeui 
contractors also revealed an unusual degree of participation in 
the process of selecting the eventual Spacelab builder. In 
addition, NASA participation gave the agency, an early look at 
the contractors, and presum£d3ly facilitated communications dur- 
ing subsequent phases of Spacelab development. 

As this round of contractual arrangements neared conclusion, 
other international liaison efforts focussed on the probable 
scientific users of Spacel€d> and attempted to pl 2 ui the lab's 
equipment and functional layout. At JSC, the Life Sciences 
Division played an important role in coordinating these activitiej 
receiving continuous data from within NASA and from potential 
customers outside the agency. The Life Sciences Division took 
a special look at such things as: dimensions; power; heat 're- 

jection; computer services; communications; rescue, EVA; crew 
size; operational modes and responsibilities. The Europeans 
had to consider customers from their side and work out appro- 
priate arrangements with NASA, since Spacelab functions inevit- 
ably involved capabilities and systems aboard the Shuttle Orbiter. 
As a clearing house, both sides agreed on the need for a Joint 
NASA/ESKO Spacelab Users Requirements Group, which held its first 
meeting on 17 July 1973. And all the while, the minutiae of a 
major international manned space program, in which the Europeeui 
contribution rode as a passenger aboard a NASA vehicle, generated 
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a growing steam of information requests from the Continent. 
Queries ran the gamut, from miscellaneous interface issues 
(electrical, environmental, thermal, and computers) to routine 
information exchanges about the qualification status of elec- 
trical parts.®® 

L:ite in August 1973, the UOU draft began to circulate, of- 
ficially entitled: 

Memorandum of Understanding Between the National Aero- 
nautics and Space Administration and the European Space Re- 
search Organization for A C(K>perative Program Concerning 
Development, Procurement, and Use of a Space Laboratory In 
Conjunction With the Space Shuttle System. 

As the document spelled out responsibilities, ESRO agreed to 
develop "mannable laboratory modules and unpressurized instru- 
ment platforms (pallets)." For ESRO, this mean one flight 
unit, one engineering model, ground support equipment, spares, 
drawings, and so on. NASA, for its part, agreed to maintain 
liali^on personnel in Europe, and to accomodate European personnel 
posted to the U.S. NASA further agreed to provide general tech- 
nical and managerial consultation, and to "monitor ESRO tech- 
nical progress in selected areas as defined in the Prog'^'amae 
plan." In this regard, NASA made sure that crucial aspects of 
Spacelab development met its rigorous standards, which had, after 
all, become honed on successive generations of manned spaceflight 
operations, from Mercury to Gemini, Apollo, Skylab, and the 
Apollo-Soyu? Test Project. As partners in a new enterprise, the 
Europeans also picked up valuable engineering and management 
expertise from NASA. 
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In fact# potential European benefits were specifically 
singled out# permitting acquisition of U.S. aerospace technology# 
As phrased in the MOD# NASA would not infringe any U.S. pro- 
prietary rights# but did agree that ’’ESRO will have access to 
technology# including know-how# available to NASA and needed to 
accomplish successfully its tasks under this cooperative 
programme. • .* As a quid pro quo # NASA received the seune 
prerogatives regarding ESRO activities. ESRO agreed to retain 
a capability to manufacture additional Spacelab units# and main- 
tained a monopoly in this regard# since NASA agreed not to in- 
augurate any independent Spacelab development, unless ESRO failed 
to produce. Furthermore# Europe secured an additional position 
in manned space flights with the concurrence of NASA to include 
at least one European in the first Spacelab flight crew. The 
MOU also outlined an important mechanism to coordinate important 
mcinagcrial and technical issues through the Joint Spacelab Work- 
ing Group (JSWG) # co-chaired by U.S. and European Program heads# 
assisted by appropriate technical personnel. The JSWG exchanged 
essential information concerning the Orbiter and Spacelab; con- 
tended with interface problems and answers; identified potential 
problems for each other; made sure corrective action occurred. 

Formal signing of the Memorandum of Understanding occurred 
24 September 1973 in Washington# D.C.; Dr. Alexander Hocker of 
Germany .«ioned for ESRO# and James Fletcher# NASA Administrator 
signev! Tor tlie U.S. The same day# Fletcher sent a NASA memorandum 
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officlRlly re-cbrlstenlng the Sortie Lab as the Spacelab. 

The Washington Post noted prior European failures to develop a 

reliable launcher; now the Europeans had embarked on a more coio- 

plex manned space project costing $300-400 million. Since the 

European Spacelab constituted a major element in future U.S. 

space programs, the collaboration represented "a new dimension 

38 

in the Atlantic partnership." Other news commentators recalled 

Europe's failure in mounting a significant space program, and 

seemed dubious about the on-going relationship implicit in the 

MOU. Hie U.S. remained dependent on Europe for Spacelab units, 

and this did not necessarily seem to be a good thing. In any cas< 

Europe appeared to have a firm monopoly on Spacelab manufacturing 

39 

as the only possible supplier in the near future. If nothing 

else, Europe had definitely found a secure niche in manned space 

operations for many years ahead. 

The contract award for Spacelab finally went to ERNO on 5 
40 

June 1074. The European space reorganization begun ducihg a 
meeting at Brussels on 31 July 1973, eventually resulted in the 
ll-niv>mbt'r European Space Agency, which became official on 31 May 
1975.'^^ Members included Belgium, Denmark, Federal Republic of 
Germany, France, Ireland, Italy, The Netherlands, Spain, Sweden, 
Swit:*or land, and the United Kingdom. These events put Europe’s 
mannod s'pace program on a firm basis. As Spacelab development 
moved alioad on a course parallel to NASA's Shuttle program, the 
"Atlantic partnership" looked forward to new dimensions in space 
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Gooding (^978) suggests that, in Uio oxidative weathering of primary iron 


sulfides (such as troilite, FeS, and pyrrhotite, Fe.. .S) under Mars surface 

\ 1 ~X ) 

ambient conditions, the most stable product would be FcSO|^ or, given sufficient 

-7 

water vapor (2.5 x 10 atm at 240 K) , its monohydrate szomolnokite , FeSO^^ H^O: 
FeS + 20^ FeSO^: AG^^gg = - 174.08 kcal ( <) 

FeSO^ + H^OCg) FeSO^ * = ‘ 

but that the equilibrium O 2 pressure for the reaction 

2FeS + 3 / 2 O 2 Fe^O^ + 2S (3) 

is nearly the same as that for reaction (1), so that reaction ( 3 ) or the even more 
exoergic 

2FeS + 7/20^ Fe^O^ + 2 SO 2 ; -^G°298 ^ 

might provide alternative product sets. The oxidation of iron sulfides would 

be an important provenance for sulfate in the Martian regolith, and therefore, 

hypothetically (Toulmin et al . , 1977; Clark et al., 1976) for the cementing agent 

of duricrust. The assemblage (Fe^O^ + 2S + 5/20^) is about 220 kcal higher in 

Gibbs- free energy than 2FeSOj, at 240 K and the ambient Martian oxygen fugacity of 
“5 

1 X 10 , so that once FeSO^^ is formed there is no possibility of its conversion 

to hematite and elemental sulfur under these conditions. Finally, the combined 

oxidation and carbonation reaction 

2FeS + 1/20_ + CO ^ FeCO + + S; AG° = - 20.1 kcal (5) 

.j 298 

would lead to spontaneous formation of siderite- 

However, in kinetically hindered reactions it is the relative rates of 
alternative pathways and not thermodynamics alone that determines tiie product set 
that is first produced, and if hematite and sulfur were once formed they could 
persist raetastably for long times under a strictly gas-solid kinetic regime. The 
stabilities of FeCO and FeOO. would depcr,d on tli--' availabllily of liquid (or 
other condensed, e.g., adsorbed) M^O to serve as a medium for further rejctions. 


The present work w.is undertaken to investigate the surface oxidation of a set 
of primary iron sulfides (troilite, pyrrhotite, and pyrite, under Mars 

surface ambient conditions, and in the presence of a strong flux of ultraviolet 
light, in the assumption that attack by photochemical! y generated oxidants such as 
0 ( 'd) , 0^, and (in the presence of wat«^r vapor) HO^ and ^^ 2^2 dominate gas- 

solid oxidation kinetics. After irradiations under a variety of conditions of 
temperature and humidity, the surface chemical state of sample chips as measured 
by the X-ray photoelectron spectroscopy (XPS) of Fe, S, C, and 0 was compared to 
that of unexposed reference chips. Oxidation of sulfur from -2 (sulfide) to +6 
(sulfate) states is accompanied by an easily detected positive shift of about 6 to 
7 eV in its 2p electron binding energy. Relative iron and sulfur intensities give 
evidence of the extent to which sulfur anions are replaced by oxide or carbonate 
in the weathering of these minerals, while the presence of carbonate and iron 
oxides is detectable from characteristic shifts in the C(1s) and 0;1s) electron 
binding energies. 

EXPERIMENTAL 

Samples were placed in A1 containers on a Cu heat sink within a 3*7 liter 
stainless steel chamber (Fig. 1). Sample temperature and chamber humidity were 
crntrolled by the circulation of thermostatted coolant through a coil beneath the 
heat sink; temperature was measured by stainless steel jacketed thermocouples at 
three positions: (1) At the point on the heat sink where the coolant coll enters; 

2 ) at the surface of a sample chip, in the UV beam; and (3) on the chamber 
baseplate, or immersed in the liquid water reservoir when present. 

Light from a 50~watL deuterium lamp Wras admitted to the chamber through a 

0.3 cm thick optical quartz window, alter traversing an air path of 8 cm. The 

emission spectrum of this source at the quartz window i: shown in Figure 2. The 

3 

volume of the beam inside the chamber was 50 cm . 


Polycrystalllne bulk mineral samploo ( id<^ntlty verified by powder X-ray 
diffraction) were sawn with a diamond blade, using purified isopropanol as 
lubricant, and cleaned by ul trasonicat ion in high-purity trifluorotrichlorocthane, 
followed by a bakeout. The face to be irradiated was then polished in an ftr 

atmosphere with SiC paper to produce a fresh surface. Samples remained in an Ar 
or CO^ atmosphere from the termination of irradiation until introduction into the 
sample chamber of the X-ray photoelectron spectrometer. 

The irradiation atmosphere was ”high-purity” carbon dioxide containing 
350 ppm 0^, except in the lowest-humid ity run, when CO^ with certified and 0^ 
concentrations of 20 ppm was used. In all but this latter run, the H20(g) compo- 
sition of the gas phase present was controlled bv varying the temperature of the 
heat sink, in the presence of a supply of H^O (1 r s) . The variation of this 
parameter is summarized in Table I. The atmospheric leak rate of the assembled 
apparatus was about 8 x 10 torr/day, and thus the total atmospheric oxygen 
introduced through leaks during a 100-hour irradiation did not exceed the partial 
pressure of 0^ (0.13J of ambient pressure) at the surface of Mars. 

X-ray photoelectron speutroxopy (XPS) was performed at Texas A&M University 
on a Hewlett-Packard 5950 photoelectron speqtrometor using monochromatized 
radiation. Data reduction employed the Surface Science Laboratories ESCA Data 
System, Revision E. 

RESULTS 

At the time of this report, the analysis of only one sample set (cold, humid; 

see Table I) has been completed. The wide-scan XPS spectra (Figure 2) show 

prominent peaks for carbon and oxygen in addition to those for iron and sulfur, 

indicating considerable alteration of the surface from the stoichiometric bulk 

compositions. Hign resolution scans of the sulfur (2p) , carbon (Is), and oxygen 

2 - 2 - 

(Is) indicate the presence of sulfur as S and ' ; carbon as C(IV) (either CO^ 


or CO^ (ads) and other compounds possibly includinc COG and (Blacl:burn et al . , 
1Q80))and, except in the case of oxygen as a metal (probably iron) oxide. 

Thus, an assemblage of FeSO^, FeCO.^, and original mineral substrate 

would account qualitatively for the appearance of the XPS spectra. The assumption 
that only these substances are present will permit calculation of the mole fraction! 
of each present in the layer (about ? twn deep) analyzed by XPS, and such calcu- 
lations will be performed automatically as data sets are recorded. The results of 
these calculations will then be correlated with conditions during irradiation (Table 
It is emphasized that such a correlation will have only conjectural significance, 
since we have no firm evidence that the four components listed above are the only 
ones present in significant amounts at the irradiated surface; indeed, the XPS 
spectra not yet obtained may lead to a revision of the basis set of compounds 
assumed. (For example, there are unassigned peaks in the C(1s) spectra of most 
samples . ) 


Table I. 
Run 

1 

2 

3 

5 


Environmental Parameters During Irradiation. 


HgO 

^Sample 

torr 

^CO^, torr 

Rubric 

246 

249 

-4 

1.6 X 10 

8 

Cold, dry 

249 

256 

0.53 

8 

Cold , humid 

305 

309 

35.6 

40 

Warm, humid 

300 

297 

26.7 

38 

Warm, wet 


233 


760 

Cold, dry, dark 


(No UV irradiation! 


Notes: (a) Temperature of the coole.st point in the irradiation chamber. 

(b) Temperature of in-beam thermocouple touching sample chip. 

(c) As referenced in the text of this paper. 
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CAPTIONS FOR FIGURES 


Figure 1. 


Figure 3. 


Irradiation chamber. All construction is of stainless steel except for 
the Cu gasket and heat sink (S) and quartz window. 

B: UV beam volume. 


TC: Thermocouple leads, 


Figure 2. Emission spectrum of deuterium arc lamp. 




Typi cal post-irradiation XPS spectra, Fe^^ ^^S, out-of-beam sample. 


A: Wide-scan spectrum sliowing peaks for Fe (2p), 0 (Is), C (Is), and S 
( 2 s and 2 p) electrons. 

B: High-resolution 0 (Is) spectrum. The low-binding energy shoulder at 
530 eV indicates iron oxide(s) present. 


C: High-resolution C (Is) spectrum. The peaks at 288 and 28^.4 eV are 


2 - 


assigned to CO^ and hydrocarbon, respectively; the 286 eV component 


is not identified. 


D: High-resolution S (2p) spectrum. All peaks are doubled because of 

2 

the multiplicity of the P , ,5 state produced by photoelectron 

emission. The doublet at about 169 eV iniicates sulfur as the 


lower-binding energy doublets indicate tv.’o chemical states of sulfur. 


2 - 2 - 
probdbly S ana 
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Abstract: 

Several Antarctic soils were studied. Water extracts were analyzed for 
sodium, potassium, magnesium, calcium, chloride, and nitrate ions using flame 
photometry, atofnic absorption spectroscopy, and ion selective electrodes. 
Significant differences were noted in soils from different regions. 'Jpwar'^ 
ionic transport above the permafrost layer with almost no evidence of migration 
below this line was observed in core samples. Very high ionic concentrations 
of sodium and chloride were observed in the Prospect Mesa Formation of Wright 
Bailey. Calcium chloride was the most predominant species in soils from the 
Don Juan Pond. The high concentrations of salts reflects the aridity of the 
dry valley system. 

The dry valleys of Antarctica provide an excellent terrestrial analog of 
the surface of Mars. Studies such as this one provide important information 
to aid in understanding the processes operating >.ri.thin the Martian regolith. 
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DETERMINATION OF WATER-SOLUBLE IONS 
IN SOILS FROM THE DRY VALLEYS OF ANTARCTICA 

The dry valleys of Antarctica furnish one /f the best terrestrial analogs 
for the surface of Mars. In many places surface features are similar. Both 
Antarctica and Mars have low humidities, desiccating winds, and an oxidizing 
environment. Even though the air temperatures of both are seldom above freezing 
even in summer, the soil experiences a diurnal cycle of freezing and thawing 
(1). It is possible to increase our understanding of changes occurring within 
the Martian regolith by studying processes operating in the cold, dry environ- 
ment of Antarctica. 

During the austral summer of 1979-80, a group of Antarctic soils was 
collected. The primary purpose of the project described here was to determine 
the abundances of the water-soluble ions magnesium, calcium, potassium, sodium, 
chloride, and nitrate in soils from different regions and at various locations 
and depths within a given region. These data not only give soil composition 
but also provide information which is useful in studying processes such as 
weathering and ionic transport. When combined with additional studies of the 
soils, they provide a detailed picture of the soil chemistry in a cold, dry 
environment. 

To retard ionic migration, soil cores were stored below -10®C, The cores 
were cut into individual samples which were then air dried, weighed, and extracted 
with water. After filtration, the samples were analyzed. 

An Orion Research lonalyzer 901 and the approfrlate ion selective electrode- s 
were used for determining calcium, chloride, and nitrate in all samples. Sodium 
and potassium were determined by ion selective electrodes for all samples except 
those in the four detailed core studies (DJ 33» DJ 207^. WV 52, and DQ 35). Fo^' 
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those samples, sodium and potassium were determined by flame photometry using an 
Instrimientation Laboratory 6^43 flame photometer. A Pei*kin-Elmer 373 atomic 
absorption spectrophotometer was used for all the magnesium determinations. 

Because of the aridity of the Antarctic environment, water-soluble salts have 
not been leached out of the soils and are often present in high concentrations. 
Salts accumulate within rocks by natural desert processes, are transferred from 
rocks to the soil during weathering, and by the action of .4nd and by slow partial 
leaching eventually accumulate in low-lying areas such as basins and lakes. 

fiost of the samples for this study were obtained from Taylor Valley and 
Wright Valley, two ice-free valleys in Victoria Land (see Figure 1). Taylor 
Valley runs east-west and is bovnded on the west by the Taylor Glacier and on the 
east by McMurdo Sound. Hills and mountains on the north and south separate Taylor 
Valley from the rest of Victoria Land. Taylor Valley serves as a basin for a much 
smaller region than that around Wright Valley. This was reflected in the data 
which showed much lower ionic abundances in Taylor Valley soils than in the other 
soils. 

The salt layer in Lake Hoare, a dry lake bed in Taylor Valley, was very thin, 
and concentrations were low. The permafrost level was near the surface. Even 
greater ionic depletion was observed in surface soil from an island in the lake. 
Figure 2 shows variations with depth in a core from Taylor Valley. Note from the 
scale of the graph that all concentrations were low. 

Patterned ground consisting of a system of polygons is common in many areas 
of Antarctica, including Taylor Valley. These features apparently result from 
seasonal expansion ard contraction. During the winter, fractures occur as the 
permafrost contracts due to the cold temperatures. By summer these cracks may 
become filled as wind and melting ice carry surface material into the fracture. 
When the soil expands during the summer, the weaker zone at the fracture is pushed 


upward. Repetition of this process through the years causes the observed patterns 
(3). Several areas of Mars show similar patterns except that the scale is very 
different (4). A sample of Taylor Valley soil taken between the polygons was low 
in all ions checked. If Martian polygons were formed by the same mechanism 
responsible for the patterned ground in Antarctica, it is likely that the Martian 
surface soil between the polygons would show similar ionic depletion. 

Wright Valley is one of the major valleys in Victoria Land. It reaches from 
the Polar Plateau to the Wilson Piedmont Glacier where it is filled by the Wright 
Lower Glacier. It is bounded on both north and south by mountain ranges. Tie 
western section of the valley is separated by the Dais Plateau into the upper and 
lower forks. Concentrations of ions in soils from Wright Valley were much higher 
than those in Taylor Valley soils but varied dramatically with depth and from one 
location to another. Figure 3 is a comparison of chloride ion concentrations in 
four different cores from Wright Valley. 

Two Wright Valley ponds gave widely divergent results. The Don Quixote Pond 
in the North Fork had relatively high concentrations of most ions near the 
surface, with sodium and chloride ions predominating (see Figure 4). Concentra- 
tions decreased dramatically to a depth of approximately 3 to 4 centimeters, where 
they started to level off to low values changing very little throughout the rest 
of the core. Upon visual examination, the core appeared homogeneous below 
approximately 3 centimeters. Analysis of the white crystalline encrustations on 
the surface of the pond gave results which were very similar to those from the top 
of the core. All ions checked, except potassium, were present in relatively high 
amounts, indicating the presence of several differenc minerals. The data showec 
that salts from upper-lying areas were transported into lower areas such as ponds 
and lakes. As an illustration, ionic concentrations near the surface of Don 
Quixote Pond were more than tenfold greater than the corresponding concentrations 
at a point 20 meters above the pond. 
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The Don Juan Pond in the South Fork remains unfrozen because of its high salt 
content (2). Two core samples were obtained fron. the soil of the Don Juan Pond. 
Sample DJ 207^ was taken from the center of the pond, and sample DJ 33 was 
obtained near the periphery. In the two core samples, ionic concentrations varied 
in a similar manner with depth. However, concentrations in the central core were 
much higher, indicai.ing a movement of the salts to the lowest part of the pond. 

The increase^ abundance of sodium ions near the surface of the central core shows 
that sodium salts are transported more easily than other salts. In both cores the 
preponderance of the sodium ions was near the surface, whereas calcium ion concen- 
tration fluctuated considerably down the core. Figures 5 and 6 illustrate that 
both sodium chloride and calcium chloride are present in the Don Juan samples with 
calcium chloride significantly predominating below the surface layer. The data 
verify that upward ionic migration occurred. As surface water evaporated or 
sublimed, the salts were left behind. From these results, it appears that sodium 
chloride has a greater upward mobility than does calcium chloride. 

Core 52 was obtained from the southeastern end of Wright Valley. Ionic 
concentrations were generally much lower than those of the Don Juan samples. As 
shown in Figure 7, the molar concentration of calcium ions did not correspond to 
any of the other species, predicting the presence of an additional calcium salt. 
Visual examination indicated the likelihood of a small amount of calcium carbonate 
and a much larger amount of calcium sulfate. The broad maximum in calcium ion 
concentration occurring at the 8 centimeter level indicates that the upward 
migration of calcium ions is hindered, probably because the calcium salts present 
in this sample are only slightly soluble in water and are not transported upward 
as readily as more soluble salts. 
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Samples were obtained from a large, 1 meter deep pit the Prospect Mesa 
Formation in Wright Valley. Figure 8 shows the ionic abundances. The molar ratio 
of sodium and chloride ions is very close to 1:1, indicating the presence of 
sodium chloride (see Figure 9). A salt layer occurred just below the surface 
again indicating an upward mobility of sodium and chloride ions. Concentrations 
of these ions decreased down to the permafrost level at a depth of 35 to 40 
centimeters. Below this level, concentrations of all ions remained essentially 
constant, indicating that little or no ionic transport occurred in the permafrost 
zone. The depletion of ions in the surface layer is indicative of weathering 
action, probably due to wind. 

In summary, all samples examined contained water-soluble salts, reflecting 
the aridity of the region. Calcium chloride was most abundant in samples from the 
Don Juan Pond, whereas sodium chloride was present in very high concentrations in 
the Prospect Mesa Formation of Wright Valley. Considering all ions checked, the 
soil pit in the Prospect Mesa Formation was the richest area examined. Abundances 
of both potassium and nitrate ions were low in most areas but were significant in 
some of the samples from the Prospect Mesa pit. Magnesium ion concentrations were 
generally low throughout the dry valleys. A few samples were checked for fluoride 
ions; the results were so low that further determinations seemed impractical. 

Movement of salts into low-lying areas was verified. Upward ionic migration 
was evident in all core samples. Of all cations observed, sodium showed the 
greatest degree of migration. 

The chonistry of soils from the dry valleys of Antarctica is both interesting 
and unusual. A detailed knowledge of these soils may prove to be quite valuable 


in planetary studies. 
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I. Introduction 


The Microwave Scanning Beam Landing System (MSBLS) is a navigation 
system which provides distance, azimuth, and elevation coordinates to the 
Shuttle orbiter during final approach and landing (Harton 1979), Azimuth 
information is transmitted to the orbiter via a pulse-encoded, vertically 
polarized, planar vertical beam sweeping through a 30® transmission sector 
centered on the runway centerline extension; vertical coverage of the azimuth 
beam extends to 23®. Similarly, elevation data are transmitted to the orbiter 
via a pulse-encoded, vertically polarized, planar horizontal beam sweeping 
trhough a 30® transmission sector; usable azimuth coverage of the elevation 
beam extfends to approximately 50®. For range determinatioi-. the orbiter, 
after detecting a solicit signal from the ground station, returns an in- 
terrogation signal. Upon receiving the interrogation signal, the ground 
station sends a reply signal to the orbiter. The orbiter then determines its 
range by measuring the time spacing between interrogation and reply. The 
entire sequence, range-azimuth-elevation, is repeated every 200 milliseconds, 
and at a frequency near 15.5 GHz. (wavelength 1,94 cm). 

At a fixed tracking level of about -4 dB, the azimuth full beam width 
is -v* 2,4? Since the ground station changes the pulse encoding (pulse-pair 
spacing) every 0.125®, as the beam sweeps past the orbiter, a distribution of 
pulse pairs is received. In computing azimuth, the orbiter averages the pulse- 
pair spacings of all pulse-pairs in the received distribution. Elevation 
signal reception and decoding is similar to that for azimuth, except that the 
-4dB beam width is 1.4°, and the received pulse distribution contains fewer 


The HSBLS is expected to be particularly useful during periods of de- 
creased optical visibility such as at night or those resulting fron the 
presence of rain. fog. or clouds in the line-of-sight between the ground 
station and arbiter. Certainly, there will be times when meteorological 
conditions are sufficiently poor that an arbiter landing would be un- 
acceptably risky. Consequently, it is important to be able to predict 
the signal attenuation, beam distortion, and coordinate errors which re- 
sult from operation of the system in poor weather conditions. 

Radar literature abounds with studies of microwave attenuation due to 
rain and fog. Robel (1978. and references therein; 1981) has appealed to 
those studies in estimating the most adverse meteorological conditions under 
which the MSBLS will function reliably. He deals with the problem of 
attenuation exclusively. One can anticipate, however, rain and fog conditions 
that will lead to systematic elevation and azimuth error as well as the 
random error expected from a weakened signal at the arbiter. The systematic, 
or bias error results from beam distortion as it traverses asymmetrical 
weather configurations such as those which could occur in fog. rain or cloud 

conditions. 

The main physical processes giving rise to attenuation of MSBLS 
signals are absorption and scattering of microwaves from distributions of 
water droplets. In scattering. , the frequency of radiation after encounter 
with a drop is the same as the frequency before the encounter; only the 
direction of propagation changes. In .absorption, the microwave energy goes 
into heating the drop, which later radiates isotropically with a black- 
body (Plankian) distribution of 


frequencies. Since scattering is a narkedly directional phenomenon, one 
expects beam shape distortion in scanning bean systems like the MSBLS to 
depend on the relative importance of scattering and absorption. These two 
effects must be evaluated if one is to determine whether or not beam distortion 

is an important source of error. 

In proceding. we utilize a fully general theory of scattering and absorption 
(Hie 1908; Debye 1909), valid for all drop sizes from fine mist to heavy rain. 
Availability of high speed, high capacity computing machines renders this 
approach preferable to those involving various limiting approximations for 
drops of different sizes. Since the general theory lends itself quite readily 
to attenuation calculations, we use them as much as possible as a check to 

a 

ensure that the progracroing Is correct. Necessary formulae are collected 
In Section II. Attenuation, scattering and absorption efficiencies are 
discussed In Section III, and finally, a general discussion and recommen- 
dations appear In Section IV. A listing of the computer program appears in 

the ap» ‘ - iix. 


II. Formulae 

Mie theory gives, among other things, the far-tield solution to the 
problem of scattering and absorption of plane-polarized E&M waves of sinu- 
soidal time variation. VII thin the Mie formalism, one can determine in- 
tensity and polarization of microwaves scattered from a single spheroidal 
drop; the generalization to drop-size distributions is straightfonvard 
We do not wish to reproduce the theory herein; interested reao . s are 
referred to Van der Hulst (1957). Our preference is to extract only those 
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formulae necessary for our computations: 

n 

i a: radius of sphere 

x: wavelength of microwaves 

1: real part of complex refraction index, m 

k*: imaginary part of complex refractive index, m 





X - aiTA/x, 

/n\, » >1 - X ^ . 

Csca* cross section for scattering (units of area) 
Cabs* cross section for absorption 
Cbk 5 cross section for back-scattering 

^sca* efficiency factor for scattering (dimensionless) 
Qabs* efficiency factor for absorption 
Qbk • efficiency factor for back-scattering 
Qext: efficiency factor for extinction (attenuation) 



(3) 


(4) 

f 

(5) 

Qevt ■ Qsca. Q& ^ 

(6) 

*■ >vi>, 

(7) 

^o(-i) = COi ^ y 5,^ ^ 

(8) 




(1) 

( 2 ) 


( 10 ) 
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Jo W* 5 /tv ^ C#l X, (11) 




[: 

l«(V)^ 2v-| 

•HV 4! J 

ReU<’^J 



r 

^IW 

#(« 

L ' 5 <’■> 





^ t ] Re [ (>;] - f?e [ ,, ('»')] (14) 


(Jsc*. = Ib^T] (15) 

'll*-/ ' 

> 1 -/ 

(?b.. ‘ -p i (’7) 

Tf (c «5 &) = o^ (18) 


fij (ccy 6 ) ‘ I, 

■"■> ( C <>5 9J -- 3 cos e ^ 
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Tn^fci e). (C5 6)^ (19) 

Tv fc®* ^) 'O, 

(cej e) r /^ (20) 

2l (ccs ® i c^l 

(coi G) - c©5 sf-nXv (cos e)--rr;^.j^ (c.j (21) 

-(a»i-i)ji»s*^[7i:;,., (fci 6)J+ e)^ 


' where 9 is the scattering angle, measured from the 
direction of the forward beam. 

, 5, (tf) = ^ U«i a) 1 

'H • / J ^ 

A 'T»t4. J 


(22) 

(23) 

(24) 


+ Ra^.b,-] 

/vi 


( 25 ) 


where the asterisks indicate complex conjugate. 


/ 


^ For the case of {/ polarization, i ,e . , the electric field 

vectors are parallel to the plane defined by the incident, and 
scattered wave propagation vectors. 



A =• * ' (26) 


where il’ ^ 

^ * 

r: distance from scattering drop to observer 

I : incident intersity 
0 

I: intersity at r 

S 2 ( 0 ): complex amplitude function (equation 23) 

And for the case of 1 polarization, i .e . , the electric field 
vectors are perpendicular to the plane defined by the incident, and 
scattered wave propagation vectors. 
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( 28 ) 


For unpolarized incident radiation. 


^ j [|s,(«)r4 |sj6jr~ 

I. ‘ k.’-r 


( 29 ) 




• The standard attenuation formula for the case of 
single scattering of drops, each having {the same) radius, a, is 

- A/j, 9 7ra"-<?ex<- 

1(0^ j(°j e. , (30) 


N^: number of drops per unit volume^ 

g: attenuation path length . 

Rain 



We must relate to the measurable quantity p. 

p: rainfall rate (nin/hr) 

D: drop diameter (cm) 

V ; terminal velocity of drop (m/sec) 
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ji- ^'*n/s€c)^ 

= ^7r>/5*AT A/y 3)* 


V can be expressed as a function of D via s polynomial fit to the 
data collected by Medhurst (1965). The specific fit is 


V* (D) =■ 5 ^ , 


v/here 


0 -0.02366 

1 ’ 47.38010 

2 -81.56769 

3 38.93966 

4 26.05586 

5 -21.69620 . 

Combining equations 30, 32, we find 

Qcrl 

xw = xf=;e 


where g is in km. 


Generalization to Size Distributions 


Equation 34 becomes 


^ I o o 




©y*' 

-h 

'D.- 
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where = percent of total volume containing drops of diameter 
between and + AD. 

Specific drop-size distributions used are those of Laws and 
Parsons (1943) as tabulated by Medhurst (1965). In general, rain drop- 
size distributions depend on rain rate ( p ) in such a way that as p 
increases, the distribution becomes wider and is shifted toward larger 
drop size. Figures 1-7 are the distributions and our truncated polynomial 
representations of them. For rain rates within the range 2.5 fc p (rren/hr) 
i 150, linear interpolation between the two distributions which bracket 
the rain rate of interest is done; for values outside of this range, the 
extrema aTone are used. 



Figures 1-7: 


Truncated polynomial representation of rain drop-size 
distributions for rain rates ranging from 150 nri/hr to 2.5 mn/hr. 
See text for details. 
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Figure 4 


Figure 5 







FigurG 6 Figure 7 




Land Fog (Radiation Fog) 


Again, Z(t) ^ T. C , 

but now the measureable quantity which is related to is the 

optical visicility, V. For landfog. 


(30) 


M = 0.00332 V-‘*-533 (36) 

(Koester and Kosowsky 1970), where M is the total water content (gm/m^). 
Converting water content to we find 


X 
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XO') 


-y 


T(o) e ^ 


I su- 


where D = 1.0 X 10”^ cm = average droplet diameter^ 
and V is given in km. 


(37) 


Note that equation 37 has a "canonical” size distribution built into it. 

Sea Fog (Advection Fog) 


Formalism is the same as for landfog, but with 


M = 0.0156V-'' •'•3 


(38) 


(Koester and Kosowsky 1970). Equation 38 is valid for V Z 100 m. 
For lower visibility, one must use equation 36. Equation 30, 38 
combine to give 

. i.-it %io'^ 3 ^ 


(39) 


..3 

where D is now 2 X 10 cm. 

One expects true fogs to have properties intermediate between 
radiation and advection fogs. 




f • 
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Clouds 


The simplest cloud can be viev/ed as high-level radiation fog. 

For the present, we adopt Equation 37 to describe the attenuation of alj_ 
clouds, but not without realizing and appreciating the general complexity 
of true clouds. These systems can contain rain, fog, hail, partially 
melting hail, snow, and partially melting snow. In the future, our 
"cloud'' subroutines could be modified to take account of these possibilities. 


III. Attenuation 

' a) Program Verification 

As a check that major parts of the computer program are working 
properly, sample calculations (at a MSBLS frequency of 15.5 GHz) for 
attenuation due to rain and tv/o types of fog were made. Temperature 
dependent complex refraction indices of Robel (1981) were used in all 
computations. In Figure 8, our results for attenuation y^. rain rate are 
plotted for a range of temperature. 


( 




Figure 8: Attenuation as a function of rain rate. 


These results are in excellent agreement with those of Robel (1978), 
especially since different drop-size distributions and different theoretical 
formalisms were used. Our attenuation - rain rate relationships shov/ a 
slight curvature. This results from the convolution of extinction efficiency 
factors with rain-rate dependent drop-size distributions in the full theory. 

In Figure 9, we plot attenuation v£, optical visibility for tv/o 
kinds of fogs at three temperatures. 




M A 




it 



Figure 9: Attenuation as a function of fog visibility 

for two types of fogs. 


r> 


Agreement between results from the full theory and the Rayleigh scattering 
approximation of Robe! (1981) is fine. 

Finally, we compare our full -theory result for rain attenuation 
with data obtained during static . ts of MSBLS performance In rain (Harton 
and Robel 1980 1980 b). The data have been normalized to units of 

decibel/kilometer, and are plotted in Figure 10 against rain-rate. 
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Figure 1C: Static MSBLS test data. 


Open circles correspond to azimuth beam (4.S7 km path length), while filled 
circles denote elevation (1.03 km path length). Superimposed on the data 
are the theoretical curves of Figure 8. Considering that the data were 
taken in July at KSC, where temperature usually ranges between 25®C and 35®C 
agreement is quite good. 

b) Breakdown into Scattering and Absorption Components. 

Q-values (see Equation 30) for absorption, scattering, and 
extinction (Qext " Qabs + Qsca) are plotted against drop diameter in 
Figures 11-15. 
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One can see that absorption doninates the total attenuation for D £ 1 inn. 
For the drop size range 5^0 (y^)<65 of a typical fog (Koester and Kosowsky 
1970), the ratio Qsca/Qabs is never greater than 2,3X10"2. For rain which 
follows the drop-size distribution of Figures 1-7, the corresponding maxinum 
value is 1.6. Table 1 gives the values of Qsca/Qabs(at T = 30®C) for the 
sizes corresponding to fmax if’ Figures 1-7. 
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Table 1 


Rain Rate 

^^^max^ 

Qsc»/Q«bs (T.30°c) 

150 mm/hr 

0.30 im 

0.28 

100 

0.27 

0.21 

50 

0.23 

0.08 

25 

0.20 

0.06 

12.5 

0.17 

0.06 

5 

0.14 

0.06 

2.5 

0.12 

0.06 


Since for fog, Qsca/Qabs is always less than about 25:1, it is probably un- 

a 

likely that scattering of flSBLS microwave radiation from fog will result in 
serious beam distortion. For rain, the situation is quite different, and 
scattering-induced bean distortion could be a significant source of error. 

IV. Discussion and Recommendations 
a) Fog 

For a given visibility ,advection fog is 3.5 times more effective 
than radiation fog in attenuation of 15.5 GHz microwaves. Except for brief 
afternoon periods, heavy fog is the rule and not the exception during 
summer months at VAFB (Robel 1981). At VAFB, this author has experienced fog 
with optical visibility £ 20 m, A 40 Ion path length through radiation fog 
of 20 m visibility at T=0®C results in 12 dB of MSBLS microv/ave attenuation 
for radiation fog. The 12 dB figure is only a rough lower limit to the true 
attenuation, since it is likely that the fog contains an advection component. 
A priori, there is no way to know what the particular composition of VAFB 
fog is. One can only set limits to attenuation at ^ixed temperature and 



visibility. Depending on (1) the specific nature of VAFB fon. and (2) the 
applicability of complex refractive indices for pure witer, microwave 
attenuation could approach or exceed 40 dB over a 40 km path length. One 
would like to be able to predict signal attenuation and link margins more 
precisely than indicated by our 28 dB range (factor of 630). Consequently, 
it is reconwended that simultaneous measurements of optical visibility and 
MSBLS microwave attenuation be made for VAFB summer fogs of differing 
visibility and temperatures. These data would be used in setting realistic 
zero-points in Figure 9. 


b) Rain 

* Ue have detennined that microwave scattering is non-negl igible for 
intermediate and heavy rains (Section III). Under these conditions, coordinate 
errors due to scatterino-induced beam distortion could contribute significantly 
to the total error budget. It is recommended that the formalism of Section II, 
and the computer program listed in the appendix be used to evaluate the 
specific magnitude o chis source of error. The additional work vs straight- 
forward, and involves the inclusion of beam profile and weather system 
geometry. The scattering phase functions are already included in the program. 
Even though we expect relatively unimportant scattering effects of light 
rain, fog, and clouds, the specific magnitudes of the errors could be 
calculated with the same beam geometry developed for rain. 

c) Closing Note 

This author expects to make considerable progress in beam profile 
inclusion during the remaining two weeks of his summer fellowship. Infon.ial 
notes kept during this tine will be made available to anyone requesting them. 
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Appendix 


Following is 


a listing of the FORTRAN computer program used for 


attenuation and scattering calculations. 
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i: 

2J 

3t 

4: 

5? 

6: 

7t 

8: 

9t 


lo: 
11 : 
12: 
13: 

14: 


is: 

1 

16: 

17: 

800 

is: 

19: 

110 

20: 

120 

21: 


22: 

23: 

24: 

130 

25: 

26: 

27: 

140 

28: 

29: 

30: 

150 

31 : 

32: 

33: 

34: 

210 

35: 

36: 


37: 

38: 

39: 

160 

40: 

41 : 

230 


DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 


42 : 
43: 
44: 
45: 
46: 
47: 
48: 
49: 
5o: 
51 : 
52: 


170 


LMNDFOGf SEMFOGt OR CLOUDS? <0fl»2r3)') 


20 

180 


IMPLICIT COMPLEX(YrZ) 

DIMENSION ZM<30) »ZB<iO) rZMN<30> »ZBN(30) 

P(30) »T(30) »D(36) »ANG(36) 

ZS1<36» 12) tZS2(36» 12) »S<36»12> »SlAS(36f 12) 

S2AS(36f 12) f ANGL(12) 

DL< 12) rXN(12) »QSCA(12) »QEXT ( 12 ) » QBK ( 12 ) » QABS < 12 ) 
DATA DL/.05r . 1 » . 15 » . 2 f .25 » . 3 f . 35 » . 4 f . 45 » « 5? .55 » . 6/ 

1 = 1 

MM=36 ‘ 

PI=3. 14159265 
TP360=2.*PI/360. 

DO 1 M=lrMM 
ANG<M)=10.*<M-1) 

D<M)«=ANG(M):<fTP360 _ _ . 

•CONTINUE ■ 

CONTINUE 

URITE<6»110) 

FORMAT ('’ ENTER NUMEiER OF TERMS (NN) ) 

R£AD(5f 120)NN 
FORMAK ) 

URITE<6r 130) 

FORMAT(' RAINf 
READ(5f 120)NRFC 

WRITE<6» 140) 

FORMAK' ENTER WAVELENGTH (CM)') 

REA1:K5» 120)U 
URITE(6rl50) 

FORMAK' ENTER COMPLEX REFRACTIVE INDEX (A»B)') 

READ <5f 120 )Z1 ... 

URITE(6f 210) 

FORMAT(' ENTER PATH LENGTH (KM)') 

READ(5f 120)PL 
IF(NRFC.EQ. 1 )G0 TO 20 . 

IF(NRFC.EQ.2)G0 TO 30 
IF(NRFC.E0.3)G0 TO 40 
WRITE(6f 160) 

FORMAT(' ENTER RAINFALL RATE (MM/HR)') 

READ(5» 120)RR 
UIRITE(6»230) 

FORMAK' SINGLE SIZE OR DISTRIBUTION? (0»1)') 
READ(5»120)NS 
IF(NS.EQ. 1 )G0 TO 820 
URITE(6f 170) 

FORMAT(' ENTER PARTICLE DIAMETER (MM)') 

READ (5t 120) DM 
DM=DM/10. 

DIAM=DM 
GO TO 820 
WRITE(6.. ISO) 

FORMAT(' ENTER VISIBILITY OF LANDFOG (KM)') 

READ(5» 120)VIS 
GO TO 820 


9 
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I 


r-* 

1 

I 
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^rs 
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53: 

30 

54: 

190 

55: 


56: 


57: 


58: 

40 

59: 

200 

60 : 


6i: 

820 

62: 


63: 


64: 


65: 


66: 


67: 


68: 


69: 


70 : 


7i: 


72: 


73: 

4 

74: 


75: 


76: 


77: 


78: 


79: 

o 

80 : 

70 

81 : 


82: 


83: 


84: 


85 : 


86: 


87: 

71 

88: 


89: 


90 : 

4 

91 : 


92: 

5 

93: 


94: 


95: 

6 

96: 

7 

97: 


98: 


99: 

220 

100 : 


101 : 


102: 

999 

103: 


104: 



URITE(6fl90) 

rORhAK' ENTER VISIBILITY OF SEAFOG (KM)'> 

READ(5»120)VIS 

IFCVIS.LT.O. 1 )GO TO 20 

GO TO 820 

WRITE<6»200> 

FORMATC' ENTER VISIBILITY OF CLOUDS <KM>') 

READ<5» 120)VIS 
CONTINUE 
DO 70 I=lfl2 
DM*DL < I ) 

CALL SIZE<RR»DH»XN» I > 

IF(NS.EQ.O)DM=DIAM -- - 

IF(NRFC.EQ. 1 )DM=1 .E-03 

IF < NRFC . EQ . 2 ) DH=2 . E-03 

IF < NRFC . Ed . 3 ) DM=1 . E-03 

X=PI*DH/U 

Y=Z1#X 

CALL ZACALC< Y»2ArNN) - . . 

CALL ZDCALC(X»ZDf ZBOfNN) 

CALL ABCALC(2ArZB»ZB0»Zl »X»ZAN»ZBNf NN) 

CALL GC ALC < X f Z AN » ZBN , QSCA r QEXT » GBK » GABS » NN » I ) 

CALL ANQCAL<X»ZAN»ZBN»QSCA»ANGL»NN» I ) 

DO 2 

CALL PTCALC<D»F>T»M»NN) 

CALL SCUMP(UrZANfZBN»P,T»ZSi»ZS2»SFpI fM»NN»SlAS»S2AS»I ) 

CONTINUE 

CONTINUE 

IF(NRFC4EQ.1)G0 TO 4 
IF (NRFC. EG. 2) GO TO 5 
IF<NRFC.EG .3)G0 TO 6 
IF(NS.EQ.l >G0 TO 71 

CALL RAIN(RR»PLfGEXT»DMfATNrATNrBf I) 

GO TO 7 ... 

CONTINUE 

CALL RAINXN(RRfPL»QEXTfDLf ATN»ATNDB»XN) 

GO TO 7 

CALL LNDFGG<DMfVIS»PL»GEXTf ATNf ATNDBr I > 

GO TO 7 

IF(VIS.LT.O. DGO TO 4 

CALL SEAFOG(DM»VISfPLfQEXTf ATN»ATNDB» I ) 

GO TO 7 

CALL CLOUDS(DhfVISfPL»OEXTfATN»ATNDB» I > 

CONTINUE 

URITE(6» 120)PLf ATN»ATNDB 
URITE(6»220) 

FORMAT(' REPEAT? ( 1=YES f 0=NO ) ' » / ) 

READ(5» 120)NREP 

IF(NREP.EG. 1 )GQ TO 800 

CONTINUE 

CALL EXIT 

END 
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1 i 

21 

31 

41 

5! 

61 

7! 

8 ! 

91 

101 

111 

121 

,131 

141 

151 

16; 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 
.32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 




1 

2 

6 


3 

4 


7 


SUHROUTINE SIZECRRf liMf XN»I ) 

DIMENSION XKIO) »X2(10)fX3<10) »X4(10) rX5<10) rX6(10> >X7(10) 
DIMENSION Xrt(lO) fXB(lO) »R<7) »CUT<7) »XN<12) 

DATA Xl/3.90906»-1 *2236E02f 1.6357E03»-5. 13497E03» 
15.4409E03»-1 .463E03»0. rO. »0. »0./ 

DATA X2/-1 .5255E01»6.38345E02»-9. 14623E03f6»9662E04» 
l-2.69783E05»5.42977E05f-5.4S427E05»2. 16777E05f0. »0*/ 

DATA X3/-9.9461 f5.81251E02»-l .2434E04» 1 .38727E05» 
l-7.68E05»2.29E06»-3.7041E06f3.27524E06»-1.161E06»0./ 

DATA X4/2.3709E01 »-9.B352E02f 1 .42472E04»-7.68455E04f 
11.91541E05 f- 2.25279E05»1.01277E05»0. fO. rO./ 

DATA X5/6.5918E01 »-2.8817E03f 4.53181E04»-3. 10239E05» 

11 . 10129E06r-2, 13855E06»2. 16497E06»-8.9636E05»0. »0./ 

DATA X6/1 .E02»-4.9157E03f8.63391E04»-6.93122E05» 
12.90943E06r-6*68496E06»7.98558E06»-3.88423E06»0. »0*/ 

DATA X7/-5.4314E01» 1 .72908E03»-1 . 1 4533E04 » 2 . 79273E04 » 
i-2.29091E04»0. »0.»0«»0.f0./ 

DATA R/150. »100. »50. »25. »12.5»5. f2.5/ 

DATA CUT/. 53 » . 58 » . 53 » . 53 f . 48 » . 48 f .38/ 

IF(RR.GE. 150. )G0 TO 1 
IF(RR.LE.2.5)G0 TO 4 

IF(RR.LT. 150. . AND.RR.GE. 100. )G0 TO 7 
IFCRR.LT. 100. . AND.RR.GE.50. )G0 TO 11 
IF(RR. LT. 50. .AND.RR.GE. 25. )G0 TO 14 
. IF<RR.LT. 25. .AND.RR.GE. 12. 5)G0 TO 17 
IF (RR.LT. 12.5. AND.RR.GE. 5. ) GO TO 20 ‘ ' 

IF<RR.LT.5. . AND.RR.GT.2.5)G0 TO 23 
DO 2 K=1.10 
XA<K )=X1 (K) 

CONTINUE 

CT=CUT(1) 

XN(I)=0. ' - ■ - 

DO 3 K=l»10 

XN( I )=XN(I )+XA<K)*DM*!j(’<K-l ) 

IF(DM.GE.CT)XN(I )=0. 

CONTINUE 
GO TO 999 

DO 5 K = 1 »10 ' ■ ■ - . - . . 

XA(K:=X7<K) 

CONTINUE 

CT=CUT(7) 

GO TO 6 
A=150. 

B=100. ....... 

DO 8 K=l»10 
XA(K)=X1 (K) 


47 

48 

49 

50 

51 

52 


8 


• XB<K.V=X2(K) 
CONTINUE 
GO TO 50 
A=100. 

B=50. 

CT=CUT<2) 
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53 1 


DO 9 K=l»10 

54: 


XA(K)=X2<K) 

55: 


XB(K)-X3(K> 

56: 

9 

CONTINUE 

57: 


GO TO 50 

58 : 

14 

A=50* 

59: 

- 

B=25. • •• “ 

60 : 


CT=CUT(3) 

61 : 


DO 15 K=l»10 

62: 


XA(K)=X3(K) 

63: 


XB<K)=X4(K) 

64: 

15 

CONTINUE 

65: 

• 

GO TO. 50 - . - - 

66: 

17 

A=25. 

67: 


B=12.5 

68: 


CT=CUT<4) 

69: 


DO 18 K=l»10 

70 : 


XA<K)=X4<K) 

7i: 


XB(K>'=X5(K) 

72: 

la 

CONTINUE 

73: 


GO TO 50 

74: 

20 

A=12«5 

75: 


3=5. 

76: 


CT=CUT(5) 

7?: 


DO 21 K=l»10 

78: 


XA(K)=X5(K) 

79: 


XB(K)=X6(K) 

80: 

21 

CONTINUE 

8i: 


GO TO 50 

82: 

23 

A=5. 

83: 


B=2.5 

84: 


CT=CUT(6> 

85: 


DO 24 K=lflO 

86: 


XA(K)=X6<K) 

87: 


XB(K)=X7<K) 

88: 

24 

CONTINUE 

89: 

50 

CONTINUE 

9o: 


XN1=0. 

9i: 


XN2=0. 

92: 


DO 51 K=lrlO 

93: 


XNl=XNHXA<K)*Dh*)»<(K-l ) 

94: 


XN2 = XN2 + X3<K)*rih*>K(K-l > 

95: 


IF<DM.GE.CT)XN1=0. 

96: 


IF<DM.GE.CT)XN2=0. 

97: 

51 

CONTINUE 

98: 


S=<XNl-XN2)/(A-FO 

99: 


XN( I )=XN2FS*(RR-B) 

loo: 

999 

CONTINUE 

loi: 


RETURN 

102 : 


END 



i 


I 


I 

i 



t 



I 


i: 

2{ 

3! 

Ai 

5 : 

6i 

7; 1 

8: 

9: 

10: 

iij 

12: 
13: 2 

14: 

15 1 
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SUBROUTINE ABCALC ( ZA . ZB » ZBO » Z1 » X » ZAN r ZBN i NN ) 

IMPLICIT COMPLEX<Y»Z) 

DIMENSION ZA(30) f ZB(30) »ZAN(30) »ZBN<30) r ZXA ( 30 ) » ZXB < 30 ) 

DO 1 N=1»NN 
ZXA<N)=ZA(N)/Z1+N/X 
ZXB < N ) =Z l:4cZ A < N ) +N/X 
CONTINUE 

ZAN<1 )=(ZXA< 1 )*REAL<ZB(1 > )-REAL(ZBO) )/(ZXA<l)*ZB< 1 )-ZBO) 

ZBN<1 )~<ZXB(1 >*REAL<ZB(1) )-REAL(ZBO> )/(ZXB<l)!l!ZB(l)-ZBO) 

DO 2 N=2»NN 

ZAN(N)=(ZXA(N>*REAL<2B<N) ) -REAL < ZB < N-1 ) ) )/(ZXA(N)*ZB(N)-ZB(N-l 

ZBN<N>=<ZXP(N)*REAL<ZB<N> )-REAL(ZB(N-l) ) )/(2XB(N)*ZB<N)-ZB(N-l 

CONTINUE 

RETURN 

END 


IJ , SUBROUTINE ZACALC ( Y » ZA t NN ) 

21 TMPI THTT rnMPI 



21 

IMPLICIT COhPLEX< Y»Z) 

1 

31 

DIMENSION 2A(30) 

i 

41 

ZA0=CC0S<Y)/CSIN< Y) 


51 

ZA( 1 )=-l/Y+l/< 1/Y-ZAO) 

61 

DO 1 N=2fNN 

71 

ZA(N)=-N/Y+i/(!J/Y-ZA(N-l> > 


81 1 

CONTINUr 

1 

91 

RETURN 

1 

i 

> 1 

10 : 

END 

11 

SUBROUTINE ZBCALC ( X » ZB t ZBO » NN ) 

1 

1 

21 

IMPLICIT COMPLEX<Y>Z) 


31 

DIMENSION ZB<30> 

f 

41 

ZBO=CMPLX(SIN(X) »COS(X) ) 

51 

ZBM=CMPLX(CQS(X) »-SIN(X) ) 

{ 

61 

ZB( 1 )=ZB0/X-2BM 


71 

ZB<2)=ZB< 1 >*3/X~ZB0 

) 

81 

DO 1 N=3»NN 

t 

91 

ZB<N)--(2*N-1 >5^2B(N-1)/X-ZB(N-2) 

i 

101 1 

CONTINUE 

i 

111 

RETURN 


121 

END 
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H SUBROUTINE SCCnP(U,ZAN.ZiN,P,T.2Sl,2S2,S,PX,n.HH,SlAS,S2AS.I ) 

2l ir.PLICIT CONPLEX(V,Z) 

3t OinCNSION ZAN(30}.ZBN(39),P<3O),T(39).ZSU36,12>,ZS2(36.12) 

41 OIRENSION S(36.t2).SlAS(36,l2),S2AS(36,i2> 

5« 2Sl<H,n-0. 

6t ZS2(ft.I)-e. 

71 00 1 N«1,NN 

• t ZSKN.X J»25l<n, 1 )*<2«N»1 )»(ZAN(N )tP(N )*ZBN(N )*T(N ) UCNt(N»l) ) 

9t zsacn,! ).ZS2tH,l )*(2tN+l )t(ZBN(N)*P(N)*2AN{N)»T(N» )/(N«(N*l ) t 

let 1 continue 

lit S<n. I )• .s*( .5*U-'PI )»«2*(CPBS(2S1(H, I ) >x>2^CAB$(ZS2<n«I ) )*ta> 

I 2 t siAS(n,i)>CABS(ZSi<n, n US2 

13t S£AS(n.I )-CAES(ZS2(n.l ) )tt2 

14f RETURN 

ISt END 






i: SUBROUTINE RAIN(RRfPLrQEXT»DM»ATNr ATNDB» I ) 

2J DIMENSION VA(6) »QEXT< 12) 

3J DATA VA/-.O2366»47,30f-81 .568f38.94»26.056»-21 .696/ 

4: V=0. 

5: DO 1 K=i»6 

6J , U=U+UA(K>*DM)«(*<K-1 > 

7J X ’ CONTINUE 

8J ATN=EXP(-RR!»cpL*QEXT( I >/(24.*V)KDM) > 

9J ATNDB=10.*AL0G10(ATN> 

lO: RETURN 

11 : END 


I ■! 


i: SUBROUTINE RAINXN ( RR » PL » QEXT » DL . ATN f ATNDB » XN ) 

2t DIMENSION DL( 12) r GEXT<12) »U< 12) r XN< 12) »ONUD< 12) 

3J DIMENSION UA(6) 

4J DATA VA/-.02366»47.38f -81 .568»38.94f 26.056»-21 .696/ 

S: SUM=C. 

61 DO 2 J=l»12 

7J V(J)=0. ■ " * - • - - -- 

8*. DO 1 K=l»6 



9*. 


U< J)=U< J)+UA<K)*DL( J)**(K-1 ) 

lot 

1 

CONTINUE 

ii: 

r> 

CONTINUE 

12 : 


DO 6 J=l»12 

13 : 


XN( J)=XN( j)/ioa. 

14 : 

6 

CONTINUE 

15? 


DO 3 J=lrl2 

16: 


QNUD < J ) =-OEXT < J ) )BXN ( J ) / < V ( J ) *DL ( J ) ) 

17 : 

3 

CONTINUE 

18: 


DO 4 J=1.12 

19: 


SUM=SUM+QNVD( J) 

20 : 

4 

CONTINUE 

21 : 


P0U=PL*RR*SUM/24. 

22 : 


ATN=EXP(-POU) 

23 : 


ATNDD=10.+AL0G10(ATN) 

24 : 

120 

FORMAT < ) 

25 : 


RETURN 

26: 


END 




r 


I 
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i: 

2: 

21 

.4: 

s: 

6 : 

7: 

8t 

9: 
lot 
u: 
12t 1 

13: 

14 : 

is: 

16 : 

17: 

18: 


SUBROUTINE QCALC < X » ZAN» ZBNt QSCAt QEXT »QBK »QftBS » NNt I > 


IhPLICIT COMPLEX<Y.Z> _ 

DIMENSION ZrtN<30)fZBN<30) 

DIMENSION OSCA< 12) »QEXT < 12) »QBK< 12) »QABS( 12) 

SUM1*0. 

SUM2=0. 

ZSUM3*0. 

DO 1 N=lrNN 

SUM1=SUM1+(2*N+1 )*(CABS<ZAN<N) )**2+CABS<ZBN(N) )**2) 


SUM2=SUM2+<2HNT1 )*REAL(ZAN<N)iZBN<N) ) 

ZSUM3=ZSUM3+(N+.5)*<-l )**N*(ZAN<N)-ZBN<N) ) 

CONTINUE 

QSCA( 1 >=2.*SUM1/X**2 

QEXTX I >=2.5*‘SUM2/X*#2 

QBK < I ) =4 . *CABS ( ZSUM3 ) **2/X#*2 

QAEtS(I)=QEXT<I)-QSCA(I) 

RETURN 

END 


1: SUBROUTINE ANGCAL(Xr ZAN»ZBN»CJSCA»ANGL»NN t ) 

2 : , IMPLICIT COMPLEX<YrZ) 

3 : DIMENSION ZAN(30> »ANGL< 12) »QSCAa2> »ZBN(30) 

4 : ANGL=0. 

S: DO 1 N=lfNN 

6: ANGL( I )=ANGL( I > +N)J{ ( N+2 ) / ( N + 1 ) *REAL ( CONJG < ZAN < N > )*ZAN(N+1 ) + 

7 : 1C0NJG<ZBN<N> )*ZBN(N-n ) )-f (2*N+ » V <N* ( N f 1 ) )4c 

.8:, 2REAL(C0NJG(ZAN<N) )*ZBN<N> ) 

' 9 : I CONTINUE 

10 : ANGLCI)=ANGL<I)*4./X**2 

11 : ANGL( I )=ANGL( I )/QSCA< I ) 

12 : RETURN 

13: END 


1: 

2: 

3: 

4: 

5: 

6 : 

7: 

. b: 

9: 

10: 

11: 

12: 

13: 

14: 1 
15: 
16: 


SUBROUTINE PTCALC ( D » P » T r M » NN > 

DIMENSION P(30) »T(30) »D(36) 

E=COS(D(M) ) 

G=SIN(D(M) )*SIN(D(M) ) 

P0=0. 

T0=0* 

P( 1 )=1 . 

T(1)=E 

P(2)=3«*E 

T<2)=3,*C0S(2.*D(M) ) 

DO 1 N==3»NN 

P<N)=P<N-1 )*E*<2*N-1 )/<N-l )-N*P<N-2)/(N-l ) 

T < N ) =E* < P < N ) -P ( N~2 ) > - < 2*N- 1 ) *G*P ( N-1 ) +T < N-2 ) 

CONTINUE 

RETURN 

END 
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SUBROUTINE LNDFOG(DMrUISrPL»OEXT» ATN»ATNDB» I ) 

DIMENSION QEXT(12) 

C®A*98E""04 4 K~tQ \ /FiM \ 

(MN=EXP(~C*PU:<caEXT<I>iVIS**<-1.538)/Dh) 

ATNDB=104*AL0G10(MTN) 

RETURN 

END 


SUBROUTINE SEAFOG( DMrUIS»PL»QEXT »ATN» ATNEiBr I ) 

dimension QEXT(12> . _ 

C=2.3^E~03 

ATN=i: XP < -C*PL*QEXT < I ) *OIS*# < -1 . 43 ) /DM ) 
ATNDD=10. *AL0G10( ATN) 

RETURN 

END 


: ' ' SUBROUTINE CLOUDSCDMt UIS»PL?QEXTi ATN» ATNDB» I ) 
: DIMENSION Q£XT(12) 

: C=4498E-04 

: ftTN=EXP<-C*PL)<cQEXT( I )*uis**(“l .538)/DM> 

: ATNDB==il0.:4:AL0G10(ATN) 

: RETURN . . 

: END 
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ABSTRACT 

Data obtained with automated anthropomt < ric data acquisition equipment is 
subject to short term errors. These errors are due to random reflections 
of light, masking of the light rays and other types of interference, optical 
and electrical. These signals are impossible to eliminate from the initial 
data produced by the television cameras. There is a need to ascertain which 
data values are erroneous and to replace them with corrected values, if poss- 
ible. This is primarily a software problem and requires a digital computer 
to refine the data off line. This paper discusses the use of tne least squares 
method for this purpose. 


Anthropometry has been defined as the practice of measuring the parts and 
proportions of the human body. In particular, NASA's interest in anthro- 
pometry arises from the need to know, quantitatively, the limits of physical 
performance imposed on an astronaut by a space suit at a spacecraft work 
station. The specific interest of this paper is the arm reach envelope of 
an individual as shown in Figure 1. From a simplistic standpoint, all that 
one would need to determine the arm reach of a subject would be a tape 
measure, a pad of paper for recording data and the subject in an environment 
of interest. As a practical matter, the time required to measure and tabulate 
this information becomes a significant value if a detailed picture is needed. 
F'jrther, the data would need to be entered into computer files for use in 
work station studies. 

There are techniques that involve high speed motion pictures which are 
projected so that the data can be extracted manually, frame by frame. This 
minimizes the time required of the subject but still requires extensive time 
for technician to collect data. 

Electronic means are now available to expedite this work. An automatic 
anthropometric data acquisition system is in use at Johnson Space Center. 

The system, developed by Southv/est Research Institute of San Antonio, TX. 
utilizes three modified television cameras in conjunction with a microprocessor 
controlled interface unit and a Digital Equipment Corporation (DEC) MINC-ll 
minicomputer. (See figure 2.) The interface unit is synchronized with the 
television cameras and determines the coordinates of a bright light source as 
it appears in the televised picture. The minicomputer stores the coordinate 
data on a floppy disk in a format that is acceptable to the mainframe computer 
in the Data Processing Laboratory. The large computer, using the PLAID program. 



generates the human subject's arm reach envelope and displays it on a 
video monitor or draws it on a sheet of paper on a plotter. (See figure 3.) 

The reach data files are also available for work station design studies. 

One of the difficulties with the television aided system is that errors 
occur in the data acquisition process. A small light bulb is attached to 
the arm of the person on whom reach data is desired. In the darkened room, 
the television cameras can easily pick up the bright point of light. The 
two principal sources of error are due to reflections of the light from 
surrounding objects and due to the subject inadvertantly blocking the camera's 
view of the light. In the latter case, the computer enters zeros into the 
data set when a particular camera fails to receive an image of the light. 

A third type of error, much less in evidence, can be most reasonably ascribed 
to electrical interference coming from sources external to the system. 

This type of disturbance, in so far as data correction is concerned, is treated 
in the same fashion as are unwanted light reflections. 

Each television camera sees only a two dimensional view of the subject. Thus 
a minimum of two cameras would be required to gain the information needed for 
determining a three dimensional reach envelope. Three cameras, viewing the back, 
side and front of the subject, provide some data redundancy so that even though 
one unit may not provide correct data, another may supply it. All three cameras 
shown in figure 2 are capable of indicating the vertical displacement of the 
light. The front and back cameras both can provide side to side motion data. 
However, only the side viewing camera can determine back and forth motion. 

Thus while some positional errors can be detected by comparing similar data 
points, this has limited potential. Other measures must be taken. 



The approach that seems to .hold the greatest promise is based on the fact that 
the subject cannot move his arms or any other part of his body at an infinite 
rate. That is, true motion is continuous in time and space without any sudden 
excursions. The type of errors that appear in the data from this system have 
this characteristic. They render values for position changes that are not 
physically possible in the short interval between sampling times. 

The system is calibrated so that the horizontal and vertical coordinates are 
set at 256 points for each axis. Thus the actual distance represented by a 
unit distance either vertically or horizontally will depend on the focal length 
of the camera lens and the subject to camera distance. In the system used at 
the Johnson Space Center Anthropometric Measurements Laboratory, the height 
and width of the field viewed is about 2 meters. Thus one unit change in a 
data value represents .78 cm in actual displacement. Since a camera is polled 
20 times per second, the time between two consecutive data values is 50 milli- 
seconds. Thus if two consecutive data values show a change of 50 units, this 
y*epresents a speed of 8 m/sec. Since most measurements are made with slow 
subject motion, a 50 unit change reflects an error and not actual displacement. 
If, to further belabor the point, the actual speed during the previous inter- 
val had been 1 m/sec. and the data indicated a speed of 8 m/sec. in the following 
interval, this sudden velocity increase would require an approximate acceleration 
of 140 m/sec. 2. The force required to move one's hand with this acceleration 
is quite beyond human capability. 

Having established the point that a subject will move in a continuous, uniform 
manner, how do we not only detect abrupt changes in data values with the computer 
but also make good approve i mat ions of what the values actually were? The approach 



used here fits an equation to a limited number of data points adjacent to 
the value In question. A predicted value Is determined and compared to the 
suspect value. If the difference Is greater than a predetermined tolerance, 
the original value Is replaced with the predicted value. The program then 
advances to the next new data point and repeats the above testing procedure. 

Using the least squares curve fitting technique, straight line and parabolic 
approximations have been used for the predictive equation. Ostensibly, the 
parabolic curve should provide a better fit than the straight line. However, 
computational Instability has been found to occur more frequently with a 
second degree polynomial equation where the -gaps due to erroneous data are 
large and the adjacent valid data points are subject to moderate, random 
error or “noise". Thus the straight line approximation based on a few data 
points is felt to be the better choice. 

An additional constraint that has been placed on this development of an al- 
gorithm for error correction is that It Is to be eventually used with a micro- 
computer of limited speed and memory. Thus the method should be straight for- 
ward and amenable to machine language coding. 

The development of the least squares error method of curve fitting will not 
be discussed here since it is found in most texts on elementary numerical 
analysis. The interested reader may find an article by F. R. Ruckdeschel 
helpful. (1) A more thorough treatment is offered by C. Daniel and F.S. Wood. (2) 


A modification to the measuring technique Is suggested. A11 three cameras 
are physically situated In a plane parallel to the floor. This means that 
the vertical displacement Is measured by all three cameras. One horizontal 
axis Is monitored by two of the cameras and the other horizontal axis Is seen 
by only one. Thus when the latter camera has Its view blocked by subject 
motion* then the information from the other two cameras Is not useful. All 
three measurements are required. 

The reconinendation Is that two of the cameras be kept In the present position 
viewing the front and side of the subject but the third camera should be 
placed so that It has an overhead view. This would mean that each of the 
three coordinates would have double redundancy. Loss of view of the light 
by the cameras is one of the major problems in using the system. Since 
celling height is not great in the laboratory, either a wide angle lens 
could be fitted to the overhead camera or else a large mirror could be pos- 
itioned overhead and the third camera could be placed in a more convenient 
position. 

An alternative scheme would be to use a fourth camera that could be under 
the control of the microprocessor and its output could be used in place of 
that of the camera that has Its view blocked. 
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Abstract : 

Banvis A, 5 and 7 false color composite photographs were prepared using 
cuiltlcmporal data from LANDSAT scenes acquired during April 1977 and 
April 1981 on coraputer compatible tapes, and these color composites 
vorc coraparod with band 7 black and white photographs prepared for the 
ontlro county. Delineations of soil boundaries at the soil association 
level were achieved using LANDSAT spectral reflectance data and slope 
naps for a portion of Carroll County, Missouri. Forty two spectral 
reflectance classes from April 1977 LANDSAT data were overlaid on 
digitized slope maps of nine USCS 7.5 minute series topographic quad- 
tangle slope maps to achieve boundary delineations of the soil associations. 
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A nocd exists in the state of Missouri to use LANDSAT data 
to aid In Soil Survey mapping. The Missouri Depai. ' ■ r of Natural 
Resources has responsibility for conservation programs In the state 
and is committed to accelerating the National Cooperative Soil Survey 
Program. Research to date has shown strong correlation between the 
spectral reflectance characteristics of soils and some of the soil 
prv'perties such as texture, organic matter content, drainage, color, 
etc.' Millc LANDSAT classification of soils will not alone distinguish 
between different soils which exhibit similar spectral signatures, 
it Is possible to prepare a product which will aid the soil m<:pplng 
process by delineating patterns, and possible divisions of soil asso- 
ciations. 

A party leader has been assigned by the state soil scientist of 
Missouri to begin napping this summer a progressive soil survey 
of Cari'v'll County, Missouri and Saline county, Missouri next year. 

The stale soil scientist has indicated that it would be of consld^’rable 
value to tho party leader if he could be provided maps based on LANDSAT 
analysis to help In the progressive soil surveys of Carroll and Saline 
counties, Missouri. Approximately two thirds of Missouri still needs 
a detailed county soil survey to assist land managers in their dec- 
isions, A roconi letter from tho state conservationist of Indiana 
t<* the state conservationist of Missouri Indicated that the detailed 
Co^mly Sell Survey of Jasper County, Indiana was completed with a 
savings v'f IS nan-years lime. Since the Missouri Department of Natural 
Reseurees would like to complete in the 1980-1990 time frame modern 
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detailed covmty soil survey reports for the entire state it appears 
that LANDSAT analysis to aid in soil survey mapping would help them 
to acheive this worthwhile goal. 

Stu dy Area 

Carroll County, Missouri was selected as the study site because 
(i) preparations are underway to procure 7.5 minute USGS orthophoto- 
quads for use as a base for conducting a detailed soil survey of this 
county by the party leader, (ii) this county contains more Missouri 
River bottomland than any county in the state (A), (iii) and the large 
tracts of farmland In the bottomlands are more likely to show relation- 
ships with the LANDSAT spectral reflectance data chan land subdivided 
into smaller units, 

Carroll County, Missouri is situated in the northwestern part of 
the state, 66 nilcs oast of Kansas City and 211 miles northwest of 
St, Louis. It is bounded on the north by Livingston County, on the 
vest by Caldwell and Ray Counties, and is separated from Saline and 
Lafayette Counties on the south by the Missouri River and |?rora Charjton 
Cc'unty on the c.ist by the Grand River. It has an area of 703 square 
miles with A49.920 acres in the county. The last Carroll County soil 
sui vey published was in 1912 (A), A soil survey of Livingston County 
was published in 1956 (3), one for Caldwell County was published In 
1974.0)', one for Lafayette County was published In 1975 (2), and a 
Si'll inventory of Chariton County from digital analysis of satellite 
and sc.»nner data was published In the Soil Science Society of America 
Jv^uinal In 1977 (5), These publications were useful as indicators of 
tho soil series to be expected In Carroll County. The soils of the 


I 
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area formed In deep loess and glacial till In the uplands and in 
alluvial sediments In the floodplains. 

Data Set Formulation 


UNDSAT-2 MSS data collected on 26 April 1977 were utilized as 
the main data source for this study and LANDSAT-2 MSS data collected 
23 April 1981 were utilized to prepare band 7 Polaroid photographs 
of Carroll County. These scenes were selected because the data were 
(1) of high quality, (il) acquired when most row crop cropland is in 
a bare i;oil state, and (ill) free of interfering atmospheric and sur- 
face conditions (l.e., clouds, haze, saturated soils, and standing 
water.) 


The LANDSAT MSS data were geometrically corrected (l.e., rotated, 
deskewed, and rescaled to an approximate scale of 1:24,000.) That por- 
tion of the scone containing the eastern part of Carroll County (Big 
Crock and Hurricane Creek Watersheds) were selected for further pre- 
prv'<Ci. ssing. These data were registered to ground control points 
selected from USGS 7.5 minute topographic quadrangles. This procedure 
produced a data set of an exact scale of 1:24,000 and registered geo- 
gr.iphlc points in the data to their exact ground position. Using these 
same quadrangles slope maps produced by USGS were digitized and sub- 
sequently overlaid and registered to the LANDSAT data. Six physio- 
graphic position boundaries (water. Bottomlands, 0=5Z slopes, 5-9% 
slopes, P-14% slopes, and greater than 14% slopes) were digitized and 
overlaid on the LANDSAT-2 data for the Aprlll 1977 scene. Forty two 
spectral reflectance values were obtained as an overlay and a computer 
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printout obtained for each of the six physiographic boundaries. 
Methods 

Topographic data (e.g., physiographic position) were used to 
deltiuvn. bouidarl.s at the soU association level, and was helpful 
in separating sails which were spectrally slnllar but genetically 
different. Alphanunerlc 1:24,000 spectral maps representing the 
42 classes Inventoried were produced for the following 9 quadrangles: 
Coloma, Tina, Hale, Bogard, Standlsh, Bosvorth, Brunswick West, 

Carr, llton East, Miami Station. Comparisons were made between the 
aps |r..luced b^ this technique and the soil association boundaries 
observed for the adjacent counties of Laf.iyette and Caldwell. 

Col.ar i::frared phottagraphs were obtained from the photographic 
services laboratory of J.SC/NASA at Houston, Texas for mission 2l3 
which Included the portion of the Mlsscurl River basin that forms 
the southern boundary of Carroll County. Attempts were made to celate 
colors and tonal differences in the photogiaphs to the published soil 
survey of I..ifavette County which contaiits .an aerial photographic base 
and cowpare these tones with similar features across the River 'In- 
Carroll County. A portion of the hott.amland was enlarged to an 
approximate scale of 1:24.000 using I.ANDSAT-2 data and the General 
Electric Image 100 system in the Scene Analysis Branch of the Earth 
Res..urces Research Plvision at JSC/NASA. This enlargement was sub- 
sequently photographed With the Color Polaroid Camera device so that 
Individual pixels were readily visible and could be compared to tonal 
differences seen in the Color Infrared Photographs of Mission 2l3 
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vhlch was filmed 11 September 1972. 
Results and Discussion 


The 42 spectral reflectance classes contained within the Inventory 
of Carroll County are shown in Table 1 for bands 5 and 7 with the 
character that was printed on a line printer for the 6 physiographic 
boundaries. Table 2 :-howa the soli association that were predicted 
using these data. Table 3 provides Information about the taxonomic 

characteristics of these soil series, and the parent materials of the 
soils. 


The Haynic-Leta- Waldron soil association Is found only In '^he 
Missouri River bottomlands. The Kcnneboc-Nodaway-Colo-Zook soils are 
found In the bottomlands along upland streams along with the associated 
Blackoar. Bremer, Booker, Moniteau, and Otter soli series. It is a 
relatively easy task to define the soil boundaries for these two 
associations using the computer line printer products, but it is more 
difficult to attempt separations of the upland soil associations. 

The GnmJy-Lagonda, the Marshall-Higglnsville-Macksburg, and 
the Knox-Marshall-Kennebec are upland soil associations. These asso- 
ciations could not be delineated using the computer printouts because 

llurc is an overlap of slopes on the landscapes where these soils are 
fv'und. 

Plgital mapping of soils using LANDSAT data is based primarily upon 
the surface loflectance properties of the various soils. Such distin- 
guishing soil characteristics as parent material, profile development, 
and landscape position are not observable with current satellite systems. 




Often widely different soils exhibit similar spectral responses and 
cannot be differentiated from one another using satellite data alone. 
However, appropriate ancillary data used In conjunction with the LANDSAT 
data can greatly increase the Infornationnl content of a spectral soils 



TABLE 1 


6 


MEAN VALUE OF SPECTRAL REFLECTANCE ASSOCIATED WITH EACH CHARACTER FOR 
BANDS 5 and 7 FOR THE 42 CLASSES 


CHARACTER 

BAND S 

BAND 7 

X 

29 

29 

B 

24 

30 

C 

23 

44 

D 

25 

24 

E 

29 

46 

F 

33 

32 

G 

25 

48 

H 

23 

21 

I 

26 

42 

J 

19 

50 

K 

28 

^ 5 

L 

20 

57 

N 

41 

42 

N 

36 

38 

0. 

29 

35 

P 

22 

26 

Q 

45 

44 

R 

34 

46 

S 

26 

15 

T 

17 

13 

U 

25 

37 

V 

22 

50 

W 

21 

18 

X 

13 

2 

Y 

32 

39 

Z 

32 

27 

1 

50 

47 

2 

56 

52 

3 

17 

16 

4 

19 

3 

5 

17 

67 

6 

18 

40 

7 

20 

22 

8 

19 

15 

9 

17 

55 


18 

18 

% . 

* 

25 

11 

* 

24 

5 

i 

24 

56 

$ 

•»*> 

38 

4 

40 

48 


34 

22 
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rREDlCTED"SOIL ASSOCIATIONS IN THE QUADRANGLES OF CARROLL COUNTY 


QIL\PRANGLE 


SOIL ASSOCIATIONS 


nOGARD 


nOSWORTH 


URUNSWICK WEST 


CARROLLTON EAST 


COLO.'LA 


HALE 


Ml.LMI STATION 


SIANPISH 


TINA 


GRUNDY-LAGONDA* 

MARSIIALL-HIGGINSVILLE-MACKSBURG** 

KENNEBEC-NODAWAY-COLO-ZOOK*** 

KNOX-M,\RSHALL-KENNEBEC**** ***** 

GRUNDY-LAGONDA 

HARSHALL-HIGGIHSVILLE-MACKSBURG 

KENNEBEC-NODA17AY-COLO-ZOOK 

HAYN I E-LETA-WALDRON** ** * 
KNOX-IL\RSHALL-KENNEBEC 

IL\YN I E-LETA-WALDRON 

KNOX-MARSHALL-KENNEBEC 

Kt-'NNEBEC-NODAWAY-COLO-ZOOK 

GRUNDY-LAGONDA 

MARSIIALL-HIGGINSVILLE-MACKSBURG 

KENNEBEC-NODAWAY-COLO-ZOOK 

I 

GRUNDY-LAGONDA 

MARSHALL-HIGGINSVILLE-MACKSBURG 

KENNEBEC-NODAWAY-COLO-ZOOK 

HAYN lE-LETA-WALDRON 

KNOX-RARSHALL-KENNEBEC 

KENNEBEC-NODAWAY-COLO-ZOOK 

knox->larsiw;ll-K£nnebec 

GRUNDY-LAGONDA 

MXRSHALL-HIGGINSVILLE-MACKSBURG 

KENNEBEC-NODAWAY-COLO-ZOOK 

GRUNDY-I-/\GONDA 

JLARSILALL-HIGGINSVILLE-MACKSBURG 

KENNEBEC-NODAWAY-COLO-ZOOK 


* Clhor levies: ADAIR, AILMSTER, LAMONI, LINEVILLE, SEATIOUR 

** v^thor possible series; LESLIE, MCGIRK, MINDEN 

*** Other possible series: BLACKOAR, BREMER, BOOKER, MONITEAU, OTTER 

**** Other possible series: HIGGINSVILLE, MINDEN 

***** Other possible scries: DOCKERY, HODGE, MODALE, MYRICK, RAY, 

SARPY, WAUBONSIE 
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r AKKNT MATERIAL AND TA^XONOMIC GROUPS OF THE MAJOR SOIL SERIES 
FOR THE PREDICTKD SOIL ASSOCIATIONS IN CARROLL COUNTY 


SOIL SERIFS 

SOIL ORDER 

' " ■ ■ "'1 ■ 

SUBGROUP 

PARENT MATERIALS 

COLO 

MOLL I SOLS 

CUhWLlC ILNPLAQUOLLS 


ALLUVIUM 

GRUNDY 

MOLL I SOLS 

AQUIC ARGIl'DOLLS 


LOESS 

HAYNIE 

ENT I SOLS 

TYPIC UniFLUA'ENTS 


VERY RECENT ALLUVl 

HIGG1KSV15.LE 

MOLL I SOLS 

AQUIC ARGIl'DOLLS 


LOESS 

KENNEBEC 

MOLL I SOLS 

CmaJLIC HAPLUDOLLS 


ALLUVIUM 

KNOX 

ALFISOLS 

MOLLIC HAPLUDALFS 


LOESS 

LAGONDA 

MOLL 1 SOLS 

AQUIC ARGIL'DOLLS 


LOESS, GLACIAL DK1 

kFTA 

MOLLISOLS 

AQUIC fluvFntic HAPLUDOLLS 

ALLUVIUM 

MACKSBURG 

MOLL 1 SOLS 

AQUIC ARGIUDOLLS 


LOESS 

MARSHALL 

MOLLISOLS 

TYPIC HAPLUDOLLS 


LOESS 

NODAVAY 

KNTISOLS 

TYPIC UDIFLUVENTS 


ALLUVIUM 

WALDRON 

MOLLISOLS 

AQUIC FLUX’ENTIC HAPLUDOLLS 

RECENT ALLUVIUM 


ZOOK 


MOI.MSOLS 


•CllhR’UC IL\rL.\QUOLLS 
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Introduction ! 

The purpose of this investigation is to verify certain 
performance criteria for a solid state power combining antenna 
proposed for use on the solar power satellite. 

Background ! 

The reference system solar power satellite concept (1) 
was developed around the use of the Klystron as the device utilized 
to convert DC power produced by sblar cells into microwave energy 
for transmission to ground based receiving antennas* The improve- 
ment in performance and reliability of modern microwave transistors 
offers the potential for the development of a solid state micro- 
wave power transmission system. The obvious benefits of longer 
life, smaller size and ease of manufacture provides adequate 
incentive to pursue the development of a solid state system. 

Under NASA sponsorship several studies have been performed to 
explore the potential for solid state microwave power transmission 
(2), (3).. (4), (5)» While the power handling capabilities of 
microv/avG solid state devices have improved significantly, it 
is generally agreed that some form of power combining will be 
required to attain practical levels of power in a solid state 
system. 

The realization of efficient transmission of microwave 
power dictates that a low loss power combining technique be 
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utilized. A prototype system has been designed, fabricated and 
delivered to J3C under contract by the Boeinti Corporation. The 
system incorporates a four feed micros trip antenna, a stripline 
antenna phasing network and provides for the integration of four 
transistor amplifiers. The antenna and feed network have been 
designed to operate in S^Band at 2,45 GHZ 

The power combining technique utilized is considered to be 
somewhat unique in that the power combining and radiating functions 
occur in a single cavity (2), Each amplifier is directly coupled 
to the cavity via four microstrip feedlines which are mounted on 
a ceramic type dielectric substrate. The substrate is backed by 
a light weight aluminum cavity which sums the power of the four 
sources. On the opposite side of the substrate is mounted two 
radiating slots which are each fed on the ends by two narrow slot- 
lines, The antenna behaves like two one-half wavelength slots 
coupled via a common cavity, thus the radiated energy is suraned 
to yield a single radiated phise. 

The original concept is illustrated in figure 1. The foil 
pattern for the unit delivered to JSC contains certain design improve 
ments and io shown in figure 2. Proper feed of the antenna requires 
that two of the RF inputs be 180^ out. of phase with t..e remaining 
two. Therefore, an important component of the antenna system is 
a balanced feed network that provides four 0° - 180° equal amplitude 
outputs. The stripline feed network consist of two 0° - 180°' rat- 
race type hybrid circuits connected to the output of a two-way in 
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(a) BREAK-A-WAY VIEW 



FIGURE 1 POWER C OMBINING MICROSTRIP SLOTLIN E ANTENNA 
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CAVITY "E" FIELD PROFILE 
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FIGURE 2 COPPER METALIZATION PATTERN FOR FOUR FEED MICROSTRIP ANTENNA 
(Two patterns superimposed) 
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phase power divider all integrated into a single stripline circuit. 

The remaininj; elemente of the test conr:!"-*: of four 

Silicon bipolar type transistor amplifiers each capable of producing 

•« 

1/8 watt for a combined total of 1/2 watt. ‘ 

The silicon devices were supplied as a means of testing the 
power combining technique. Other studies are concentrating on the 
development of satisfactory amplifyin;; techniques for SPS usage 

(l)i (3)» (^)f (5)* The interconnection of the antenna, amplifiers 
and feed network is illustrated in figure 3« 

The solid state antenna module has been evaluated for gain, 
pattern and efficiency on the antenna range with and without 
amplifiers by the builder. However, it was desired to obtain addi- 
tional information concerning the performance of the power combining 
prototype. The following test were accomplished! 

(1) Frequency response of the antenna and feed network 
with and without amilifiers over a band of frequencies 
ranging from 2,2 GHZ - 2.7 GHZ. 

(2) SV/R was determined for several frequencies. 

(3) Output power of the antenna was measured with one 
feedline disconnected, 

(4) An investigation of the noise "threshold'' of the 
antenna and amplifier system was attempted, 

(5) Phase tracking and phase jitter of the antenna feed 
network and amplifier system was measured. 
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INPUT 9 2.45GHz 


FIGURE 3 POWER COMBINING ANTENNA. FEED NETWORK AND 
POWER AMPLIFIER BLOCK DIAGRAM 
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Procedure and Results i 

The majority of the testin/?; was ac Ir.hod hi. the sm-'ll 

anechoic chamber located in the JSC antenna test facility. The 
phase tracking and jitter test were accomplished utilizing the 
LINCOM statistical loop analyzer system located in the ESTL facility 
at JSC. 

The test configuration for all of the measurements other than 
the SVJR and phase jitter test is shown in figure 4. The separation 
of the two antennas is approximately three meters which should 
provide adequate near field correction(2). The input power wae 
maintained at a constant level, while the signal generator was 
varied over the desired frequency range. 

The output signal was measured with the microwave power meter 
and recorded. Since the gain of the receiving antenna varies with 
frequency it was necessary to correct the raw data for the effects 
of the receiving antenna before plotting the results. Table I 
and figure 5 depicts the data for the antenna and feed network 
alone while Table II and figure 6 chows the results with the ampli- 
fiers installed. 

In reviewing these results, it would appear that the cavity 
associated with the power combining and radiating circuits is not 
a particularly good filter. 

The half power points (-3 db) occur approximately at 2.3 GHZ 
and 2.6 GHZ, These points are somewhat difficult to determine 
because of the erratic nature of this data. This would indicate 
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Standard Gain 
Antenna 


Figure Test Configuration For 
Frequency Response Measifrements 
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TABLE I 

(Antenna and Feea Network Only) 


Frequc'iioy PO 

dbin 

Receiving 

Antenna 

Gain 

Gain 

Normalized 
to 2.45 GHZ 

SP 3 

Antenna 

Re G pons e (db.m) 

22 P S 

- 7.2 

16.2 

f- 0.7 

-6.50 

22 SO 

- 5.1 

16.35 

-t- 0.55 

- 4.55 

227 S 

- 5.8 

16.44 

+ 0.46 

’ 5.34 

2300 

•■* 5.8 

16.50 

+ 0.40 

- 5.40 

232 '’ 

- 3.5 

16.55 

+0.35 

-3.15 

23' 0 


16.65 

+0.25 

-4.35 

23 ’.' 

- 3.0 

16.70 

+0.20 

- 2.30 

4 . ' ' < 

- 2.95 

16.30 

+0.10 

-2.35 

o ' » ( 


160 32 

-t- 0.08 

-2.62 

O '• r N 

4 ’r ' 

- 2.5 

16.90 

-0 

-2.50 

Zm-V ^ 

- 3.1 

16.95 

-0.05 

-3.15 

2 :' 0 v’ 

- 2.4 

17.05 

-0.15 

-2.55 

or; . 

- 4.0 

17.10 

- 0.20 

- 4.2 

r! f ' 

- 3.7 

17.15 

-0.25 

- 3.95 

2 : ' 

* 5.3 

17.20 

-n.30 

-5.6 

2(5 ' ■' 

- 6.5 

17.25 

-0.35 

-6,35 

** t 

-5.6 

17.30 

- 0.40 

-6.00 

a'- v' 

-7.0 

17.35 

- 0 . 45 

-7.45 

V 

-6.7 

17.40 

-0.5 

-7.20 
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, ^ TABLE II 

(Antnnna and Food Network with Amplifiers connected) 


Frequency 

PO 

dbm 

Receiving 

Antenna 

Gain 

Gain 

Normalized 
to 2.45 GHZ 

bPo 

Antenna 
ResDonse ^dbrnl 

2200 

<-30 

16.2 

+0,70 


2225 

-21.3 

16.25 

-♦• 0.65 

-20,65 

2250 

-18 

16.35 

+0.55 

-17.45 

2275 

-13.8 

16.44 

■♦' 0,46 

-13.34 

2300 

-13.3 

16.50 

-fO .40 

-12.90 

2325 

- 11.7 

16.55 

+ 0.35 

-11.35 

2350 

-13.8 

16.65 

•► 0,25 

-13.55 

2375 

—12 . 4 

16.70 

•♦•0.20 

-12.20 

2400 

- 7.5 

16.80 

+0.10 

- 7.40 

2425 

- 4.0 

16.82 

+ 0,08 

- 3.92 

24 

- 1.5 

16,90 

lo 

- 1.5 

2475 

- 3.6 

16.95 

-0,05 

- 3.65 

2500 

- 7.2 

17.05 

.-0.15 

- 7 o 35 

2525 

-11.5 

17.10 

-0.20 

-11.70 

25J0 

-13.0 

17.15 

-0,25 

- 12.75 

25 "’ 5 

-17.2 

17.20 

-0.30 

-16.90 

2600 

-18.3 

17.25 

-0.35 

-17.95 

2625 

-19.7 

17.30 

-o, 4 o 

-19.30 

2650 

- 26.3 

17.35 

-0.45 

-25.85 

2675 < 

-30 

17.40 

-0 , 5 



11 




OBlGU'irX 

OF POOR 

a bandwidth of roughly 275-350 yieldinfj a low Q for the cavity. 
Since the IMS band where the 3P3 will likely operate *1 j orily 100 M}IZ 
wide* it would appear that some form of filter would be required 
if an antenna similar to this one is used. The results for the 
antenna with the transistors installed was improved considerably. 

The approximate half power points occur at 2.420 GHZ and 2.485 GHZ 
yielding a bandv/idth of 65 I4HZ. The improved performance with the 
amplifiers is attributed to the fact that each of the four amplifiers 
were gain and phase matched to operate at 2.45 GHZ by the builder, 
thus some filtering is supplied by amplifier action. 

Additional questions concerning the filtering capabilities 
of the power combining cavity v/ere raised by swept fx-equency 
measurements performed on the solid state module utilizing the 
Hewlett Packard network analyzer. The results were obtained in 
the form of Smith chart plots. These plots indicate a number of 
frequencies outside of the IMS band where the VSIVK at the input 
to the module is 1.5 or less indicating that energy is being ab- 
sorbed into the antennE. system and would possibly be transmitted 
by the antenna. Figures 7 and 8 are examples of the plots obtained, 
v/hile Table III is a list of some of the out of band low V3V/R 
points. 

During the testing of the antenna and feed networki prior to 
the insertion of the amplifiers, one leg of the feed network was 
disconnected from the antenna and the effected parts left in an 
open circuit condition to simulate the effect of the failure of 
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TABLE III 


Frequencies with VSWR less than I.5 
2.^5 GHZ 
3.32 - 3.37 GHZ 
3.75 GHZ 
3.88 - 3.98 mz 
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Pig»jre 7 VSWR Measurements between 2,2 and 2,7 GHZ 
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one transistor. The output power was measured at the center 

frequency of 2.45 The obnorvod pov/or v/.'ui ':>[•! *ib 

compared to the results with all fcedlinos connected. The expected 
result when there is a 25J^ loss of power is about -1.25 db. This 
observation was discussed with the builder of the antenna (6) and 
it was concluded that the reflection coefficient introduced into 
the cavity from the open fecdline probably would account for the 
additional reduction in output power. 

The final tost conducted in the anechoic chamber involved 
an attempt to check the “noise threshold” of the power combining 
module. In this test the SPS antenna with amplifiers was placed 
directly in front of the receiving antenna while the signal source 
was attenuated by a 0 - 99 db step attenuator. On the detector 
side I the HP 432A pov/er rooter v/as replaced first with an HP 435A 
power meter equipped with the HP 8434A diode type sensor and later 
with a Stoddard sensitive RFI communications receiver. The 
minimum sensitivity of the HP 8484 a sensor is -65 dbm while the 
Stoddard rccrivor and its connecting cable was calibrated to a 
noise level of -83 dbm. Neither instrument was able to detect 
a measurable level of thermal noise. It has been determined that 
the addition of a high gain amplifier of known noise performance 
was needed to accomplish this testing function. 

The phase tracking and jitter test on the antenna feed net- 
work and the silicon transistors was conducted utilizing the 
LINCOM statistical loop analyzer (SLA). 
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The SLA provides for automatic acquistion and tracking of 
two signals. In Uio tracking mode it can bo prograr.uned to acquire 
up to 1000 samples over a period of one minute. The signals are 
analyzed for a variety of information including phase difference 
and phase jitter. The test configuration utilized is shown in 
Figure 9* 

Since the SLA requires an 8,5 I®IZ input signal on each of its 
inputs I it was necessary to mix the SPS output of 2.45 GHZ with 
the output of a 2.4415 GHZ reference oscillator to obtain the 
desired 8,5 MHZ SLA input signals. The mixers are followed with 
8,5 MHZ bandpass filters. As noted on the diagram, the unused 
amplifiers were connected to 50 ohm terminations. 

The signal sources are synthesized type signal generators 
with stable phase tracking and jitter characteristics. This is 
indicated by a measured mean jitter of 0,827° for the two oscill- 
ator mixer combinations. 

The SPS antenna feed network and amplifier system was connected 
to the SLA and the two 180° ports (1 - 4) were selected for com- 
parison, Later the cables were rotated and the 0° ports (2-3) 
were tested. Following this, port one was selected as a reference 
and was sequentially compared with each of the remaining ports. 

The data for these measurements are plotted in figure 10, The 
data that are plotted represents a meam of all data runs. 

The RMS phase jitter plot has been corrected by subtracting 
the RT-IS value of the oscillator-mixer phase jitter. As noted by 
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Transistor 



PiguTG 9 Test Configuration 
For Phase Jitter Measurement 
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MKAN VALUE OP AI/IPLIl lER JITTER 


Amplifier Pair 

Mean Jitter Amplifier 
with Oscillator Mixer 
Noise Included _ 

RMS Value of 
Amplifier Jitter 

1 & 4 

1 . 72 ° 

1 . 51 ® 

2 & 3 

1 . 17 ° 

0 . 82 ° 

1 & 2 

H 

• 

00 

o 

1 . 68 ° 

1 & 3 

1 . 07 ° 

0 . 68 “ 

Oscillator-Mixer 

0 . 827 ° 



•rms with 
mixer noise 



o 


_ I --J U- 

1&2 1&3 

Amplifiers 


Q 


J„ 

2&3 


rrms 

amplifiers 

only 



Figure 10 Amplifier Phase Jitter 
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SUA MODEL. ±Oi-e 


P.HftSIL.LB R QR->i!IT.H R S LIP MJi 


— TEST CONDITIONS -- 

THE LOOP TYPE UNDER TEST IS CM CARRIER 
THE TRACKING BANDMIDTH • A 200.00 HZ 
TEST PARAMETER - 1 TO 4 H 

SIGNAL POWER TO NOISE DENSITY RATIO - 0.00 OD-HZ 

DOPPLER OFFSET - 1.0 KHZ 

SYMBOL RATE - 1.0 KSPS 

ACCURACY OF TEST RESULTS - 95.52 S 


-- TEST RESULTS -- 

DIFFERENTIAL PHASE MEASUREMENT • 57.75 DEG 

STANDARD DEVIATION OF THE PHASE ERROR JITTER • 1.14 DEG 

AVERAGE CLOCK-WISE SLIP RATE - 0.00000000 SLIPS PER SEC 

AVERAGE COUNTER CLOCK-WISE SLIP RATE - 0.00000000 SLIPS PER SEC 

AL AVERAGE SLIP RATE - 0.00000000 SLIPS PER SEC 


Figure 11 Sample printout from the Statistical Loop Analyzer (SLA] 
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the graphs t the amplifier feed network appears to stable in terms 
of ph'i:50 jittor. Tho mean variation notod for tho 180® por^.-j (1 - 4) 
was 1«51°* These data were taken from 12 data runs* which includes 
12 1 000 samples* 

The mean variation noted for the 0® ports (2-3) involving 
2000 samples was 0.82®. Figure 11 indicates a sample of the SLA 
data output format. 

Conclusions and Recommendations i 

This exercise involved a series of test on a prototype power 
combining antenna proposed for use on the. solar power satellite. 

The primary purpose was to investigate the filtering capabilities 
of the antenna along with the phase jitter of the antenna feed 
network. The effect of transistor failure was briefly examined. 

In addition, an attempt was made to determine the noise threshold 
of the system. 

A summary of the findings followi 

1. The frequency response test indicate that some 
filtering will be needed to achieve the band- 
width necessary for operation in the IMS band. 

2. The limited swept frequency testing that was 
accomplished indicated several frequencies 
outside of the proposed band of operation 
where the VSWR of the antenna feed system was 
less than 1.5« This testing included only a 
few data runs and did not include any harmonic 
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frequencies* The implication is that energy 
wouli be abnorboci into the ant<'>r.na at 

these frequencies and probably transmitted* 

With one feedline disconnected the observed 
power drop was -2*3 db. The normal loss 
expected when the input is reduced by 255 ^ is 
roughly -1.25 db* This would indicate that 
the one open feedline causes enough mismatch 
to degrade the performance of the power 
combining cavity beyond a 25^ loss. The 
assumption is that the normal failure mode 
for a transistor is an open circuit and that 
its geometry with respect to the input port 
introduces the necessary reflection coefficients 
to further degrade the signal. 

The phase tracking and jitter measurements 
indicate that the feed network is quite stable 
with a low level of phase jitter. The use of the 
statistical loop analyzer greatly enhanced the 
acquistion of this data. No effort was made 
to correlate the absolute phase difference 
between the ports because of the difference in 
the lengths of the interconnecting transmission 
lines , 

The attempt to measure the noise threshold of 
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the SPS modulo was inconclusive due to a failure 

to flovise a cuit.itlc teatin^ proc^^durn. Thn no 

signal noise power level is apparently less than 
the -83 dbm sensitivity of the Stoddard RPI 
receiver available in the antonna testing facility* 
Recommendations for furth 1 • development include 1 

i* Gallium arsenide MSSFETS suitable for use in 
the SPS application need to be designed and 
developed. Most researchers propose a 
switching mode of operation (class DorE) (5)« 
Hov/ever, the current generation of MSSPETS are 
not well suited for operation in this mode (4). 

2. “Off the shelf gallium arsenide devices 
capable of operation at 2 - 5 watts have 
been recently adapted to operation at the 
proposed SPS frequency (4), It is recommended 
that these devices be incorporated into the 
prototype antenna system. This would provide 
for a realistic module pov/er of 8 - 20 watts. 

VJhile these devices operate in class AB and 
are somewhat short of the desired efficiency, 
they represent a reasonable increment of 
improvement. 

It is recommended that the full power prototype 
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be tested for the followinf^t 

a. Gain* directivity and 'pattern at 

full power. 

I 

b* Harmonics t out of band freq.uencies » 
and noise levels to determine the 
level of filtering required. 

Improvements be made to the microstrip antenna and 
cavity as noted by its builder (2), Improvements 
in the Q of the cavity system could greatly reduce 
the need for filtering 

Since class E and other switching modes of 
operation require the use of matching networks 
to optimize;? the performance of the transistor (?) 
it is reconunended that suitable designs be 
developed and built. Further, it is recommended 
that these be implemented in microstrip with an 
eye toward the ability to mass produce. Ideally 
any filtering requirements set forth in step 3 
could be implemented in the matching network. 
Additional analysis is needed concerning the 


effect of amplifier failure on the power output. 

The geometry and design of filtering, phase 
matching and impedance matching networks possibly 
could be used to minimize the effects of mismatches 
caused by the failure of an amplifier. 
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NONLINEAR COUPLING BETWEEN THE IN-PLANE AND 
OUT-OF-PLANE MOTION ABOUT A STATIONARY POINT IN 
GEOSYNCHRONOUS ORBIT 
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The equations of motion of a satellite In a coordinate system rotating with the 
Earth were recently given in reference 1. In this system when the oblateness 
cocrfloiont and cllipticity coefficient (J 22 ) the geopotenlial are 
given, It is a relatively simple matter to obtain the equilibrium points of the 
problem. This is done in a manner quite analagous to that of finding the 
equilibrium points of the restricted three body problem (ref.2). There are 
four equilibrium points which are found to occur in the equatorial plane. Two 
of the points are located over the ends of the minor axis of the elliptical 
figure of the Earth's equator and two of the points are located over the major 
axis. The r^dii to there points are very nearly the same or the radius of a 
satellite orbit in a 24-hour period about the Earth. 

The stability of our orbit In the neighborhood of these four points Is found by 
linearizing the equations of motion about the points. This linearization by 
definition results in solutions which are valid in our infinitesimal region 
about an equilibrium point. This type of stability analysis follows that done 
In reference 2 for the restricted three body problem. 

In this note the range of validity of the solutions in the vicinity of the 
llbratlon points will be extended to included significant out-of-plane 
oscillations by the satellite. This will permit study of the stability of 
inclined geosynchronous satellites. 
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2.0 THE DIFFERENTIAL EQUATIONS 

Following Band (1), the differential equations of motion of a satellite in a 
gravitational field with Jj (oblatencaa) and Jpg (ellipticity) perturbation forms 
Included in the potential function in a coordinate oystem rot.ating with the 
Earth are 


X - 2wy s ♦x 


y ♦ 2wx = ♦y 


( 2 . 1 ) 


V * , 

t 

where the force function, ♦ , is defined as 

» t » 

♦ = 1/2 w2 (x2 + y2) - u , 
and 

^2 . 

U - - _ 4. — (3 ain2 A - 1) 

r r3 

^22 

. cosg 4» cos*^ (X-A 22 ) » 

r3 

is the geopotential function, 


( 2 . 2 ) 


(2.3) 


where 


M = the gravitatiotial constant of the Earth, 

62 = 1/2 M J2(«.)2, 

622 - 3 m Jzz (R.)* > 0 (aln u J22 < 0 ), 

Jg is the oblateness coefficient, 

Jp2 is the ellipticity (of the Earth at the equator) coefficient, 
Rg la the Barth's radius, 


2 
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r * (x2 ♦ y2 ♦ 22)1/2^ (the radius to the satellite), 

^ is the geographic latitude of the satellite measured from the 
equator of the Earth, 

f 

X 10 the geographlo longitude of the satellite measure froai the 
Greenwich meridian, and 

X22 is the longitude of the major axis equatorial ellipse of the 
Earth (a constant) 


The quantities ^ and x ®re related to the rotating coordinates by 


.2 (x 2 ♦ y 2)^^2 

Sin AS-, cos M z — — — — — 
r ^ r 


and 


sin X 


y 

(x2 + y2)1/2 


cos X = 


x 

(x2 ♦ yi)l7i 


( 2 . 5 ) 


The values of 
experimentally, 
is of order 10“ 


«^2» ^22 > and X22 ane numbers that have been determined 

If turns out that Jg is of order 10“3, Jg2 Is of order 10"^ 
e 22 is of order 10“5 (by non-dimensionalizing) , and X22 Is 


9 


approximately -25.3. 
tem is such that x 
meridian. 


By equation (2.5), the orientation of the coordinate sys- 
is in the equatorial p).ane and passes through the Greenwich 


By equations (2.2) and (2.3), 


« = 1/2 0)2 (x2 + y2) (3 3i„2^ . 1) 

r r3 


£22 

♦ — - cos 2 $ cos 2 (X-X22) 


( 2 . 8 ) 


Then 
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f « 


1/2 (1)^ (x2 ♦ y2) 


♦ 


V 

- ♦ 


r- 


— ( €2 ♦ £22 2 ( X - X 22 ) ) 

r3 


(3 C2 ♦ 622 ®®s 2 (X - X22)) 

I • pD V 

The point (s) (x,y,z) = (a,b,c)| where C|b,c are constant are called 

equlllbrluni (stationary) points, and nay be computed by noting that, since 

i 

X X a, 
y * b, 

Z S c. 


y 


z s X s y 


s z s 0, 


SO that equation (2.1) becomes 


$x (®»b|C) - C, 
^y (a,b,c) — 0, 
♦2 (a,b,c) = 0, 


( 2 . 8 ) 


By taking partial dorivatives of ^ in equation (2.7) with respect to x,y, and 
z, respectively, and using the abbreviations 


W (x,y,z) 



— ( E2 ♦ 622 2 ( X - X 22 ) ) » 

r5 


2 £22 

Z (x,y,z) = sin2 (X - X 22 )i 

r5 

r 

Y (x,y,z) s — (3 62 622 2 (X - X 22 ))» 

r5 


(2.9) 


4 
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In light of equations (2.8), it turns out that 


f s (a,b,c) = 


z2 *2 ’ 


« = 


C 5 


(a I b I c ) 

^2 (a,b,c) 


xW ♦ yZ (1 - f (a,b,c) = 0. 

z2y 

yW - xz (1 - — y (a,b,c) * 0. 

- 1 

-Z (o^ - W) + - (2 - 5--) Y (a,b,c) = 0, 

5 p2 


( 2 . 10 ) 


which Is satisfied if all of 


W (a,b,c) s 0 
Z (a,b,c) = 0 


( 2 . 11 ) 


and 


Z = c = 0 

are satisfied. 

By equation (2.9), 


Z (a,b,c) s 0 


gives 


sin 2 (X - X22> = 0 


so that 


TI 3it 

X - X22 = 0 »““* ““f 


( 2 . 12 ) 


5 
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and thus 

1 

cos 2 (X - X 22 ) s 

t 

Note that a and b may be computed from equation (2.5)* 

Also, by equation (2.10) 

j W (a,b,c) = 0 

through equation (2.9) gives 

uj^p = y/p^ (3/p*0(e2 1 £22^ » 

where 

p = r = (a2 

thus enabling p to be computed. (It turns out that p = 42,400 Km.) 
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3.0 LINEARIZATION 


Motion in a neighborhood of the four libration (stationary) points may be 
investigated by Introducing displacements Hi C « libration point 

a,b,c by • 


X * a 4 

y s b 4 Til (3.1) 

Z * C 4 t . 


Consider the Taylor series expansion of the force function, about (a,b,c), 


o o o o 1 0 o 
« (x,y,z) =4) 4 ♦♦yn 4^2^) 4 

o - - 0,0 o o 

4 -♦yyn^)^ ”*• ■►♦xy^n 4$x2i^C ■►♦ysfic) • 

By equation (2.1), if a second order Taylor series approximation is used, 
following equations are obtained 

c - a)n = 'H74.; (C,n,c)'f 

i? 4 = <ty = 4>y 4 (4> 4 4> yyT) 4 

= 4*y 4 74>y (K ,n,C 

c = 4*2 = 4*z4 (<t J5X^ 4 4 zyH + ♦ 22^ ) = ♦ 2 4 V4» z (C 
0 0 0 

Since 4x=4>y=4'z = 0, 


(3.2) 


the 


(3.3) 
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equation ( 3 > 3 ) reduce to 


ff - 2un = ^ ^ xy*^ ^ XT?" " X 

^ ♦ a4 ^ (5,n,C)T (3.^) 

^ ~ ^ zy[* * ^ zy^ * ® z » 

(C S auii , n ♦ a*) 4 * , C) « H<>(C,n, 5 )’^ , a Unear equation in 5 ,n ,C , 


where 


o ^ 0^0 
^xx ^xy ^xz 


.0 *.0 

H® 5 ♦yx ^yy ^yz 


« 


o 

zx 





o 

zz 


ia the mission matrix of ^ at (a|btC), the stationary point. 

By using equations .(2.2) and (2.9) i the elements of H® are shown in appendix 

A to be 


♦ jjjj s aa2 - 3b2 , 

♦yy = ab^ - Ba2, 

® “5 

'^zz = - “ Y 

o « *® 

®xy • ^ 

^xz = 0 = ®zx* 
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wh«r« 

a - 3M/p5 4. (15/P7)(e2 ♦ e^g), 

8 s ♦ HC22/P7, (3.6) 

Y . (2/p5) ( 3^2 1 ^22>- 

Thus, by equations (3.^) and (3.6), 

C - JiWn - Oa2 - 8b2)C ♦ ab(<* B'n 

n ♦ 2w^,= ab (a + 8)C ♦ (Cib^ - Ba2)n (3.7) 

C - (0)2 . Y)C. 


Note that equation (3*7) states that the out-of-plane motion (z-directlon) 
is completely uncoupled from the In-plane motion (x-y plane, or equatorial 
plane). Bond ( ,pp.i2-17) discusses the stability of this motion which 

allows a solution to the differential equation (3.7). 
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H.O QUADRATIE EFFECTS 
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An In section 3.0, consider a third-order Taylor series approxlaatlon to the 
force function given by equation (3.2), i.e., 

f 

, oooolo.o^o., 

♦ (x,y,z) s ♦ ♦ ♦ taft ) ♦ ♦ ♦yyH^ ♦ 

0 0 o^1o-o, 

2^xyCn ♦ 2^xzSC ♦ 2fyjnc) •K^xxj^^ * ♦yyyh^ 

I 

■*• ♦zziiC^ 38xxyC^ * S^xxzC^C ♦ 3^xyyCn^ 


O - 0,0- 

* 3#yyzq% 3*xzzSC^ 3®yzzOC^ 
o 

64>xyzSTic )■►••• • 


o o o 

since = 4>y = 4*2 - using equation (4.1) and substituting into equation 

(2.1), the following equations are obtained: 


5 - 2u)q = $x = ^^xx5 ^xyh ♦xzC) ^^xxx5^ ^xyyH^ 

0,0 o o 

^XzzC*^ 2®xxy^n + 2$xy2^»; + 2$xyzHC) 

= v<tx (c.n.c)'^ (c.n.c) Hx (KfnfO’^, 

. 0 O 0 1 o ® o 

n = 4>y = f4*yx^ ♦yyh + 4*yzC^ ■*■ "I ^^yxxS** * ♦yyyh^ 

0,0 o o 

■*" "^yzzC*^ 2^yxyCn 2^yxzCC ■♦■ 2<^yyzlC 

= 7<i)y (c.n.c)'^ ♦ (?.n,c> hJ (C.n.c)’^. 

o 0 o 1 o , o , 

C = ♦z = '*zxC ♦zyH ^zzC) ■*■ “* ^^zxxC*^ ■*■ ♦zyyH'^ 

0,0 0 o 

■♦■ ®zzzC^ 2<t25xy^n 2®zxzCC ■♦■ ^^gy^nC^ 


(M.2) 
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0 

wxx 

^wxy 

wxz 

0 

wyx 

wyy 

wyz 

,0 

wzx 

^wzy 

^wzz 


Inhere the index variable w s x>ytZ. 

Note that equation (4.2) pass on a quadratic fora not found In equation (3*4). 

By u<)ing equations (2.2) and (2.9) (see the appendix), the following equations 
are obtained (this is similar to the analysis performed in equations (3.4) 
through (3.7)). 


C - 2wf| s 


(Oa2 - Bb2) 5 ♦ ab (a ♦ B)n ♦ 



o ♦ 


— ^(b20 - o2)) 

p2 


K2 


♦la o ♦ — (a2 - 4b2) B - 5b2v + 2 B n2 

2 p2 


1 5 3 

♦ — a (o + ““ Y)(»2 ♦ b Cl 4- (b2 — 4a2) B 

2 p2 p2 

- 5a2v ♦ 2B Cn, 


(4.3) 


1 5 

n + 2w^ = ab (Cl + B)C -4 (Cib2 Ba2)h ♦ -b a ♦ — (b2 ♦ 4a2)B 

2 p2 

3 5 

- 5a2v +20 52 + «b a + -. (a20 - b2v) n? 

2 p2 (4.4) 

1 5 3 

+ - b (Cl + — Y)^2 + a a + — 

2 p2 p2 


ia2 - 4b2)B ~ 5b2v 
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C ■ (4i) 2 « Y)C ♦ (a ♦ -5 Y) (at ♦ bn)C, 

p2 


(*1.5) 


where 


V 


W 7 

-- 1 . (C 2 ♦ €22). 
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In 1 I by using a second-order Taylor series expanded about a stationary 
pointy Bond found the following linear second-order differential equations 
governing motion about the stationary point (see oq.(3.7)): 

C - “ 8t>2)5 + ab <a + 0>n» 

T) - = ab (a + 0 >5 ♦ (otb^ - 0a2>n, 

5 = (*u)2 - Y )^ . 


By using a third-order Taylor series expanded about a stationary point 
been shown here that the following quadratic second-order differential 
governing motion about the stationary point (see eq.i^ 3), (*1.4), and 
(4.5) } : 


, it has 
equations 


C “ = (aa^ - 0b2)5 +ab(a*»-0)n - aa+-- (b% 

- ^ _ 3 (.2 - I|b2)g . ^ 

2 

1 5 , 3 

+ -a(a + --y)^2^5 ^ ^ . 4a2)3 - 


^ I 


n aoC = ab (a + 0) ^ + (ab2 - 0a2)p 

1 3 

+ - b a ♦ — (b^ - 4a2)0 - * 2Q ^2 


3 5 15 

+ - b a + -- ( - b^y ) ri^ + “b(ot + — 

2 ^ ^ ^ 


+ a a + (a^ - 4b2)g - Sb^y + 20 ^ri , 


C = (ti)2 - y) C (a + -- Y> (a^ bn)^ . 
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6.0 CONCLUSIONS 

Note that the first two terms of the ln>plane equations (eqs. (H.3) and 
are thr same as in the linear case (the first two equations of equation (3.7)). 
Also, the first term in the out-of-plane equation (equation (^<.5)) is the oaine 
as In linear case (the last equation of equation (3.7)). The remaining terms 
in equations (i|.3) through (il.5) are quadratic terms. Then, it is necessary to 
go to a third-order expansion of the force function ^ to observe coupling 
between the In-plane and out-of-plane motion, i.e., so that the in-plane and 
out-of-plane motion do not seem to be independent of one another. 

The next phase of study will be to consider only the coupling effects of the 
In-plane motion in the out-of-plane equation (equation (i|.5)). Note that the 
.out-of-plane equation has a perturbation from the nonlinear terms which is of 
the order of the frequency, (j). 
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Figure 1.- The geometry of the rotating system. 
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TEMPORAL CHANGES OF THE GLOBAL REFLECTANCE 
OF A WHEAT FIELD AS A FUNCTION OF 
DAILY SOLA.R IRRADI/ANCE 


Prof. Guy A. Franceschini, Ph.D. 
Department of Meteorology 
Texas A&M University 


Supervisor: Dr. David E. Pitts 

Supporting Research Branch 
Earth Observations Division 

Abstract : Based on in situ measurements of incident and 

reflected solar irradiation over a wheat field, 
daily values of the surface reflectance, a scene 
signature, were determined, for a crop year. 
Diagnoses of these data reveal the character of 
the signature, and its changes with time, crop 
stage, and the magnitude of incident irrad lance. 
The latter varies inversely with cloud cover. 
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1. Introduction 


This is a diagnostic study of the characteristics of thermal-IR 
signatures of a field planted in spring wheat during the 1978 crop 
year. To date, and to my knowledge, the greater effort in dealing 
with crops and satellite-derived data has involved shortwave radiation, 
i.e., solar radiation, of which 99% is in wavelengths shorter than 
4.0 iJm. This concentration of effort is primarily a consequence of the 
large variations in optical properties, viz., reflectance, associated 
with soil and crops. Such variations are wavelength dependent, and are 
a function of crop geometry, stage of maturity, and stressing which is 
associated with atmospheric and edaphic conditions. In addition, the 
nature, condition, and moisture of the exposed soil also influence the 
surface reflectance to a great degree. Hence, spectral measurements of 
surface-reflected shortwave radiation make it possible to characterize 
soil moisture, and crop extent and conditions, as well as to estimate 
ultimate possible crop yield. 

In considerations of thermal infrared radiation, i.e., that within 
the terrestrial waveband which extends from 3.0 to 80 pm, we essentially 
are dealing witli a thermodynamic coordinate, temperature, of the target. 
Ideally, most terrestrial surfaces behave nearly as 'black bodies' in 
this longwave spectral domain. However, except for limited wavebands, 
the atmosphere also behaves as a black body primarily because of carbon 
dioxide, water vapor, methane, nitrogen-oxygen molecules, and clouds. 
Consequently, to view the surface with satellite-mounted sensors on 
essentially cloud-free days, we must select a spectral region over which 
the atmosphere is nearly transparent, e.g. , 8-13 pm. Unfortunately, the 
atmosphere, even in this waveband, behaves as a somewhat 'dirty window' 
having a transmittance less than one. Nonetheless, meteorological sat- 
ellites effectively and routinely use this waveband to delineate surface 
(land and ocean) and cloud top temperatures, i.e., corditions associated 
with we:|l defined horizontal temperature gradients. 

A 4uestion which should be asked is: Can such infrared spectral 
information be used tc characterize croplands? In this study, we have 
examined spectral radiance data obtained with sensors mounted on an 
aircraft flying at a low altitude over an area planted in spring wheat. 
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OF 


Because atmospheric interference is minimized due to the low altitude, 
the data may be considered a ground truth for satellite measurements. 


2 . Theoretical considerations 

I 

a. The nature of radiation 

According to Planck's law, the monochromatic radiance, , from a 
blackbody may be expressed as a function of wavelength, X ym, and the 
temperature, T °K. For the units chosen, 

= (ci/A^) (exp ( c 2 /XT) - l) (yW/cm^ sr ym) , (1) 

where cj = 1.19096 x lO^'^; and C 2 = 1.43879 x 10‘^. For given 

values of and A, the effective or radiation temperature of a target 
may be computed using (1) . The Planck curve, as given by (1) , for a 
blackbody at a temperature of 300®K is shown below. 
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As seon in the figure, the modal wavelength, i.e., the wavelength 
of maximum radiance, is near 9.66 pm. The total waveband of interest in 
this study, 10.4 - 12.6 pm as indicated in the figure, involves nearly 
14% of the total radiance. In addition, within this spectral band, the 
value of varies inversely with wavelength. Observed values in this 
band range from 750 to 1130 SU (1 SU = 1 Standard Unit = 1 pW/cm^sr pm) 
with associated effective temperatures from 290°K to 310“K, respectively, 
over the limits of the total waveband. 

b. The energy budget 

Changes in temperature of a target depend upon its net gains and 
losses of energy. A physically sound approach to the problem is given 
by the energy equation which is a statement of the boundary condition 
at the interface between the atmosphere and its subjacent volume. The 
condition states that the energy flux is continuous across the interface, 
i.e., the net flux on both sides of the boundary are equal. In terms of 
energy fluxes, and with the convention that fluxes toward the interface, 
on both sides, are positive, the energy equation may be symbolically 
expressed as 


N + Q L + H = 0. (2) 

In (2) , the first three terms refer to fluxes on the air side of the 
interface, and the last term, H, refers to the flux on the opposite 
side, i.e., the side of the subjacent volume, o.g., soil and/or plants. 
Individ t.AC.4 1 terms are: N, the net radiation, both short- and longwave; 

Q, the sensible heat flux; L, the latent heat flux associated with the 
phase change of water; and H, essentially the sensible heat flux below 
the interface. 

Individual terms of (2) may qualitatively be expresses as • 

) 

net radiation, 

N = (1 - r)I + (Fa - F) , 

where I is the solar irradiance, r is the associated surface reflectance, 
and the longwave fluxes. Fa and F, are 
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Fjj “ T^** (downward from the atmosphere) 

F oc (upward from the surface target) ; 

the sensible heat flux, 

Q “ (Tg - T)v", 

in which v is the wind speed, and n is an empirically derived exponent 
often taken as n * 2? 

the latent heat flux which may produce surface evaporative cooling, 

L “ (6g - e)v*^, 

where e^ is the vapor pressure of the air, and e is the vapor tension 
of the target surface; and finally, 

the sub-interface flux, 

H oc (T - Tg) , 

in which T^, is a near-surface temperature of the interior of the target. 
Consequently, the energy equation shows the interrelations and feedback 
mechanisms between the transfer processes. More importantly, for our 
purpose, it indicates the dependence of target temperature, T, on the 
atmospheric elements, i.e., temperature, wind speed, and moisture, as 
well as the moisture state of the target, both soil and crop. 

In particular, the heat transfer, H, within crop and soil, which 
must balance all the other in-air fluxes, represents the energy which 
enters the surface, and is consequently responsible for target temper- 
ature changes. For a thin near-surface layer of the target, extending 
to a distance zj from the interface, 

(9T/3t) “ (Hi - H) , (3) 

where Hi represents the flux across the level Zi, and depends on the 
temperature gradient at that interior level. In (3) , a prognostic 
equatio^, we see the dependence of the rate of surface temperature 
change on the atmospheric fluxes, since, from (2), H = - (N Q + L) . 

c. Areal averaging 

In principle, (2) applies to each different element of the surface. 
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crop and .soil, separately. However, for this study, we will treat the 
target arc. u, as a collective system with representative temperature 
T, as an area average defined by the mean-value theorem: 

m 

T » (l/u)/Tdu = (1/u) I ^ u . ( 4 ) 

u 1*1 ' ' 

In (4), T represents the variable temperature of in target elements, 
Since a target area usually consists of soil and crop, a differentia- 
tion of the contribution of each to the upward flux of longwave radia- 
tion, F, is important. If c represents the fractional area of the 
target which is covered by the crop, i.e. , ground cover, then in a 
manner analogous to (4), area-weighting each contribution leads to 

F = cF^ + (1 - c)Fg, (5) 

vfhere subscripts c and s refer to crop and soil, respectiv’ely , and 

F <x T** a 



where is the measured monochromatic radiance of the collective 
target, as given by (1) . 

Since airborne sensor measurements slur the successive readings 
over different wavelengtli intervals, a further area-weighting of the 
individual readings is possible. For example, if represents the 
fractional target-area overlap successive readings, a different 
space-average value of target temperature, <T^>, at time t, may be 
defined as 

<Tt> = (At-jTt-j + Tt + At+kTt+k)/(At_j + 1 + (6) 

I 

where: the T values are obtained from radiance data assigned to the 
central wavelength of the interval; and j and k represent the interval 
of time, respectively, between values preceding and following that at 
time t. The values of A^ in (6), pertinent to this study, are treated 



6 


opjdnv. 

Of poo^i 



quality 


in the section on target geometry. However, due to time limitations, 
such area-averages of target temperatures were not calculated. 


3. Observations 


There were five observation periods, each on different days of the 
1978 crop year. Measurements were made of: spectral radiance over a 
time interval of 48 - 56 s, once each day; atmospheric elements; soil 
field moisture; and certain agronomic information. Each is discussed 
below. 

a. Ancillary data 

Supporting simultaneous atmospheric measurements of wind, pressure, 
temperature, relative humidity, and cloud cover were made during each 
observation period. Agronomic data included: ground cover (including 
weed coverage) , maturity stage, plant height, and dates of planting 
and harvesting. Field moisture information consisted of quaiitive 
estimates, e.g. , dry, damp, and wet. The latter condition, although 
not reported, was assumed after viewing photographs taken during the 
observation periods. Information on soil type was not available. 

b. Radiance data 

The radiance sensor of the airborne Field Spectrometer System, 

FSS, has a field of view of 2u.00®, and a scan rate of 1 s, i.e., all 
channels were slur-sc^uned sequentially every second. In all, there 
are 42 channels from 10.4 to 12.6 pm, and each covers an interval of 
0.05 pm. Average values of monochromatic radiance were available for 
five wavebands involving different numbers of channels, hence, differ- 
ent timi intervals, and different areas viewed. Details are tabulated 
;i - 

below, where AX is the w; '.ength interval, and X is the median or 

central wavelength of f iterval. 



7 


OF PCGu 


Band 

AX 

(ym) 

No. Channels 

X 

(ym) 

Scan Time 
(s) 

1 

10.40 - 

• 10.80 

8 

10,50 

0.191 

2 

10.80 - 

■ 11.30 

10 

11.05 

0.238 

3 

11.30 - 

■ 11.80 

10 

11.55 

0. 238 

4 

11.80 ■ 

■ 12.30 

10 

12.05 

0.238 

5 

12.30 - 

■ 12.50 

4 

12.40 

0.095 


Target geometry 

Flight altitude was 60.96 m (200 ft). Consequently, with a 20.00® 
field of view, the target diameter was 21.5 m. At a flight speed of 
96.52 km/h (60 mi/h) , the distance between initiation of successive 
1-s scans was 26.82 m, and the average length of the area scanned was 
approximately 1.6 km. Hence, the area scanned was approximately 3 ha. 

Since the distance t 'tween successive scans of any specific wave- 
length band is greater than the target diameter, there is no overlap of 
related targets. Howsver, target areas of sequential but different 
wavebands do overlap. In general, if we assign the average radiance 
value of a waveband to the central time, hence, position, of that 
waveband, the overlap area of two successive and related targets 
may be expressed as 

A = {ird^ - 2L d^ - - (2d? ) sin’l (L /d) }/4 , (7) 

\ \ \ \ 

where d is the target diameter, and is the variable distance between 
successive target centers. Values of and the related fractional 
target overlajj, A^ associated with successive average radiance readings 
are given below. Fractional target areas from (7) may be used in (6) 
to evalliate <T^>. Again, due to time constraints, it was not possible 
to complete the evaluations. 
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B 


Bands 

L 

i 

(m) 

A 

1 

(%) 

1 to 2 

5.747 

74.58 

2 to 3 

6.386 

71.77 

3 to 4 

6.386 

71.77 

4 to 5 

4.470 

80.19 

5 to 1 

3.832 

83.01 


5. The site 


The observation site is located in the wheat drylands of Hand Co. , 
SD, in field No. 281, having an area of 128 ha (316 acres). Specific 
location is: 44° 30’ N, and 99° 0* W. The field was planted in spring 
wheat, variety Fortuna HRS, on 26 April 1978, using the R-and-R method 
with 7 in. between rows. The crop was harvested before the last period 
of observation on 16 August 1978. Plants had attained a height of .81m 
at which time the ground cover was 30%. 


6. Source of data 


Data were obtained from the Laboratory for Applications of Remote 
Sensing (LARS) of Purdue University, West Lafayette, ID. LARS was 
responsibre for taking the observations, as well as determining the 
average radiance values - 


7 . Discussion 

I 

Data are presented on the following five figures, one for each 
observation period. Each contains a time reference, all available 
ancillary data, and graphs of average monochromatic radiance values 
for each of the five wavebands. The latter are plotted versus time 
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and approximate travel distance. Included in the i.jraphs are scales 
for temperature, one for j»ach waveband. These effective temperatures 
were calculated using the F’lanck relation as given by (1). Local mean 
solar time for the site is ^approximately 6.6 h earlier in the day than 
the indicated GMT. 


a. Observation number 014, 780515, 17:13 GMT 

On this day, 19 days after planting, the crop had not emerged, but 
there was a reported ground cov«»r of in'*, due to weeds. The wind was 
strong from the bSL during this late morning period, and the relative 
humidity was only 40%. Both conditions would have tended to enhance the 
sensible and latent heat fluxes, but the latter must have been minimal 
since tlie soil was reported as being dry. 

Variations of radianc'e valu<>s for all bands were in phase, and 
showed a range of 61 to 87 8iJ with B"5 (Band 5, 12.40 \im) and B-1 
(Band 1, 10.60 urn), respectively, which correspond to ranges of T, 
the effective temperature, of 5.1° and 5,5°K. The average range, AT, 
for all bands was 5.3°K. 

In considering the average values of T, i.e., T, for all bands, 
shown by horizontal lines in the figure, the maximum differe'nce, dT, 
between values was 2.5°K, and was found between data for B~1 and B-5. 

B-1 had the highest, D-5 the lowest, viz., 307. 9°K end 305. 4°K, respec- 
tively. A similar result appeared in all the {leriods, and may be due 
to tlu^ method of determining tlie average radiances. For bands 2, 3, and 
4, dT < 0.2°K. Of the three, B-2 had the highest temperature, 306. 5°K. 

For all bands, the values of T were higher than air temperature, 
295. 1°K. Differences ranged from 10.3°K to 12.8°K. Although excessive, 
such vertical differences are not surprising for the following reasons: 
the 1% cloud coverage pemitted maximum insolational heating; the soil 
was dry and essentially bare, thu.s, evaporative cooling was a minimum; 
the time was during the U|)swinq of the diurnal temjierature cycle when 
air tem|ieraturo lags that of t''e surface; and the strong winds induced 
sufficient tui'bulence to concentrate the vertical temperature gradient, 
as shown by the large differences, in the near-surface layer of the 
atmosphere . 


OWGINA*- 
OF POOR quality 
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MONCX:ilROMATlC RAOIANCF,, , AND EFFECTIVE TEMFERATURE, T 

DaUVTimo: 7H0 'j 15/17 : li GMT Observation No. 01^ 

Air Temperature, 295.1'’K- Wind, 1 ;)!>■’ 32 km/h; Relative Humidity, 40% 
(iiound Cover, 10 % (weeds); Maturity Stage, no emergence; PLHT, 0 dm 

Field Moisture, dry. 
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b. Observation number 094, 780002, 15:55 GMT 

The maturity stage of the crop, on this 37th day after planting, 
was reported as tillering prebud. Plant height was 0.1 m, and ground 
cover was 10% wj th no weeds present. Although the field moisture was 
reported as dry, antecedent precipitation and field photographs taken 
during the period indicated a possible wet condition. The wind was 
light to moderate from the west; and relative humidity was 60% which 
was more in keeping with a non-dry surface. Cloud coverage was only 
5%. 

As with OBNO 014, radiance variations were in phase, and all mini- 
mum values corresponded to the wet and dark areas that were seen in the 
photographs. Ranges were excessive, the largest encountered in this 
study, and varied from 80 to 110 SU for B-5 and B-1, respectively. The 
associated ranges in T were 7.3"K and 7.7“K, with an average for all 
bands of AT = 7.6'’K. 

Values of dT, the interband maximum difference in average values 
of temperature, was only 1.8“K, and was found between data for B-1 
(T = 295. 1°K) and B-5 (T = 293. 3°K). For the central bands, 2, 3, and 4, 
dT _< 0.6°K; and for B-2, T = 294. 6®K, which was highest of the three. 

Althougli all values of T were higher than air temperature, 290. 9“K, 
the differenr o were much smaller than in OBNO 014, and ranged from a 
minimum of 2.4°K to a maximum of 4.2°K. The lower temperatures were 
presumably a consequence, in part, of evaporative cooling of the wet 
surface. Also, advection associated with the westerly wind may have 
contributed. The latter two comments are speculative, since the true 
situation was not investigated. 



ORIGINAL PAGE 13 
OF POOR QUALITY 

MONOCHROMATIC l^VDIANCE, , AND EFFECTIVE TEMPERATURE ,T 

Date/Time; 780602/15:55 GMT Observation No, 094 

Air Temperature, 290. 9°K; Wind, 267° 14 km/h; Relative Humidity, 60% 
Ground Cover, 10%; Maturity Stage, tillering prebud; PLHT, 1.0 dm 

Field Moisture, wet. 
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c. Observation number 243, 780709, 19:01 GMT 

This period was shortly after local mean solar noon on the 74th 
day after planting. The crop was fully headed; plant height was 0.76 m; 
and the ground cover was 30% with no weeds present on the dry field. 

The wind was light to moderate from the north; and the atmospheric 
relative humidity was 50%. Sky cover was 10%, but should not have 
interfered appreciably with solar heating. 

As with the two preceding periods, the radiance variations were in 
phase, but showed smaller ranges. These were from 38 to 73 SU, again 
associated with B-5 and B-1, respectively. Corresponding ranges of 
temperature, AT, were 3.4®K and 4.8°K, respectively. However, the 
minimum range, viz., AT = 3.2°K, was with B— 2, not B-5 as in the two 
preceding observation periods. The average range of temperature for 
all bands was 3.9°K. 

The interband maximum difference between average temperatures, dT, 
was 2.3°K between B-1 with T = 300. 1°K, and B-5 with T = 297. 8“K. As 
in the other periods, minimum differences were found between bands 2,3, 
and 4, viz., dT < 0.5°K; and of the three, B-2 had the highest T, viz., 
299. 4°K, as in OBNO 014. 

The average effective temperature, T, associated with each of the 
hands was again higher than the air temperature, 296. 5°K; and the range 
of differences was from 1.3°K to 3.6°K. These differences were smaller 
than those found in the preceding two periods. 
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MONOCHROMATIC RADIANCE, , AND EFFECTIVE TEMPERATURE, T 

Date/Time: 780709/19:01 GMT Observation No. 243 

Air Temperature, 296.5'>K; Wind, 005° 16 km/h; Relative Humidity, 50% 
Ground Cover, 30%; Maturity Stage, fully headed; Plant Height, 7.6dm 

Field Moisture, dry. 
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d. Observation number 332, 780726, 15:10 GMT 

On this day, the 91st after planting, the crop, at a height of 
0.81 m, was beginning to ripen, and covered 30% of the ground. The 
field was dry and weedless. A very light westerly wind, relative 
humidity of 54%, and a clear sky characterized the atmosphere. 

Radiance variations for all bands were again in phase; and ranges 
varied from 55 to 84 SU, with associated AT values of 4.4®K to 5.5°K, 
respectively. However, although the maxima were found with B-1, the 
minima were in the B-4 data. Nevertheless, values for B-5 were quite 
close, viz., 5 7 SU and 5.0°K. The average temperature range values 
for all bands was 5.0®K. The light wind, with a minimum of turbulent 
mixing, may have been responsible for these larger values, i.e., as 
compared with those of OBNO 243 which had the same percentage value 
of ground cover. 

Between B-1 and B-5, the interband maximum difference in average 
temperature, dT, was 2.3®k, as in the preceding period. Again, the 
highest temperature (302. 0°K) was associated with B-1; the lowest 
(299. 7®K) with B-5. For bands 2,3, and 4, dT £ 0.5°K, with B-2 
having the highest of the three, viz., T = 301. 2"K. 

For each band, the average temperature was higher than the air 
temperature, 299. 3 °K; and differences, which ranged from 0.4'*K to 
2.7®K, were smaller than in the preceding three observation periods. 
This may have been a consequence of the very light winds which may 
have minimized the near-surface vertical temperature gradient because 
of reduced mechanical mixing. 
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e. Observation numlL5er 423, 78081fi, 15:11 GMT 

At this time, 112 days after planting, the wheat crop had been 
harvested, and the remaining stubble was 0.25 m high, having a ground 
cover of 10%. The field was damp and weedless, but very likely had 
some scattered remains of the harvested plants. Wind was moderate and 
southerly, sky was clear, and the relative humidity was 45%. 

Except for minor deviations, the variations of radiance values in 
all bands were in phase, and showed the smallest ranges, viz., 22 SU to 
45 SU, of all the periods in the crop year. The associated temperature 
ranges were 1.9°k for B-4, and 3.1“K for B-1, respectively. The average 
AT for all bands, as well as for B-5, was 2.5°K, which was the smallest 
of all periods. 

The maximum interband difference of average temperature, between 
B-1 (T = 296, 9“K) and B-5 (T = 295. 3®K) was 1.6°K, and consequently the 
smallest of all periods. For bands 2,3, and 4, dT ^ 0.3“K, This value 
compares favorably with the 0.2®K associated with the period before the 
crop emerged when the field was sparsely covered with weeds, i.e., OBNO 
014. 

Unlike all other cases, the average effective temperature of the 
surface, associated with all bands, was lower than the temperature of 
the air, which was 298. 2®K. Differences ranged from 1.3®K with B-1 to 
2.9°K with B-5. Since the soil was damp, we might surmise that this 
may have been due to evaporative cooling associated with the moderate 
southerly wind. 
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8 . Summary 


P^^tinent details of this data diagnosis are summarized in the 
following figures, in which data for each observation period are 
aligned vertically. Fig. A presents field moisture and agronomic 
features, i.e., ground cover and plant height as a function of the 
number of days after planting. Fig. B shows temperatures of surface 
(only the extreme values) and air, as well as AT and maximum dT values 
as a function of time of day. The last. Fig. c, presents atmospheric 
data, i.e. , relative humidity and wind (numbers near the plotted wind 
speed points represent direction in degrees), plotted versus date. 
Since the characteristics previously discussed are evident if the 
figures, no further comments are needed. 
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FIGURE CAPTIONS 


Sponge moisture variable. ^ 

Texas sponge values: 1980 (Appendix) 

Normal annual Texas sponge values " 

Texas sponge values: April, 1980 " 

Normal Texas sponge values: April " 

Texas sponge values: July, 1980 " 

Normal Texas sponge values: July " 

Texas sponge values : October, 1980 " 

Normal Texas sponge values: October '' 

Example of a band of 5 NOAA/AVHRR Scanlines and the 5x5 pixel sampling 
grids along the Texas transect. 

The normalized difference (ND) regression for 3 sample dates in 1980 
along Texas transect. 

The mean normalized difference (ND) regression line for 3 dates in 
'i980 along Texas transect. 

Long-term (1941-1970) sponge normal values and regression from 25 
stations on Texas transect. 

Normalized difference (ND) values vs. sponge indices for 12 sites 
(3 dates) along- Texas transect. 

Mean normalized difference (NO) regression (see Figure 8) vs. sponge 
line (see Figure 9) on the Te^as transect. 

Vegetation - Sponge Index (V5I) values along the Texas transect and 
suggested limits of the vegetation regions. 


MAP CAPTIONS 



1. 8-inch sponge normals for Texas: 1941-1970 J 

^Except as noted in Table 2 (Appendix). 

2. Normal annual average of 8-inch sponge - State of Texas. 

3. Normal January average of 8-inch sponge - State of Texas. 

4. Normal average April 8-inch sponge - State of Texas. 

5. Normal July average 8-inch sponge - State of Texas. 

b. Normal October average 8-inch sponge - State of Texas. 

7. 12 sample locations on Texas transect. Letters (A-D) refer to major 
vegetation regions along transect.^ 

Viodified from Gould (1975). 

8. Location of stations used to calculate sponge values along Texas 
transect. 

9. Mean 8-inch sponge values along Texas transect: 1980 
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I. OVERVIEW 


Phase 1 of this investigation has several objectives: (1) to test a new 

experimental climatological model /variable termed the sponge for potential 
biogeographical, ecological, and climatological applications (the sponge 
is a measure of moisture availability based on daily temperature maxima 
and minima, and precipitation); (2) to investigate the feasibil ity of 
utilizing NOAA/AVHRR meteorological satellite data for vegetation classi- 
fication; and (3) to initiate a vegetation gradient model that utilizes 
climatological (i.e., sponge), biological, and NOAA data that is ultimately 
applicable to global vegetation stratification and monitoring. 


To accomplish the initial objective, mean monthly and annual sponge values 
are calculated for 75 Texas locations for the "normal" period of 1941-70. 
Similar values are also computed for approximately 25 stations along an 
east-west transect across Texas for 1979 and 1980. Results suggest that as 
a generalized climatic index, sponge’s simplicity and sensitivity make it 
particularly appropriate for trans-regional biogeographic studies. 


The latter two objectives were approached by acquiring vegetation, clima- 
tological (sponge), and AVHRR pixel data (channels 1 and 2) for 12 locations 
along the east-west .Texas gradient. The normalized difference (NO) values 
.for the AVHRR data when plotted against the vegetation characteristics (bio- 
mass, net productivity, leaf area) and the sponge values suggest that a 
multivariate gradient model incorporating AVHRR and sponge data may indeed 
be useful in global vegetation analysis. 


II. THE SPOMGE: A NEW ECO-CLIWATOLOGICAL VARIABLE 


The planetary distribution of types of natural vegetation is largely a 
function of climate, most especially of spatial variations in energy and 
moisture budgets. As a general rule, climate is recognized as the pre- 
eminent contro'' of natural vegetation at subcontinental-to global-scales. 

For small areas, geologic, pedologic and other local factors may dominate 
for long, even indefinite, periods of time. 

Natural scientists such as ecologists and biogeographers have attempted to 
utilize climatic measures and indices in surveys, stratifications and 
classifications of natural vegetation. Typically, it is assumed that most 
major ecoregion and native vegetation-region "boundaries" actually represent 
climatic discontinuities or "breaks". This assumption, while not entirely 
valid, probably is reasonably accurate. 

A key problem for such scholars has been the determination of a simple but 
accurate climatic "index" (or indices) which would reflect the primary 
spatial variations in moisture and energy balances and, hence, could be 
applied to the classification of native vegetation. Climate itself, a 
complexly synergistic synthesis of many variables, does not readily succumb 
to quantitative cl assi ficaiion. [For example, what is the "real" boundary 
between a desert and a (semi-arid) steppe? In fact, of course, there is no 
abrupt statistical limit but rather a gradual transition from one type into 
another.] 

Identifying and quantifying climatic factors which explain the distribution 
of natural vegetation is even more challenging. Temperature and precipitation 
by themselves are poor descriptions of climate and, hence, explainers of 
vegetation distributions (Mather and Yoshioka, 1968). 

It is therefore essential to select and/or develop climatic Indices which 
influence vegetation growth and development. Most recent approaches to this 
problem emphasize the importance of moisture availability (e.g., surpluses 
versus deficits) and, more specifically, evapotranspiration. Unfortunately , 
evapotranspi ration is measured at very few places, and even evaporation 


itself is not widely monitored. Numerous models have been developed which 
estimate evaporation from measurements of air temperature, average wind 
speed, and net radiation (Penman, 1948; Jensen ^ a][, 1970), but their 
usefulness is constrained by the sparcity of stations which record solar 
radiation and wind speed. 

Simpler evapotranspiration schemes which require only air temperature and 
precipitation -- both commonly measured around the world -- have been de- 
vised (e.g., Thornthwaite, 1948; Griffiths, 1964; Moe, 1965; Trenchard, 
1976). Thornthwaite 's classification made use of mean monthly temperature 
and precipitation values to generate a moisture index. Because of its 
relative simplicity and accuracy, Thornthwa Ite's approach has been widely 
adopted (e.g., it is used to calculate the USDA's Crop Moisture Index). 
However, its use of average monthly temperatures somewhat limits its 
sensitivity to variations in cont’nentality and altitude. 

An alternate method of relating climate to vegetation is that of multi- 
variate discriminant analysis of climatic variables to determine their 
relative influence in a particular ecoregion (e.g., Biogeocl imatic Units 
of Vancouver Island , Klinka and Nuszdorfer, 1979). While very accurate 
for detailed, site-specific studies, the resulting multiple regression 
equations tend to be (a) cumbersome and lengthy, and (b). less applicable 
to regional and global-scale vegetation classifications. 

Thus climatologists, geographers , and ecologists interested in large-area 
comparisons have found themselves forced to choose between an approach 
which stratifies climate somewhat too broadly (e.g., Thornthwaite's) and 
another which "hides the forest for the trees", viz., too much emphasis 
on detail (e.g., Klinka and Nuszdorfer). 

Recently, Trenchard and Artley (1981) developed a new climatological/ 
meteorological variable whose simple form, minimal data requirements and 
accuracy make it an ideal candidate for application to meso- and macro- 
scale biogeographical , agrocl imatological , and ecological investigations. 


3 


This hypothetical medium is termed the sponge (see Figure 1). Sponge' 
rationale is summarized as follows (Trenchard and Artley): 

We desired a simple moisture variable with a sound 
physical basis that used common meteorological variables, 
was suitable over a broad range of climates, and appli- 
cable to a single station. The result v/as named sponge. 

Sponge is described as a simple medium with 8 inches 
of water holding capacity which is initialized half-full 
of water on 1 January.* Each day, in accordance with 
the hydrologic cycle, water is added to the medium from 
precipitation and lost through evaporation. Precipitation 
(both liquid and frozen) is added at the full amount until 
the layer is saturated. It is this sponoe like behavior 
which gives the variable its name. Any additional pre- 
cipitation is assumed to be lost as run-off or drainage. 
Evaporation occurs at a fraction of the Class A Pan rate, 
the exact proportion being the ration of the current 
contents to the total capacity of the sponge. Either 
actual or estimated evaporation pan values may be used. 

The daily contents of the sponge are defined as: 

5,- ' 5(-l * ^ *Si.l/CAP) 


Where: 

= Sponge contents on day i, in inches. 

P^. = Precipitation on day i, in inches. 

E = Actual or estimated pan evaporation in inches on day i. 
CAP = Sponge capacity in inches 
and 0 £ S. £ CAP 

When evaporation pan measurements are not available, they 
may be estimated with a divisor of 30 days to convert the 
evaporation function to a daily value. 

+ P._(EP(TX,. , TN.)*S._i/CAP*30) 


*A1 ternati vely , the final valu.» of the previous year may 
be used as an initial value, and the capacity may be varied 
for a particular region. 



SPONGE MOISTURE VARIABLE 




Where: 

SP = Pan evaporation function 
TX^ = Maximum temperature on day i . 

TN. = Minimum temperature on day i. 

Because of its simple data requirements {daily pre- 
cipitation and evaporation estimated form maximum and 
minimum temperatures), the sponge can be calculated at 
any temperature precipitation observation station. 

Long-term (1941-1970) sponge "normals" (average values) were recently cal- 
culated for all first-order meteorological stations in the counterminous 
USA (Trenchard, 1981). Some representative values are presented in Table 1. 
It is interesting to note that mean annual sponge values were found to range 
between Yuma, Arizona (0.20") and Mt. Washington, New Hampshire (8.00", or 
absolute sponge capacity). 

Sponge is concluded to present a meaningful measure of areal "environmental 
moistness". As a generalized climatic index, sponge's simplicity and 
sensitivity make it particularly appropriate for transregional biogeograph- 
ical studies (e.g., a large-area vegetation classification). 







J 

Mean January 

Mean July 

Mean Annual 

station 

Latitude (®N) 

Longitude {®W) 

Sponge Value (”) 

Sponge Value (’*) 

Sponge Value (") 

Mt. Washington, N.H. 

44.27 

•71.30 

8.00 

8.00 

8.00 

Stampede, Pass, Wash. 

47.28 

121.33 

8.00 

6.78 

7.49 

Caribou, ME. 

46.87 

68.02 

8.00 

5.47 

6.89 

Duluth, Minn. 

46.83 

92.18 

7.97 

6.01 

6.73 

Pensacola, Fla. 

30.47 

87.20 

6.98 

6.11 

6.41 

• 

Seattle, Wash. 

47.45 

122.30 

8.00 

2.98 

5.67 

Cincinnati, Ohio 

39.05 

84.67 

6.89 

4.00 

5.24 

Chicago, 111. 

. 41.98 

87.90 

5.85 

4.32 

5.04 







Key West, Fla. 

24.55 

81.75 

3.77 

4.66 

4.50 -1 

Houston, TX. 

29.97 

95.35 

5.56 

3.39 

4.27 gp 

San Francisco, Cal. 

' 37.62 

122.33 

6.17 

1.24 

3.35 €5 

2 & 

Dallas, TX. 

32.90 

97.03 

3.26 

2.13 

3.06 5a 

Walla 'Walla, Wash. 

46.03 

118.33 

4.65 

1.14 

2.79 

Salt Lake City, UT. 

40.78 

111.9b 

3.03 

0.99 

2.14 

Goodland, Kan. 

, 39.37 

101.70 

1.24 

2.15 

1.60 

San Diego, Cal. 

32.73 

117.17 

2.29 

0.60 

1.37 

El Paso, TX. 

,31.80 

106.40 

0.76 

0.67 

0.64 

Yuma, Ariz. 

32.67 

114.60 

0.42 

0.04 

0.20 


Table 1. 

Selected USA 8- Inch Sponge Normals: 1941- 

-1970.^ 
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III. THE SPONGE VARIABLC AS APPLIED TO MOISTURE AND VEGETATION 
GRADIENTS IN TEXAS 


There exists in Texas the most pronounced continuous, non-orgraphic, 
intrastate climatological gradient found anywhere in the United States. 

At least four distinct,, first-order climatic types occur within the state 
(humid subtropical, tropical steppe, tropical desert, and mid-latitude 
steppe), with many more important subtypes (e.g., subtropical subhumid). 

In particular, there Is an extraordinarily steep east-west moisture 
gradient, ranging from very humid .in southeastern Texas (mean annual 
precipitation >50") to true desert in far western Texas (average yearly 
precipitation <3.0"). This gradient strongly influences ecological 
patterns, and virtually controls the regional distribution of natural 
vegetation. 

Texas is, then, an excellent natural "laboratory" to test the responsive- 
ness and usefulness of the sponge variable (e.g., with respect to the 
classification of natural vegetation utilizing satellite data). With this 
in mind, mean sponge values were calculated for various Texas locations in 
order to address these questions: 

1. Does use of the sponge portray the distribution of climates 
(especially moisture regions) in Texas better than, say, 
precipitation alone? 

2. If sponge accurately reflects the climates of Texas, can these 
values be meaningfully correlated with vegetation-index ("green- 
ness") values as measured from snace by NOAA meteorological 
satellites (see the discussion of these indices later in this 
report)? 

3. If the answer to question (2) is affirmative, can a combination 
of sponge and satellite-derived greenness indices be used to 
classify the major natural vegetation regions/types of the 
state? If so, it night well prove feasible to utilize this 
methodology for other large-area and even global-scale vegetation 
surveys and classifications. 
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Climatic Strata And Gradients In Texas Using The Sponge; Long-Term "Mormals*' . 

In order to assess sponge's potential usefulness as a climatic index in Texas, 
long-term annual and monthly sponge "normal" (1Q41-1970 mean) values were com- 
puted for 75 locations widely distributed throughout the state (see Map 1). 
Sponge values were obtained by utilizing (1) the formula presented in section 
1 of this report, and (2) mean monthly temperature maxima and minima, and pre- 
cipitation, as compiled by t^ U.S. National Weather Service (NOAA)*. 

The results of these computations are illustrated in Maps 2-6, while Tables 
2a-2e (Appendix) present average normal (1941-70) monthly and annual precip- 
itation and sponge values for each ‘of the 75 locations. They are largely 
self-evident, but several of the more intriguing aspects should be briefly 
addressed. 

Mean annual sponge values are greatest in the southeast (e.g., 5.2" at 
Beaumont) and decrease continuously to lows in the westernmost quadrant 
(0.49”) at Presidio in the Chihuahuan Desert). This is virtually identical 
to the pattern or average annual precipitation. However, sponge appears to 
more accurately portray (1) seasonal moisture changes across the state, and 
(2) the magnitude of differences in the relative moistness of the varous 
parts of Texas than does either precipitation or potential evapotranspiration 
(estimated by Thornthwaite's method; see Table 3). 

With respect to season, Figures 2-5 (Appendix) and Tables 2 (Appendix) and 3 
clearly indicate that, for vegetative activity, winter and spring are the 
wet seasons in East Texas while summer and fall are the moistest periods in 
West Texas -- quite unlike the seasonal distribution of precipitation alone, 
which is greatest in the summer throughout the state. (Thornthwaite's "Index 
of Moisture" would also reveal this aspect but, because of its reliance on 
mean monthly temperatures, with less spatial sensitivity than sponge; Carter 
and Mather, 1 966) . 


♦Daily data were simulated from monthly mean temperature maxima, minima , and 
precipitation for each station for 1941-70 using a series of harmonic trans- 
formations. For 1979 and 1980, actual daily data were used. 
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NORMAL ANNUAL AVERAGE OF 8-INCH SPONGE - 

STATE OF TEXAS 

(PERIOD 1941-70, EXCEPT AS NOTED IN TABLE 2) 
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NORMAL JANUARY AVERAGE 8-INCH SPONGE - 

STATE OF TEXAS 

(PERIOD 1941-70, EXCEPT AS NOTED IN TABLE 2) 
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Map 4. 

NORMAL AVERAGE APRIL 8-INCH SPONGE - 
STATE OF TEXAS 

(PERIOD 1941-70, EXCEPT AS NOTEU IN TABLE 2) 


102 ° 100 ° 
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Map 5. 

normal JULY AVERAGE S-INCH SPONGE - 
STATE OF TEXAS 

period 1941-70. EXCEPT AS NOTED IN TABLE 2) 
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Map 6, 


NORMAL OCTOBER AVERAGE 8-^INCH SPONGE - 

STATE OF TEXAS 

(PERIOD 1941-70, EXCEPT AS NOTED IN TABLE 2) 
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Potential 

Actual 


^ station 

Period ! 

Precipitation 

Evapotranspiration 

Fvaootranspi ration 

Sponge 

1 Paso/La Tuna 

Annual 

8.06 

38.40 

8.72 

0.59 

1 Paso/La Tuna 

Januatiy 

0.41 

OAf 

1 

0.44 

0.71 

1 Paso/La Tuna r 

April 

0.15 

2.72 

0.32 

0. 35 

1 Paso/La Tuna 

July 

1.61 

6.96 

1.76 

0.66 

1 Paso/La Tuna 

October 

0.70 

2.60 

0.84 

0.73 

IcCamey 

Annual 

12.75 

42.72 

14.28 

0.94 

Ic Carney 

January 

0.64 

0.32 

0.32 

0.98 

|r, Carney 

April 

0.77 

3.24 

0.88 

0.61 

ic Carney 

July 

1.64 

7.48 

1.60 ; 

i 

0.99 

Me Carney 

October 

1.29 

3.04 

1.20 

1.14 

leinple 

Annual 

33.37 

41.2(5 

33.68 

3.03 

ler ? 

January 

2.35 

0.44 

0.44 

4.03 

leniple 

Apri 1 

3.67 

2.92 

2.92 

3.50 

leinple 

July 

1.96 

7.40 

5.20 

2.03 

temple 

October 

2.73 

2.92 

2.92 

2.24 

“ort Arthur/Beaumont 

Annual 

54.77 

43.44 

42.44 

5.20 

Port Arthur/Beaumont 

January 

4.57 

• 

0.76 

0.76 

7.00 

Port Arthur/Beaumont 

Apri 1 

4.43 

3.12 

3.12 

5.33 

Port Arthur/Beaumont 

July 

5.71 

7.20 

7.08 

4.3C 

Port Arthur/Beaumont 

October 

3.19 

3.40 

3.36 

3.8E 


Table 3. 


Comparison of Normal Average Sponge, Potential Evapotranspiration , 2 

Actual Evapotranspiration^ , and Precipitation at Four Texas Stations. 


■P. 

Pa 

CU 


2 


E. and A. E. values extracted 
VII, Unit ed State s, 1954, 


QoLJi’ Centcrton, 


All values in inches. 


Data of the _Co^ti_r-''2. 


from: AvPraoe_^i£aJbj^>-' ter B a 1 anc e 
I'l. Thornthwaite Assoc Laboratory of Climatology, 

/ J. 
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Regarding differences in the absolute magnitude of available moisture from 
place to place (i.e., how much wetter is site *'x'' than site "Y"?), sponge 
also proves highly effective. For example, note that fieaumont's mean pre- 
cipitation is 6.4 times that of Presidio annually and, *'n April, Beaumont's 
average rainfall is 21 times that of Presidio. Sponge shows that the 
moisture gradient between these locations is actually much steeper ; Beau- 
mont's mean annual sponge value is 10.5 times greater than that of Presidio 
and, in April, Presidio's average sponge value of 0.15" is only 3 % of 
Beaumont (5.33"), a difference of 35. 5X. In other words. West Texas is 
nearly twice as ary -- compared with the humid eastern part of the state -- 
as precipitation averages alone would suggest. Considering that moisture 
availability is the primary limiting factor with respect to ecoregions and 
natural vegetation communities in Texas, it may be concluded that the sponge 
variable is an effective tool for analyzing climate-vegetation relationships. 

Recent Sponge Conditions in Texas: 1979 and 1980 . Mean monthly and annual 
sponge values for the 75 test stations, as well as an additional number of 
locations along the "Texas Transect" (see next section) were calculated for 
1979 and 1980 to assess sponge's responsiveness to inter-annual moisture 
variability. (Only briefly examined here, these variations, and their 
associations with satellite-measured vegetative index values, will be more 
intensively studied in a later phase of this research effort.) Refer to 
Tables 2 and 4 (Appendix; also Figures 2-5 (Appendix). 

It is evident that 1979 was wetter than normal in East Texas (e.g., 1979 
Xjp at Liberty - 5.31" compared with normal annual of 4.41"), while it 
was in central Texas (e.g., at Brady, in 1979 = 2.09"; 

XspHormal = 1.96") and near-normal in West Texes (e.g., at Balmorhea, 

1979, Xgp = 0.87"; X^pHormal - 0.81"). 
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Intrastate moisture conditions in 1980 were quite different than those of 
1979. East Texas was unusually dry in 1930: For example, Huntsville's mean 

annual sponge was 2.57", only SH of the long-term normal value. The differ- 
ence was especially pronounced in mid-summer, when this area experienced 
drought conditions (see Figure 4 Appendix). 

By contrast, 1980 was a relatively moist year in West Texas: Pecos and 

Balmorhea, for instance, had annual sponge values nearly 100% above their 
30-year normals (see Table 4, Appendix). 

Preliminary Assessment of Sponge 

Based on these early results, it may be concluded that the sponge is a 
useful new climatic variable for purposes of identifying and interpreting 
trans-regional moisture (and, therefore, ecological) gradients and strata. 

In fact, it may v/ell prove to be, on balance, the best such measure yet 
devised for practical large-area analysis. Accordingly, sponge is utilized 
in the following sections of this report as a generalized climatic index, 
one which is correlated with vegetative indices derived from NOAA meteoro- 
logical satellite imagery, as part of a gradient study of natural vegetation 
along a hypothetical east-west "Texas Transect". 



IV. THE VEGETATION GRADIENT UTILIZING NOAA SATELLITE IMAGERY 


Historically* two broadly conceived research methods have evolved to allow 
stratification and abstraction of plan communities, classification and 
gradient analysis (Kessell, 1979). Classification involves grouping samples 
together on the basis of shared characteristics into an abstract class of 
plant communities. Such a grouping of communities by any definition of 
shared characteristics is referred to as a community-type (Whittaker, 1975 b). 
The second method, gradient analysis, deals not with discontinuous classes 
but with continuity and gradient relationships. When the arrangement is 
along a predetermined environmental gradient, i.e. moisture, the method is 
termed direct gradient analysis . Indirect gradient analysis is the arrange- 
ment of samples along abstract axes that may or may not correspond to envir- 
onmental gradients. The process of arranging samples along one or more 
environmental gradients is called ordination (Goodall, 1954 cited in Kessell, 
1979). Since vegetation varies continuously along a moisture gradient, samples 
can indeed be ordinated. 

Frequently, the development of a useful classification system requires the 
use of ordination methods. Discontinuities in the natural vegetation are 
sought for the purpose of determining the boundaries of the comTiunity types 
recognized. These are often best determined objectively by employing the 
methods of gradient analysis and ordination. The development of a Montana 
habitat-type system {Phister, et al . T9/7 cited in Kessell, 1979) is a good 
example of the successful use of ordination in developing a classification 
system. 

Gradient modeling has been the first extensive application of gradient 
analysis to the needs of resource management information systems (Kessell, 
1979). Gradient modeling involves the linkage of a multidimensional gradient 
analysis with a remote site-specific inventory and appropriate computer soft- 
ware. Once the gradient model is complete, it can provide quantitative com- 
munity inferences (i.e., biomass, cover) if the location of each site within 
the gradient matrix is known (geographic coordinates, elevation, aspect, etc.). 


The initial step is to obtain information about the vegetation. Data on 
the vegetation can be obtained by field samples ("ground truth" studies) 
and remote methods (aerial photography and satellite imagery), Most systems 
use both. Detailed ground truth data are used to derive commun ity- types , 
whereas aerial photographs and imagery are generally used to infer the 
vegetation of unsampled areas. 

To date, a considerable number of vegetation investigations have been carried 
out using Landsat MSS imagery but very little has been attempted with the 
meteorological satellite systems, particularly the NOAA/AVHRR. Gray and 
McCrary (19S0, cfted in Gray and McCrary, 1981) obtained a high correlation 
for detection of vegetation greenness between the NOAA-6 AVHRR Large Area 
Coverage (LAC) data sets and Landsat MSS data within identical taraet areas. 

This finding led Gray and McCrarv (1981) to suggest the NOAA satellite systems 
should be used for monitoring global vegetation. One major advantage of NOAA 
over Landsat is the tremendous increase in frequency of data collection. Gray 
and McCrary (1981) anticipated that variations in the AVHRP responses will 
provide information about reactions of vegetation to moisture availability 
and thermal effects. They have demonstrated this for croplands in southern 
Texas before and after Hurricane Allen in April 1980. Since vegetation, 
particularly in regions arid and semi-arid is very responsive to moisture 
patterns, it is well worthwhile to investigate temporal changes in natural 
vegetation and how closely these relate to shifts in the AVHRP vegetation 
index. 

Quite likely, we shall ultimately discover that the success of stratifying 
different vegetation types from AVHRR vegetation indices will depend not on 
spatial distinctions but on temporal distinctions, i.e., the rate and mag- 
nitude of the spectral shift during a single season. Eventually, vegetation 
indices can be ordinated (indirect gradient analysis) and correlated to 
ground truth vegetation and climatological gradients (direct gradient analysis). 
The ultimate gradient model, incorporating both field and satellite data, may 
permit vegetation classification and monitoring of changes with minimal ground 
truthing. 


Data Acquisition and Processing 

The vegetation gradient model initially necessitated establishing a sample 
series along an environmental gradient. At NASA/JSC, we were geographically 
sitting at the eastern edge of perhaps one of the best natural east-west 
gradients in North America: Trans-Texas, along approximately 30 (®N) 
latitude. As one moves from Beaumont to El Paso, Texas, one passes through 
four major natural vegetation regions: (a) mixed pine-hardwood forest, 

(b) savannah, (c) shrubland, and (d) desert/desert scrub (r^ap 7). Pre- 
cipitation exhibits a continuum that has an annual mean > 50" {east Texas) 
toc8" (west Texas). The elevation is 0.0' at the Gulf, 5000' just east 
of El Paso. It would-be difficult to identify a better east-west continuum 
anywhere that changes gradually, yet dramatically and without any obvious 
disjunctures over a distance of approximately 750 miles. While the soil 
and geology definitely change across Texas, we do not intend to Include a 
discussion of those variables at this time. 

Our main objective has been to design a model that may ultimately allow 
vegetation classification on a global scale utilizing satellite imagery. 

IVe initially expected to use Landsat data. By a stroke of good fortune, 
we discovered that NOAA/AVHRR "Hetsat" data was not only being archived 
locally by NOAA personnel (T. Gray, D. McCrary) in a readily useable form 
but it fit our specifications perfectly. The appropriate software had been 
written by Lockheed, Inc., to be able to retrieve raw pixel data for AVHRR - 
channels 1 and 2 along specific scanlines or hands of scarlines across the 
entire state of Texas. In addition, the software provided geographic 
coordinates for each pixel. In order to obtain a specific scanline, it 
was really only necessary to provide the specific coordinates at the be- 
ginning and end of our trans-Texas transect. Since the NCAA - n series 
of satellites orbit is near-polar, sun-syncromous , and twice daily, it 
crosses a given longitude at an angle and at varying places. Consequently, 
since scanlines are perpendicular to the orbit, it was impossible to select 
scanlines that remained "isolatitude" or were exactly superimposed from one 
date to the next. It is also important to remember that NOAA scenes cover 
such an expanse that related angles to each pixel vary greatly. To permit 
comparisons, the pixel radiance values have been normalized to an overhead 


ORIGINAL PAGE 13 
OF POOR QUALITY 


Map 7. 


12 SAMPLE LOCATIONS ON TEXAS TRANSECT 


12 11 D 




LEGEND 

A - PINE/HARDWOOD FOREST 
B _ SAVANNAH 
C - SHRUBLAND 
D - DESERT/DESERT SCRUB 








It has been our intention to acquire scanlines for four cloud-free days 
during 1980 - one from each season. At this point we have only been able 
to process 3 dates (April 19, July 10, October 9); winter has been excluded 
for lack of data. We requested and received bands of 5 adjacent scanlines, 
extending essentially from El Paso to Beaumont (Figure 6). At predetermined 
locations which corresponded to our ground truth sites along the strip, we 
sampled a 25 pixel grid (5 x 5), obtaining an average grid value of pixel 
counts for each two channels. The selection of the 25 pixel sample grids 
was somewhat difficult because it was not possible to accurately ground- 
truth the transect line. The intention has been to select 12 sites, approx- 
imately 3 sites in each of the four major vegetation regions bisected (Map 7). 
Using Texas vegetation- type maps (Texas Parks and Wildlife Department - based 
on Lendsat data), original Landsat MSS scenes, aerial photos, Aeronautical 
Navigation maps (1:1,000,090), and selected vegetation references (Gould, 

1975; see Smeins, 1978), an effort was made to choose "homogeneous" natural 
vegetation sites, devoid of water, urbanization, and cultivation. The site 
locations were shifted slightly between sampling dates because the scanlines 
could not be superimposed. 

It was difficult deciding just how to initially treat the satellite data. 

Gray and McCrary (1981) have devised their own vegetation index, that they 
now rather appropriately call the Gray-McCrary Index (GMI). The GMI is 
simply the difference between the solar-zenith corrected albedo value for 
the two channels. At least initially, we are using the Landsat-derived 
normalized difference (NO) equation of Rouse, et al. (1973) and Heering, 
et al . (1975) where: 

Channel 2 - Channel 1 


Channel 2 + Channel 1 
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FIGURE 6. 


Normalized difference values were obtained for each of the twelve (1?) 
sites on the three (3) dates. The next step was to decide how to best 
statistically treat the NO values. At this stage we have simply been 
able to do some preliminary analysis of variance. Normalized difference 
values were plotted as a function of longitude, sponge and certain vegetation 
characteristics {i.e., biomass, net productivity, leaf area). 

Vegetation Regions on the Texas Transect 

The Texas transect selected for our model essentially runs from Beaumont 
(94®W) to El Paso (106°) and the 12 sample sites have been numbered east 
to west (see Map 7). Since it was not feasible to visit the transect for 
optimal site selection, it v/as necessary to utilize vegetation maps. At 
this level of our investigation, the ground-truth precision was not terribly 
critical since our initial concern has been to get a general feeling for the 
potential of NOAA imagery for global vegetation stratification. 

Gould's (1975) vegetation map and discussion of the vegetation regions of 
the state is large-scale but is the best complete map available. The 
Texas Parks and Wildlife Department is in the process of completing a state- 
wide series of land-use classification maps based on Landsat data. Utilizing 
primarily those two sources, it appears as though ou?' transect bisects four 
major vegetation regions. Table 5 enumerates those four regions from east 
to west, their approximate longitudinal boundaries on the transect, and 
document vegetation. For a more complete vegetation description, see Gould 
(1975), Texas Parks and Wildlife vegetation-type maps, and Smeins (1978). 
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VEGETATION REGION 

GOULD (1975) 
NOMENCLATURE 

L0NGITU0E(°W) 

LIMITS 

DOMINANT VEGETATION 

Pine-Hardwood 

Forest 

Pineywoods 

J3.45°-95.45° 

loblolly pine, slash 
pine, sweetgun, oak, 
elm, pecan, blackgum 

Savannah 

* 

Post oak savannah, 
Blackland prairies , 

95.45°-98.30° 

« 

grasses, oak, elm, 
hackberry, (cropland) 

Shrublan^ 

Edwards Plateau 

98.30°- 102° 

oaks, ash, juniper, 
mesquite, (rangeland) 

Oesert/Desert- 

scrub 

Trans Pecos 
Mountains & basins 

102°-106° 

creosote, tarbush, 
yucca, (rangeland) 


TABLE 5. The vegetation regions, their approximate longitudinal boundaries, and 
dominant vegetation along the Texas transect. 
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Results and Discussion 
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1 . Nonna! i zed difference (ND) as a function of longitude . As pre- 
viously mentioned it was impossible to insure superimpcsi tion of pixel 
sampling locations between dates because the satellite orbit fluctuates. 

For example, while still in pine-hardwood forest, site 1 on April 19 is 
not geographically identical to site 1 on October 9. The tabular and 
graphic summaries when the three sample dates were individually plotted 
as a function of longitude are presented below (Table 6, Figures 7, 8). 

There is generally a high correlation between MD and longitude 
(mean r = .756); MD decreases from east to west. 

The sample size of 12 is not sufficiently large enough to merit serious 
discussion as to significant differences among the three sampling periods. 
Neither can we eliminate the real possibility of cloud cover affecting 
reflectance values. Ue definitely know that there was some cloud cover 
in West Texas on October 9. What the regressions do suggest is what one 
would expect knowing the phenological nature of the vegetation regions on 
th«i transect. April and July values are high in the mixed forest because 
the deciduous trees have leafed out. By October, they have dropped their 
leaves, reducing their "greenness". Progressing westward across the state, 
■? grea'er percentage of the perennials are non-deciduous but the vegetation 
becomes less dense and more dependent' on infrequent precipitation. In 
addition, the amount of exposed ground in the desert scrub poses problems 
of separating soil spectra from vegetation spectrore, lectance (Miller, 

Lee 0., pers. comm). 

The variation between seasons appears to diminish going from east to west. 
The mean regression for all three dates (Figure 8) indicates not only that 
the greater between-date variation is in the pine-hardwood forest but also 
that April, and to a lesser extent July are well above the line for the 
three dates. We would have expected much higher values during the annual 
blooms in April in the Trans-Pecos but more extensive sampling may clarify 
this. 
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Site 


April 19 


July 10 


October 9 


Combined 


1 . 

.258 

.238 

.161 

2. 

.289 

.229 

.171 

3. 

.286 

.231 

.215 

4. 

.259 

.136 

.147 

5. 

.224 

.205 

.151 

6. 

.126 

.137 

.090 

7. 

.138 

.125 

.033 

8. 

’.062 

.141 

.101 

9. 

.038 

.159 

.112 

10. 

.044 

.029 

.051 

11. 

.046 

.005 

.044 

12. 

.035 

.001 

.063 

X 

.150 

.128 

.112 

Y 

2.526 - .024X 

2.258 - .022X 

1.23 - 

r' 

.752 

.934 

.636 


.130 

2.00 - .019X 
.750 


table 6. Indidivudal normalized difference (ND) values, means, regression 
equations, and r^ for 3 dates in 1980 along TexoS transect. 
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VEGETATION INDEX REGRESSION LINES FOR 3 DAYS IN 1980 

TEXAS TRANSECT 
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CES AND MEAN REGRESSION 






2. Sponge index as a function of longitude . For a discussion of a 
new moisture variable, the sponge, we refer you to earlier sections II, III 
of this report. Long-term annual "normals” (1941-70) were calculated for 
26 stations along or near the Texas transect. These values were plotted 

as a function of longitude (Figure 9). The results indicate an extremely 
high correlation (r = .911) between sponge and longitude; sponge increases 
(as does precipitation) from west to east. 

3. Normalized difference (ND) as function of sponge . Since both the 
vegetation index,* ND, and sponge showed similar positive correlation with 
neographic position on the transect, it seemed appropriate to interpolate 
sponge values from the 26 stations (see Table ^ in Appendix; *^aps 8,9) 
along the transect "‘ir each of 36 ND values (3 values for each of the 12 
sites). When ND v/as plotted against sponae, there was a good correlation 
(r = .777) (Figure 10). Regressing the two regressions, (NO (Figure 8) 

9 

VS. Sponge (Figure 9 ), the result is a very please (r" - 1.001) but transect- 

limited, prediction model that permits estimating longitude, sponge index, 

and ND value, requiring input of only one of the 3 variables (Figure 11). 

This model serves to illustrate the very high correlation between the 

vegetation index and the sponge. To further establish this correlation, 

if one plots the highest ND value of the three dates at each 12 stations 

2 

against the long-term sponge value for that particular month, r = .946. 

Normalized difference (ND) as a function of vegetation . As earlier 
mentioned, it was not feasible to actually ground truth the transect to 
verify and quantify the vegetation. No doubt this should be done at some 
later stage. In lieu of a better alternative, biomass, net productivity , 
and leaf area estimates from Whittaker and Likens (1 973 ) wer*" used (see 
Table 7}. Since the predominant vegetation limiting factor along the Texas 
gradient is moisture, these vegetation parameters predictably decline from 
east to west. It has been previously demonstrated that all of these para- 
meters have been correlated with spectral data (see introduction to Tucker, 
et al . 1981 for literature review). While general, when ND means for each 
site are plotted against biomass, net productivity, and leaf area mean 
values, the results are prophetic as to which vegetation characteristic has 
' the highest correlation with the satellite data (Table 8). 
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NET PRODUCTIVITY 

BIOMASS 

LAI 


qm/M^/Yr 

Kg/M^ 

M^/M 

Temperate Evergreen/Deciduous 
Forest 

1250 

32.5 

8.5 

Savannah 

900 

4.0 

4.0 

Shrub! and 

700 

6.0 

4.0 

Desert/Desert 

90 

0.7 

1.0 


TABLE 7. Moan values for world-wide estimates of not productivity biomass 
and leaf area index (Whittaker and Likens, 1973). 
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2 

Net primary productivity (g/m yr) 

.895 

Leaf area (M^/M^) 

.815 

Biomass (Kg/M^) 

.650 

TABLE 3. Coefficients of determination (r ) when mean 
.'Texas transect sites are plotted against net 
area, and biomass. 

NO values for the 12 
productivity, leaf 






33 


SPONGE INDEX (INCHES) 


LONG TERM (1941-70) SPONGE NOf JIALS FROM 26 STATIONS - 
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VEGETATION-SPONGE PREDICTION MODEL 
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STATION LOCATIONS FOR SPONGE CALCULATIONS 
ALONG TEXAS TRANSECT 

(PERIOD 1941-70 AND 1980, EXCEPT AS' NOTED IN TABLE 2) 
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The Inference is that net productivity has the highest correlation to the 
vegetation index, biomass the lowest. Several studies have recently shown 
that currently used remote sensing techniques are not sensitive to non-green- 
leaf components of the phytomass. However, there appears to be a high corre- 
lation of spectral data with green-leaf area (biomass) and net production of 
certain vegetation types (see Introduction, Tucker, et al . 1981). Our selection 
of NO as a vegetation index was largely founded on Deering and Haas' (1977) high 
correlation between Landsat-derived ND and rangeland biomass. 

In Table 7, it is* noted that of the three parameters of vegetation, biomass 
is the only one that doesn't consistently decline from east to west along the 
transect. Shrubland, consisting largely of woody perennials, does not produce 
the annual net production that a savannah, containing more herbaceous annuals 
would, but its accumulative biomass would be greater. The ND does not re- 
spond to the increase ir biomass from savannah to shrubland because much of 
that biomass is tied up in non-green components in snrubland which is not as 
true in the savannah. Because net productivity and leaf area more closely re- 
flect the actual spectral component of those vegetation regions on our transect, 
we would anticipate their previously higher correlation with ND. 

5. Vegetation-Sponge Index (V$I) . Having previously established a high 
correlation between ND and sponge, we would like to propose a new index that 
represents the multiplicative of the two variables: the Vegetation - Sponge 

Index (VSI). 


ND X SPONGE = VSI 

The mean VSI values for the 12 sample sites are presented below in tabular 
and graphic form (Table 9, Figure 1?). 


SITE 

1. 

VFRETATION REGIONS 
Pine-hardwood forest 

SPONGE 

NO 

VEGETATION-SPONGE 
INDEX (VIS) 

5.20 

.2192 

1.13 

2. 

II 

4.61 

.2298 

1.06 

3. 

M 

4.40 

.2440 

1.07 

4. 

Savannah/cropland 
■ •* - 

4.10 

.1806 

ft 

0.74 

5. 

M 

3.31 

.1935 

0.64 

6. 

Shrubland 

2.70 

.1176 

0.32 

7. 

II 

2.12 

.0989 ■ 

0.21 

8. 

11 

1.71 

.1012 

0.17 

9. 

II 

1.36 

.0695 

0.10 

10. 

Desert/desert scrub 

0.70 

.0414 

0.03 

n. 

H 

0.67 

.0315 

0.02 

19 

II 

0.63 

.0333 

0.02 


TABLE 9. Sponge, normalized difference (NO), and Vegetation-Sponge Index 
(VSI) mean values for three dates in 1980 (April 19, Ju y , 
October 9) for 12 sites along Texas transect. 
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Figure 12. 






When VSI Is plotted as a function of longitude, the separation of the four 
major vegetation regions becomes much more apparent (Figure 12). It is 
ever possible to suggest boundaries for the four classes: 

>1.0 pine-hardwood forest 
0. 4-1.0 savannah/cropland 

.05-0.4 shrubland 

< .05 desert/desert scrub 

Using these suggested limits, it becomes readily apparent that the least 
within-class variation occurs in desert and forest; the greatest variation 
is in the shrubland and savannah/cropland classes. This seems consistent 
with what we would expect. The desert and mixed forest would be more homo- 
geneous in the sense that desert has extensive bare soil during much of the 
year while the forest would have relatively little. The shrubland and 
savannah would be considerably more het rogeneous. 
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V. THE TEXAS MODEL: CONCLUSIONS AfiD PROJECTIONS 


It has been demonstrated that the sponge variable is a superior tool for 
the analysis of climate/vegetation relationships. Furthermore, NOAA/AVHRR 
satellite data proved useful for vegetation stratification. Finally, a 
preliminary multivariate model of vegetation distribution, the VSI, was 
developed based on an experimental east-west Texas gradient. The next 
stage of this research effort will Involve a more extensive analysis of 
the application of the model to the Texas transect (e.g., increasing the 
sample size), to be followed later by refinement of the model v/hich will 
then be tested against other natural vegetation regions across North 
America. It is anticipated that ultimately such a model may be utilized 
for global vegeta'.ion surveys, and as a means of remotely monitoring 
vegetation region dynamics (e.g., desertification). 
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1.38 
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1.46 
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3.98 
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1.97 
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3.70 
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qUfcUTf. 

/T.e’; tyr '» 




















































































































































































JOJJDQUX ERAME 


ORIGINAL PAGE » 
OF POOR 


Table 2(C). Mean *, x\s Precipitation and 8- Inch Sponge Homa 
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Itude 
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m: 




™nr 


SL 


aPrIL 


IS 2 . 
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T 
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"luo 
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* 

Athens" " 

32.20 

95.85 

2.51 

4.30 



3.30 

5.10 

2.58 

4.78 

4.64 

4.51 

5.17 

4.81 

3.14 

3.77 

1.41 

I. 

t 

Eoaumoi ; 

1 

3C.08 

94.10 

4.57 

7.00 

4.38 

7.47 

2.82 

6.41 

4.43 

5.33 

4.79 

5.10 

4.75 

4.33 

5.71 

t. 

I'.r'TT’ha'n ‘ ^ 

30.15 

96.40 

2.78 

4.84 

3.12 

5.06 

2.46 

4.53 

3.71 

3.97 

4.43 

4.01 

3.24 

3.29 

1.90 

?. 

r ^16 

Cuero 

29.08 

97.25 

2.03 

2.08 

2.41 

3.18 

1.41 

2.71 

2.95 

2.41 

3.79 

2.91 

3.18 

2.64 

2.22 

1. 

t 

F re^?por*. 

23.93 

95.33 

3.41 

6.03 

3.48 

5.31 

2.23 

5.55 

3.25 

4.68 

3.29 

4.36 

4.17 

3.99 

5.46 


;Hal lett: • 
j ville^^' 

29.45 

96.93 

2.25 

3.36 

2.70 

3.58 

2.08 

3.29 

3.01 

2.86 

4.65 

3.39 

3.53 

3.18 

2.22 

1 . 

I 

I'.ntsvii (' 

30.72 

95.57 

3.80 

6.25 

3.76 

6.65 

2.98 

5.85 

4.34 

5.13 

4.35 

4.70 

4.57 

3.97 

3.28 

3. 

i { 

iU berty^’ 
1 

30.05 

94.32 

4.21 

6.21 

4.18 

6.58 

2.77 

5.67 

4.16 

4.73 

4.50 

4.51 

4.33 

3.82 

4.65 

3. 

jLufki 

31.23 

94.75 

4.29 

6.32 

3.83 

6.70 

3.49 

5.96 

4.33 

5.16 

4.81 

4.73 

3.48 

3.75 

2.83 

2. 


32.53 

94.35 

4.40 

7.C5 

3.85 

7.62 

4.03 

7.14 

5.02 

5.42 

4.78 

5.52 

3 ; 45 

3.94 

2.97 

2. 

1 

h'3 1 ac^O!. 

1 

28.63 

96.613 

2.10 

3. 70 

' 

2.96 

4.10 

1.61 

3.69 

2.80 

3.29 

3.21 

3.23 

5.34 

3.90 

1.85 

3.( 

. 15 

.■• an s 

33.67 

95. 57 

2.59 

5.49 

3.52 

6.04 

3.49 

6.12 

5.59 

6.19 

5.42 

6.27 

4.46 

5.01 

3.76 

3.( 

Ceg-ji n^ 
- 

2.-. 57 

97.95 

1.72 

2.67 

2.34 

3.00 

1.75 

2.71 

3.27 

2.67 

3.18 

2.87 

3.17 

2.53 

1.89 

.{ 

L . n 

i ex.' r'Kaf ' 

33.30 

94.17 

4.09 

6.51 

4.05 

7.43 

4.04 

7.33 

3.91 

7.15 

4.81 

6.24 

3.54 

4.35 

3.16 : 

LC 

1 ii> ;n 3Cf. 1 ! 

32.38 

96.85 

2.19 

4.04 

] 

2.98 

4.59 

2.19 

4.34 

1.51 

4 . 09 j 

4.87 

4.57 

3.60 

3.87 

?.21 2 

’.4 


cecc noted. All values in inches. Data from Cl i'^iatolo oical, Sutn -arigs_ for Texas, Nat 
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3.24 
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3.17 
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4.65 

3.39 

3.53 

3.18 
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3.05 
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2.57 

3.10 

2.93 

2.62 
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1.98 
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3.80 

1.97 

3.35 

2.73 



1.73 

2.68 

3C.:.3 

. 49 

15 

4.81 

6.24 

3.54 

4.35 

3.16 

3.03 

2.71 

2. 16 
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1947 - 1966 


^1931 - 1969 
^931 - 1970 
^1932 - 1969 
^933 - 1967 
^1933 - 1968 
®19:" - 1969 
^937 1966 

®1937 - 1969 
^1938 - 1967 
’^1939 - 1968 
^^1939 - 1969 
^^1940 - 1966 
C ^^1940 - 1969 
^^1941 - 1967 
^^1941 - 1968 
^^1942 - 1971 
^^943 - 1970 
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’•^947 - 1970 
^°1948 - 1967 
^^1948 - 1969 
^^1948 - 1970 
^^1949 - 1966 
^^1949 - 1967 
^^1949 - 1968 
^^1949 - 1969 
^^1949 - 1970 
^®1950 - 1969 
^^1950 - 1971 
^°1951 - 1966 
^^1951 - 1967 
^^1953 - 1970 
^^1956 - 1970 
^^1959 - 1971 


Periods of Record Other Than 1941-1970: Tables 2(a)- 2(e). 
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I. DEFINITION OF A FUEL CELL 

A fuel cell is a device for converting chemical energy directly into 
electrical energy where the reacting materials are supplied from 
reservoirs external to the electrochemical cell. 

II. HISTORY OF THE FUEL CELL^^^ 

In 1839, Sir William Grove of England constructed the first known 
fuel cell. In 1842 he built a "gaseous-voltaic battery" (Figure 1) 
with bO H 2 - O 2 cells in series. Fifty years later the name "fuel-cell*' 
was given to the bank of cells built by Mond and Langer which operated 
with a 50% efficicjicy and produced Ih watts. However, it wasn't until 
the early 1930 's that Sir Frances Bacon in England built a 5 kilowatt 
fuel cell using nickel electrodes of dual porosity, H 2 and O 2 , and KOH 
electrolyte. This cell operated at 205°C and 600 psi. Industry (e.g., 
Allis Chalmers) began developing fuel cells in the 1950' s for terrestrial 
use and then NASA gave impetus to the technology in the 1960's with funding 
to companies such as GE and Pratt and Whitney who developed fuel cells 
for the Gemini and Apollo programs respectively. The advantages of 
fuel cells were seen to be high efficiency compared to heat engines, 

r ' * 

operation independent of environmental conditions, and potable water as 
a by-product for use by crew (for hydrogen-oxygen fuel cell). 

III. THEORY OF THE FUEL CELL 

Since the fuel cell is a voltaic cell in concept then there is an 
oxidizable material (fuel) and an oxidizing material (oxidant). Some 
fuels that have been used include hydrogen, hydrazine, methane and 
methanol; oxidizers have coninonly been oxygen or air.^^^ The main fuel 
cell developnient thus far has \ aen on the hydrogen-oxygen system and thus 



FIGURE 1 FUEL CELLS AS i'»OWN BY GROVE 



FUEL CELLS 
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the remainder of this report will deal with this system utilizing an 
alkaline (KOH) electrolyte. In Figure 2, a schematic of the fuel cell Is 
shown. In contrast to a conventional battery, note that in a fuel 
cell, the electrode materials do not react chemically but they provide 
reaction sites for the fuel, oxidizer and electrolyte as well as a 
conduction path for electrons through an external circuit. Between the 
electrodes in the cell is the electrolyte which is shown as an alkaline 
solution such as potassium hydroxide in water. The electrolyte could be 
an acidic solution also with appropriate cell materials. Through the 
electrolyte the electrical circuit is completed by. ion transoort. 
between the electrodes. One advantage of the fuel cell is obvious from 
the diagram since the use of several external storage tanks or a re- 
chargeable tank system allows for continuous operation of the cell 
without interrupting the external electrical circuit. By contrast, the 
circuit would have to be interrupted to recharge secondary batteries 
or to replace primary batteries. 

In order to understand the fuel cell operation consider the basic reaction: 
^ 2 ( 9 ) + hOig) -*■ H20(i 1 + Energy 

Figure 3 illustrates how the energy released is utilized in a heat 
engine-generator system and in a fuel-cell. The heat engine must operate 
at high temperatures to get reasonable efficiencies but, then other 
components still bring down the overall efficiency to near 30%. The 
fuel cell has a much higher overall efficiency at ordinary temperatures. 

For example the fuel cells on STS-1 were operating at about 64% efficiency 
based on data given in Rockwell Internal letter (PG-MWF-81-027). 
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FIGURE 2 SCHEMATIC OF HYOROGEN-OXYGEN FUEL CELL 


>rELECTRICAL 

WORK 



(ANODE) (CATHODE) 


Anode Reaction 
Hg + 20H” EHgO +’ 2e 


Cathode reaction 
Jsi02 + ” 2 ® ^ 2e -► 20H“ 


OVERALL REACTION 


^2 ^ ^2 ^2® * work + Heat 


FIGURE 3 BLOCK DIAGRAM OF HEAT ENGINE-GENERATOR VS. FUEL CELL SYSTEM 


Heat Engine-Generator 



Fuel Cell 


Hp + %02 H 2 O 


^Energy 

Cell 


Elefctrical 

Energy 

60t 


Nr 

Heat 

Loss 

Thermodynamic Heat Loss, TaS ■ lOX 

Current and Voltage 
Inefficiency loss * 30X 

Total Heat Loss * 40X 


* Carnot efficiency for heat engine operating Approximate Overall Efficienc y 60K 

between 900°F and 70'^F. 


Approximate overall efficiency 0.90 x 0.60 x 0.80 x 0.95 * 0.41 or__41% 





3 


4 


« 


The energy given off when hydrogen and oxygen are burned to form liquid 
water at constant pressure Is called enthalpy, aH, of this reaction and 
for Ig-mole (18g) of the liquid water formed, aH * -68.317 kcal. The 
negative sign means that the reaction Is exothermic (heat Is liberated). 
That would be the maximum heat available to operate a heat engine. The 
efficiency of that engine would be t^« basis of 

one g-mole H^O liquid formed). 

The maximum efficiency of a heat engine is limited by the Carnot cycle 
and Is given by 


nc 


T> - T, 

To 


X 100 where T 2 Is Absolute Temperature 

of high T«nperature Reservoir. 
Tx Is Absolute Temperature of 
low Temperature Reservoir, 


In Figure 3, t 2 * 900 F 
T 2 « 1359°R 
ti = 70°F 
Ti = 529°R 

Therefore - (^ . 359 , . - . ^52 9) ^ ^ , 

Other energy losses due to boiler Inefficiency, friction and electrical 
loss reduce the efficiency of heat engine-generator system to about 40% 
and In actual practice it would be less than this. 


For the fuel cell operating at a fixed temperature and pressure, the 
energy available for conversion to electricity Is not equal to AH 
but Is given by tK Gibbs free energy aG, where 
AG * aH - TaS 

and AS for H 2 (g) + H02(g)-^ H20(l) Is -0.039 ~- 

mole K 

Therefore at 25°C 

AS . . 17 f 298C0.0390) - -56.69 
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a6 Is the energy available for external work in a system other than P-V 
work and the only work we are considering ts electrical work so 
a 6 ■ -nFE 

where E is voltage of fuel cell, F is Faraday's constant and is equal to 
23.060 and n equals the number of moles of electrons involved in the 
reaction. To* find n, the reaction can be broken down into an oxidation 
half and a reduction half where by definition oxidation occurs at the 
anode of a cell and reduction occurs at the cathode of a cell. The 
overall reaction can be divided as follows: 

Oxidation (Anode) H 2 (g) + 20H *> 2H20(‘*1 + 2e 

Reduction (Cathode) %02(g) + H20(i) + 2e 20H 

Overall Reaction H 2 + hOz H 2 O 

From the half reactions n=2. since 2 moles of electrons are transferred 
through an external electrical load for each riiole of water produced. 

Since a 6 * -nFE 


then E * -56.69 kcal 

mole 
-(2 moles) ( 


-) = 1.229Y. 


This is the maximum voltage obtainable from H 2 - O 2 fuel cell in basic 
solution at 25°C. However, to compare on an equivalent basis to the 

heat engine the theoretical efficiency of a fuel cell should be given by 

_ AG X 100 _ -56.69 x 100 00 ne/ 

H T iR 

or in terms of the voltage ratio, it would be 

It '55:^ Tjm 


This is the maximum thermal efficiency of the fuel cell at 25°C and 1 


atntosphere pressure. 
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The actual thermal efficiency of the fuel cell Is 

“ 67.52E(%) where E U the measured voltage 

s 

output of the cell. If the electrical efficiency of the fuel cell is of 
interest then ^ ** 81.37E(%) where E is the actual 

voltage and 1.229 is the' maximum theoretical voltage of the cell at 25°C 
and 1 atmosphere. 

The electrical efficiency should contain a factor for current inefficiencies, 
however current efficiencies are normally > 95^ unless there are mechanical 
problems v/ith the cell such as cell shorting, diffusion of reactmts to 
opposite electrode to react directly, or leakage of reactants out of the 
cell. The voltage inefficiency has been a major concern of developmental 
work in the fuel cell. The factors involved that result in the cell 
voltage being less than the theoretical voltage of 1.229V can be summarized 
into three categories: 

(1) Electrode Activation 

(2) Ohmic Losses 

(3) Diffusion Losses 

The first effect of electrode activation can be understood from the fact 
that almost every reaction has an activation energy as shown in Figure 4. 

Even though the products are more stable (lower energy) than reactants, 
there is a certain energy barrier that the reactants must overcome in 
order to react and form products. This energy barrier is the activation 
energy. For example, the reaction of 82 ( 9 ) + ^ 02 ( 9 ) H 20 (l) is very 

exothermic (-68.3 kcal/mole HgO), and yet a mixture of H 2 (s) snd 02 ( 9 ) 
is stable almost indefinitely because of a high activation energy for this 
reaction. Heat or a spark added to this mixture would then cause an 
explosive reaction with the formation of H 20 (l). 
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FIGURE 4 ENERGY RELATIONSHIPS IN A CHEMICAL REACTION 
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Another alternative to going over the high activation energy barrier is to 
provide another reaction mechanism by means of a catalyst. A catalyst 
is a substance that alters the activation energy and thus the rate of 
the reaction but is not used «:p itself in the reaction. Most catalysts 
lower the activation energy and speed up the reaction, but there are cases 
when "negative" catalysts or inhibitors are desirable. 

Thus in Figure 4 the reaction path with a catalyst has a much lower 
activation energy than the normal path and the reaction will proceed 
faster. Increasing the temperature will increase the thermal energy of 
the reactants and will also increase the rate of the reaction since the 
reactants will have mors energy to go over the activation energy barrier. 

In the fuel cell there are two half reactions and each one must be 
catalyzed separately and each one has a separate effect on reducing the 
voltage output of the cell due to its activation. At the anode of the cell 
the following reaction takes place: 

U 2 + 20H" -> 2 H 2 O + 2e 

Many of the transition metals. Ruthenium, Indium, Platinum, Paladium, 

Osmium, Gold, Rhenium, Tungsten, and Nickel catalyze this reaction. 

For practical use Platinum-Palladium mixtures have been found to be the 

best for catalytt*: activity and long life. However, carbon monoxide and 

chloride ions have had a detrimental effect on the catalytic activity of 

Pt-Pd mixtures. The catalyzed reaction appears to be a 2-step mechanism 

as illustrated below. . . (1), physical adsorption of hydrogen on the 

metal surface followed by (2), dissociation into chemi-sorbed hydrogen 

atoms which are chemically reactive. 

H-H— ^ rH-H% H H 

Pt-Pt-Pt- -Pt-Pt-Pt-Pt- -Pt-Pt-PVPt- 

( 1 ) ( 2 ) 
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The anode reaction takes place on a catalytic electrode with very 
little voltage loss and usually It Is not a problem In the fuel cell cp- 
eratlon. At the cathode of the cell the. other half of the reaction 
Is as follows: 

HO 2 ^ 2e + H 2 O 20H’ 

This reaction is also catalyzed by many transition metals but best results 

have been achieved with a Gold-Platinum alloy catalyst. This reaction 

is more complicated than the anode reaction and has several steps In 

f3l 

the mechanism. One postulated mechanism is as follows: ' 

O 2 + Pt -> Pt(02) 

PtCOz) + e + H 2 O ptCHOz) + 0H“ 

PtiH 02 ) + e Pt(H 02 ") 

PtCHOa ) + 2e + H 2 O - Pt + 30H~ 

The cathode exhibits the greatest voltage loss due to activation. If 
the perhydroxyl Ions do not break down In the last step, the voltage 
loss is very high. Gold seems to be a very effective catalyst in the 
decomposition of perhydroxyl ions and therefore. It is used in combination 
with platinum on the cathode to reduce this activation loss. In general, 
activation losses can be reduced by better catalysts, more surface 
area of catalyst exposed for reaction. Increasing temperature of system, 
and Increasing the pressure of the reactants. 

The second type of voltage loss in fuel cells is the ohmic loss. This 
is due to the electrical resistance internal to the cell (usually called 
internal resistance of a battery). This generally depends upon the cell 
electrolyte, electrode geometry (area and distance between electrodes) 
and the temperature of the system. The use of an alkaline cell meant that 
some type of hydroxide must be used and potassium hydroxide (KOH) was 
chosen because it has a high electrical conductivity. Within limits. 
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the more concentrated solution of electrolyte also Is more conductive, but 
viscosity and phase change effects rule out very concentrated solutions 
except at higher temperatures. For cell operation at approximately 
180®F, KOH concentrations of about 30% have been satisfactory. Design 
considerations determine the electrode area for required performance. 

The last type of voltage loss Is that due to diffusion effects. Both 
electrodes Involve reactions with Ionic species. For the reaction to 
proceed It is necessary that the active catalyst be In contact with 
reacting gas and dissolved ionic species. At very high current flow 
the diffusion of Ionic species into or away from reacting sites becomes a 
current-limiting factor. 

The diffusion -limited region can be extended by Increasing temperature, 

increasing electrolyte concentration and by electrode design, but eventually 

current will become diffusion-limited. The three effects are summarized 
C4) 

in Figure 5.' ' This figure summarizes many of the concepts for the 
hydrogen- oxygen fuel cell. It gives the total heat energy possible and 
thus a scale for the thermal efficiencj^/^j and also it gives the maximum 
free energy which determines the maximum voltage of a cell and the 
asisocfated electrical efficiency based on that value, 1^^. The activation 
loss Is evident at low currents. There is a diagnostic test called the 
Taffel plot in which voltage Is plotted as a function of the log of 
current In the low current region. The slope of this plot compared with 
a norm established from other cells determines the quality of activation 
for the electrodes. The ohmic loss is a linear function of current, but 
finally the current reaches a limiting value due to diffusion of ionic 
species in the electrolyte. The linear region of the graph is the range 
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FIGURE 5 HYPOTHETICAL HYDROGEN-OXYGEN FUEL CELL 
PERFORMANCE NEAR ROOM TEMPERATURE 
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of performance utilized In the design of a fuel cell for a specific 
appi icatlon. 

IV. ORB ITER FUEL CELL 

Ktth thts background of fuel cells, the Orbtter fuel cell will be considered 
tn regard to its components, present status and future improvements. 

A. CELL COMPONENTS 

Figure 6 illustrates a single fuel cell in a schematic as well as con- 
figurational diagram.'' ‘ In actual use 32 of those cells are connected 
in series and then two of these stacks are connected in parallel to form 
the power section of the Fuel Cell Powerplant. This provides 27.5 - 32.5 
VDC with a 2 - 12 kw output. In Figure 7 this power section is shown 
with auxiliary equipment such as reactant gas regulators, water condensing 
equipment, coolant system and associated valves and pumps. This 

schematically represents the fuel cell, powerplant. However^ this report 
deals with the fuel cell itself and referring back to Figure 6, the componen 
are as follows: 

Oxygen Electrode CCathode) - Nickel screen plated with gold and 
then 20-25 mg/cm^ Gold-Platinum alloy with 20% TFE deposited on 
screen and sintered at 590°F (TFE is tetrafluorethylene and acts as 
a wet-proofing agent which prevents electrode flooding while 
providing more interfacial area between catalyst, gas and electrolyte). 

2. El ectrolyte- Potassium fiydroxide in water with an allowable concentrate 
range of 25 - 45%. 

3. Matrix ^ Porous asbestos mat to contain electrolyte. RAM refers to 
reconstituted asbestos matrix devleoped as a process by H. McBryar 
at JSC to produce a more uniform matrix with high bubble pressure 
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(prevents crossover reactions of one gas to the other electrodes). 

Also a side benefit Is that some contaminating material has been 

leached out. 

4. Hvrironen Electrode - (Anode) - Nickel screen with gold plating 

and then 10 mg/cn? of Platinum-palladium with 30% TFE is deposited on 
screen and sintered at SOO^F. 

5. Electrolyte Reservoir Plate (ERP ) - Porous Nickel plate is present to 
provide additional inventory of electrolyte and to accomodate volume 
changes of the electrolyte during use. This plate is perforated in 
order to permit passage of hydrogen to the reactive electrode surface. 

6. Separator Plates - between series cells are gold-plated magnesium 
separator plates grooved for reactant gas flow or coolant flow as 
specified by power plant design. 

7. Frame - Cells are separated with fiberglass-epoxy material and sealed 
against fluid leakage with rubber rings seated .into grooves around the 
periphery of the cell plates. 


3, STATUS'OF ORBITER FUEL CELL 

A shuttle prototype stack of 32 cells referred to as DM2A fuel cell was 

evaluated by the Thermochemical Test Branch at JSC and the following comments 

( 6 ] 

were made concerning cell performance:" 

1. The fuel cell experienced .a high performance degradation rate for the 
first 100 hours of operation; the degradation rate then tapered to a 
normal value. 

2. The Phase I performance loss was attributed to a change in catalyst 
activity. 

3. Approximately half the performance loss experienced during Phase I 
testing was regained upon startup for Phase II testing. 
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4. The performance recovery was attributed to regeneration of the degraded 
catalyst activity. 

5. Fuel cell performance during Phase II testing was erratic and unpredlctabl 

6. The erratic performance complicated the definition of the fuel cell 
degradation rate; however, the average degradation rate during Phase II 
testing was approximately 0.168 millivolts per hour. 

7. If the cause of the high initial performance degradation rate could 
be determined and alleviated, the fuel cell stack useful lifetime could 
be extended considerably. 

8. High continuous loads (150 amperes) resulted in accelerated degradation 
as much as 16.2 times that of moderate current loads. 

9. Less than half the degradation associated with high current loads was 
due to a change in catalyst effectiveness; the majority of the loss was 
due to increased internal resistance or increased electrode diffusion loss. 

10. The fuel cell required as much as two hours to achieve voltage stability 
following a step load change. 

11. The stack and individual cell voltages exhibited high amplitude damped 
oscillations with "instantaneous" load increases followed by typical 
exponential decay to the voltage value commensurate. with the new load. 

12. Purge gas flow blockage in specific cells precluded normal boot- 
strap startups following inerting after a total on load operating time 

of 3818 hours. * 

13. Oxygen dual feed mode operation after an ET of 2954 hours considerably 
extended the allowable oxygen purge interval compared to single feed mode; 
the single feed mode resulted in rapid inert buildup in select cells. 

14. Stack voltage performance degradation between on-load elpased times 
of 1920 hours and 5150 hours primarily resulted from increased electrode 
diffusion losses, with some catalyst degradation also apparent. 
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Previous work reported cell performance degradation due to anode flooding 
which seemed to be caused by a gray deposit containing silicon on the anode. ^ 
Other work had shown that CO or COg contamination of reactant gases 
result In the formation of carbonates In the electrolyte replacing the 
hydroxide lons.^®^ This Is detrimental to cell perfomance since the 
resistance of the cell Increases due to a decrease In ionic conductivity 
as hydroxide Ions are replaced by carbonate Ions. The decrease In cell 
perforamnce was evident in electrode and elect'*olyte degradation. 

However, the source of degradation was found to be contaminants In reactant 
gases, matrix degradation and frame corrosion. 

C. POTENTIAL IMPROVEMENTS IN THE ORBITER FUEL CELL 

Goals for the orfalter fuel cell enhancement program recommended by Power 
Systems Division of United Technologies Corporation are to decrease 
weight, increase operating life to 10,000 hours. Increase power output 
from 12 kW to 20 kW and reduce cost. Proposed areas of work to meet 
these goals are as follows: 

1. Fleecing electrode process 

2. , Superior Matrix 

3. New Frame Material 

4. Graphite Electrolyte Reservoir Plate 

5. Platinum-on-carbon anode 

6. Bonded cell construction 

D. BACKGROUND INFORMATION TO THESE PROPOSED CHANGES 
1. Fleecing Electrode Process 

The standard process for electrode manufacture Involves filtering a 
suspension of Tgfion and catalyst particles onto a teflon sheet, trans- 
ferring to the electrode screen, wash catalyst mixutre, dry, and sinter. 
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This process has produced electrodes with undesirable variations In 
performance. 

An improved process Involves the addition of a proprietary fleecing agent 
to the suspension before fllterlng.'^Hhese flocced electrodes when 
compared to standard electrodes have a finer and more uniform grain 
structure. This Is the result of a more unlfom distribution of Teflon 
and catalyst particles. These flocced electrodes have less variation In 
production runs than the standard electrode reducing the rejection rate 
of electrodes from 10% to almost zero. 

Test data on floccing have shown that with few excepMons Improved cell 
performance resulted with flocced electrodes. For example data from 
Navy Power Plant tests Indicated that fuel cells wjth fleeced electrodes 
averaged Id nv higher ^^ tanea.d electrodes which translated to an 
Increase of 0.1V on a 3^pi stack. Also, after 500 hour, test at. about 
175 amperm ^ the^average 

while the standard cells lost 9.9 mv per cell on the average. 

2. Matrix Materials 

Asbestos was one of the earliest matrix materials used In fuel cells and 

is still the best’ one. . .Other materials have been tried hut are 
not used for various reasons. Some of those are discussed below: 

a) Brucite (Magnesium Hydroxide) - Not available in quantity. 

b) LItofol-S-asbestos - dissolved too rapidly In hot alkali. 

c) Silicon nitride - also dissolved In hot alkali. 

d) Polybeniimidazole (FBI) - This material had excellent resistance 
to alkaline solution but NASA-Lewls research Indicated that this 
material poisoned the catalyst. 
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e) Fybex (Potassium Titanate) (PKT) - This material had excellent 
rejststance to alkali solution and tests in fuel cells wi^^h 
Fybex matrix indicated after 14 hours that the carbonate 
production was only 50% of that in standard asbestos matrix 
cell. However this material is not available. It was re- 
ported by a phone conversation with Johns Manville rep- 
resentative that Fybex was a rapid tumor producing agent. 

Since asbestos is the best material, an improved matrix may be obtained 
by washing or leaching the asbestos matrix with alkaline solution or 
possibly EDTA before use in cell. This has resulted in a 50% 
reduction in carbonate, formation when compared to the standard 
asbestos matrix cell over 1400 hours of operation. 

3. Frame Materials 

This is definitely a problem area as the fiberglass-epoxy frames 
presently used in fuel cell corrode in hot alkali electrolyte solution. 
Of the other frame materials tried such as polyphenylene sulfide, 
polysulfone, and Celanese epoxy Novolac; the polysulfcne had a very 
low corrosion rate in KOH solution. It is not easily molded however, 
and would best be used as a hybrid with fiberglass where contact 
with the electrolyte is limited to the polysulfone material (Figure 
8). Test by UTC reported a 44% lower carbonate level in the electrolyt 
after 5000 hours in cell with hybrid polysulfone-fi berg! ass compared 
to standard epoxy-fiberglass frame cell.^^^ 

4. Graphite ERP - This change compared to Nickel ERP is proposed to 
effect a weight reduction in a 32 cell stack of 8.3 lbs or about 13% 


FIGURE 8 
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weight reduction In the fuel cell stack. Graphite was used In 
early fuel cells as electrode material but they were not long- 
life cells. The graphite ERP Is In a reducing atmosphere and 
therefore no adverse ffects on perfon. nee would be expected. 

5. Platinum on Carbon as an anode has been tested and preliminary 

studies show good long-term stability {6000 hours). The platinum 

on carbon electrode has more voltage loss initially but seems 

stable and possibly improves with time although the data are 
( 9 ) 

sketchy . ‘ 

6. The bonded cell (Figure 9) is a technique to completely encapsulat 
each cell in its reactant housing with thermo-setting plastic 

bond material . This technique provides a good seal to contain 
the electrolyte and also provides a significant weight reduction, 

16.2 lbs. per 32 cell stack which is about a 25% reduction in 
stack weight. However, the fabrication technique has produced 
some problems in the cell, and further work is needed in this area. 
Also since new matei^ials are used, these could present unforseen 
problems in cell operation. Thus a thorough screening and 
selection activity must be an integral part of a new cell 
construction approach. 

RECOMMENDATIONS 

Based on the past, present and anticipated future needs of the orbiter 
fuel cell, the following recomr^^endations are made with regard to 
future fuel cell enhancement. These recoirmendations are based on the 
presumption that for the immediate future, fuel cell operating life is 
more important than weight reduction. 



FIGURE 9 
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1. Short-Range Enhancements 

The following items can be considered for the Orbiter Fuel -Cell 
Program in the time range from immediate use to a one-year testing 
and acceptance program. 

a. Electrode Floccing - This procedure is ready for 
immediate incorporation and involves no major changes in 
materials or processes. The major effect would be a more 
uniform, consistant-performing electrode. This appears to 
be almost a 0% risk of decreasing performance contrasted to 
a 90% possibility of improving performance. Standard pro- 
duction results in a 10% rejection of electrodes whereas 
flocced electrodes have a 0% rejection so far. For "" power 
plants (the present orbiter order) containing 96 cells at 
$1500 X 96 X 15 X 10% = $ 216,000. It is difficult to oive 
a dollar value for the improved performance, but this will 
proivde a higher voltage per cell and a longer life before 
falling below the minimum voltage specification. 

b. Matrix Leaching - This procedure is ready for immediate 
incorporation into the Orbiter Fuel-Cell Program. This 
procedure provides a more uniform composition of the matrix 
by leaching out detrimental impurities that are present in 
variable amounts. It is particularly effective in removing 
calcium carbonate and this has been identified as a problem 
in the deterioration of the cell electrolyte. The leaching 
agent is a solution of ethyl enediaminetetraacetic acid (EDTA) 
and is known to complex divalent metal ions such as calcium 
and iron. Data presented by Ben Gitlow from United Technology 
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Corporation in a March 3, 1981 presentation at JSC indicated 
that a leached matrix fuel cell had only 50% of the carbonate 
build-up in the electrolyte after 1390 hours compared to a 
standard production cell. This could double the life of the 
fuel cell since carbonate build-up is the principal degradation 
mode. 

Based on a presentation by Fulton Plauche' of the Power 
Generation Branch, JSC dated June, 1981 the following data 
are used for cost estimates. The 3-substack BOM fuel cell 
is estimated to have a 4000 hour life and for four orbiter 
vehicles requiring 20,000 hours of fuel cell life each, 
the estimated cost is $ 127.6 million. Flocced electrodes 
and leached matrix enhancements have a good possibility 
(80% my estimate) of achieving a 10,000 hour fuel cell. This 
would eliminate the refurbishments and rebuilds in this 
total program as well as two fuel cell changeouts for a 
savings of $ 42.1 million or a 33% cost savings in the fuel 
cell program. 

Therefore, I strongly recommend immediate incorporation of 
these two enhancements described above. 

c. New Frame Material - The present fiberglass-epoxy frame 
is subjected to long-term degradation in the hot alkaline 
solution in the fuel cell. If a satisfactory material could 
be found, then my estimate would be that with flocced 
electrodes, leached matrix and improved frame material the 
probability of a 10,000 hour fuel cell is close to 100%. 
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The most promising frame material to be incorporated in the 
fuel cell in the near future would be a hybrid structure of 
polysulfone and fiberglass. Polysulfone has been shown to 
have a very low corrosion rate (weight loss) in alkaline 
solution and also in a 5000 hour test in fuel-cells. It 
exhibited a 50% reduction in the carbonate produced compared 
to the standard cells with fiberglass-epoxy frames. 


Therefore I recommend that effort be put forth to develop 
manufacturing procedures for this hybrid frame and concurrently 
do long term testing of this frame in short stack test cells. 

The. possible savings of $ 42.1 million; has already been 
mentioned for the Orbiter Fuel -Cell Program by developing a 
10,000 hour cell and these three short-range enhancements have 
the potential to provide that life. 


2. Long-Range Enhancements 


There are several areas of investigation that could improve fuel 
cell stability and life or reduce weight and these are recommended 
for longer range study (one - three years). 

a. A gr«iphite electrolyte reservoir plate (ERP) would 
provide a significant weight reduction (13%) over the 
standard gold-plated nickel ERP. Fabrication, assembly and 
long-term testing would be required before new materials could 
safely be used in the cells. At the present time, extended 
cell life is more important than weight reduction and thus 
this enhancement is assigned to longer term study. 
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b. The platinum-on-carbon anode has been proposed for a 
long-time fuel cell. Preliminary testing on this electrode 
has shown a good stability for 5000 hours although the initial 

(5) 

performance was about 15 mv lower than the standard cell. 

This modification involves new material introduced into the 
cell, new electrode design and long term testing and evaluation 
of this modification. 


c. Finally, the bonded cell approach has the potential to 
reduce the weight of the fuel cell {2B% reduction in stack 
weight). However, the techniques of bonding a new frame 
material and all the associated problems of leakage, 
thermal expansion and corrosion problems would require long- 
range investigations. 

The problem in evaluating these long-range enhancements is to 
determine their cost effectiveness. If it were possible to 
extend the life of the fuel cell stack to 2C,000 hours is this 
as effective as extending it from 5000 hours to 10,000 hours? 
The answer is probably "no“ since the accessory section 
components are designed for 10,000 hours and will require 
servicing at that interval. Thus it becomes uneconomical to 
overdesign one part of the fuel cell power plant compared to 
the other. 


It seems reasonable to continue some long-range improvement 
work since beyond the first four orbiter vehicles is the 
possibility of an on-going effort for other orbiter vehicles 
and even other applications unknown at the present time. 


20 


ORIGINAL PAGE B 
OF POOR QUALITY 

However, I feel that the short-range enhancements have the 
potential to produce a 10,000 hour fuel cell and a major 
program decision must be made concerning the effort to 
develop the fuel cell life beyond this. 

POTENTIAL PROBLEM AREAS 

1. One of the contaminants Identified In the degradation of fuel cells 
Is carbonate which is believed to come mainly from the matrix and 
the frame materials. There is the possibility, however of almost 
any organic material in the presence of oxygen, catalyst and alkaline 
solution to be oxidized to carbon dioxide and then converted to 
carbonates. This means that any organic solvents used in washing of 
electrodes or matrix materials must be thoroughly washed away. 
Experimental work has shown that platinum and many other catalytic 
metals bond very tightly to organic molecules so that heating above 
the normal boiling point will not remove the organic material. In 
some cases very high temperatures (400 - 500°C) and a high vacuum 
are necessary to remove contaminating materials from a catalyst surface. 
This heat conditioning in vacuum or in other gas atmospheres is 
called "activating" the catalyst. After activation the catalyst is 
handled very carefully to avoid contamination before it is used. 
Therefore the use of organic solvents or other organic additives in 
processing the fuel cells needs to be carefully evaluated to determine 
possible contributions to accelerated fuel cell degradation or to 
reduced activation of the electrodes. 
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2. Another potential problem area is the introduction of new 
materials* graphite electrodes, and graphite , . . 

ERP. Poor activation of electrodes or degradation of the cell 
performance could be due to trace amounts of materials present in 
these materials. For example sulfur-containing substances are 
notoriously bad catalyst poisons (i.e. reduce activity of catalyst) 
and polysulfone which is a good frame material contains sulfur. Other 
contaminanti- in graphite or even the products of a slow oxidation 

' could produce contaminants in the cell. Therefore, poor electrode 

i 

I activity or cell performance degradation should not be dismissed 
as just a statistical probability but careful analysis of data and 
materials used as well as chemical analysis of electrodes and 
electrolyte should be performed to identify other (if any) 
contaminating materials. 

OTHER POSSIBLE ENHANCEMENTS FOR CONSIDERATION 

In the cours^of reading and preparing this study two possible enhancements 
were identified for consideration. 

1. Electrode activation is more serious at the cathode than at the 
anode. Early fuel cell work indicated that a radioactive electrode 

( 12 ) 

lowered the activation potential resulting in increased cell voltage. 

It seems possible then^ to incorporate a weakyd* source into the oxygen 
electrode to reduce the activation potential and improve the voltage 
efficiency of the fuel cell. 

2. If graphite electrodes are reconsidered, then the use of an acid 
electrolyte is possible. The asbestos matrix and fiberglass-epoxy 
frame would be less subject to corrosion in the acid medium. Also, 
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carbonate formation would not be a problem in acidic solution so 
that carbonate-forming contaminants in reactant gases as well as from 
other materials would not be as critical. 

This report is not an exhaustive review of the fuel cell work to date 
nor are all factors considered in regard to cost savings or reliability. 
Production costs of changeover to new processes are unknown and could 
not be considered. 

It was my intent to evaluate certain proposed enhancements and make a 
judgement in my opinion as to the most promising enhancements to be 
incorporated into the orbiter fuel cell and to provide an estimate of 
the cost savings and/or reliability factor. 
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